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A REAL-TIME DIGITAL PROGRAM FOR ESTIMATING AIRCRAFT STABILITY 
AND CONTROL PARAMETERS FROM FLIGHT TEST DATA BY USING 
THE MAXIMUM LIKELIHOOD METHOD 

By Randall D. Grove and Stanley C. Mayhew* 

Langley Research Center 

SUMMARY 

A computer program (Langley program Cl 123) has been developed for estimating 
aircraft stability and control parameters from flight test data. These parameters are 
estimated by the maximum likelihood estimation procedure implemented on a real-time 
digital simulation system, which uses the Control Data 6600 computer. This system 
allows the investigator to interact with the program in order to obtain satisfactory results. 
Part of this system, the control and display capabilities, is described for this program. 
This report also describes the computer program by presenting the program variables, 
subroutines, flow charts, listings, and operational features. Program usage is demon- 
strated with a test case using pseudo or simulated flight data. 

INTRODUCTION 

A computer program has been developed at the Langley Research Center to improve 
the capabilities for estimating aircraft stability and control parameters. Improved capa- 
bilities result from the combined utilization of the inherent features of the Langley real- 
time simulation (RTS) system and the maximum likelihood method. The program has been 
designed to take advantage of the integrated software and hardware features comprising 
the RTS system. This system allows the analyst to be an integral part in controlling the 
direction of the parameter identification study, that is, the analyst communicates or inter- 
acts with the program. The analyst interacts with the program by selecting the mathe- 
matical model to be used, the variables to be matched with the flight test data, and the 
parameters to be estimated. The RTS display features aid the analyst in determining 
whether any further interaction is necessary and when to stop the computer run. The 
analyst uses the results of the computer run as an aid to solve the problems of correla- 
tion, uniqueness, . and identifiability of the parameters. Parameter estimation is not a 
straightforward procedure; thus, having the capability of the analyst interacting with the 
program is highly desirable. 

*Electronic Associates, Inc. 




The objective of this report is to present the computer program which was written 
especially for the RTS system at Langley. A description of the system software is beyond 
the scope of this report; however, the RTS subroutines used are briefly described in order 
to aid in understanding the flow of the program. The program flow differs from batch 
programs in that operator action is required. In presenting the program, it becomes 
necessary to describe the program control station (appendix A) from which the analyst 
interacts with the program. The maximum likelihood estimation procedure for extracting 
the stability and control parameters was developed in reference 1. 

The program has been dynamically checked by comparing its output with the output 
of an independently written batch program. Test cases were made by use of simulated 
flight data with different measurement noise levels to check the estimation procedure. 

(See ref. 1.) Application of this program is not restricted to the aircraft example being 
used in this report. The program may be applied to other aircraft as well as other 
dynamic vehicles satisfying the assumptions of the maximum likelihood estimator. (See 
ref. 1.) The program developed has been successfully applied to the analysis of flight 
test data for generically different aircraft. (See refs. 2 to 5.) These reports on the 
analysis of flight test data also reflect the desirability of analyst program interaction. 

SYMBOLS 


Symbols on CalComp plots (figs. 2 and 4) are not standard because of the limita- 
tions of the computer. 


A 

a I 


a X,I> a Y,I’ a Z,I 


b 


( ^^t»®a,t»®r,t 

( Cm k,t. 6 e,t 


sensitivity coefficient matrix 

accelerometer vector at instrument location 

longitudinal, lateral, and normal (positive down) components 
of aj 

wing span 

rolling-, pitching-, and yawing-moment coefficients 
rolling- moment coefficient at /3 = j^, 6 a = 6 a ^, 6 r = 

pitching- moment coefficient at ar a = 6 e = 5 e> t 
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^ n ^t’ 6 a,t> 6 r,t 

C T’ C T,0’ C T B 

Ct T’ Ct T 0 ,Ct /3 t 

C X> C Y> C Z 

c X,o> C Z,o 

( Cx )a'a j t 5 6 e ,t 

c Y,o 

< ' CZ ^a,t’ 6 e,t 

c 

D e , d t 

F 

F i 

G 

g 

Ix>IY>Iz 

J XZ 

i,j,k 

% 


yawing-moment coefficient at /3 = fy, 6 a = 6 ajt , 6 r = 6 r> t 

main- engine thrust coefficients 
tail-rotor thrust coefficients 

longitudinal-, lateral-, and normal-force coefficients 
longitudinal -force and normal-force coefficients at a a = 6 e = 0 
longitudinal -force coefficient at a a = a a> t, §e = ^e,t 

lateral-force coefficient at 0 = 6 a = 6 r = 0 
lateral-force coefficient at 0 = 0t, 5 a = 6 a 6 r = 6 r ^ 

normal-force coefficient at a a = a? a t , 5 e = 6 e ^ 

mean aerodynamic chord 

blade diameters of main engine and tail rotor 

vector function defined in equations of motion 

components of F, where j = 1, 2, . . 8 

sensitivity equation matrix 

acceleration due to gravity 

moment of inertia about X-, Y-, and Z-axis, respectively 
product of inertia 
integers 
wing tilt angle 
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performance index function 


J N 


^ r b^TP 

M X> M Y> M Z 

m 

N 

P E> P T 

p,q,r 


P' 

R 

S 

T 


Tx } T Y ,T Z 

t 


ti 

u,v,w 

V 


V SS 

X,Y,Z 


x k 


coefficients in moment equations 
rolling, pitching, and yawing moments 
mass 

number of data points 

normalized throttle settings of main engine and tail rotor 

roll, pitch, and yaw angular velocities 

number of parameters 

measurement noise covariance matrix 

wing area 

flight -test time period 

thrust along X-, Y-, and Z-axis, respectively 
time 

data point time, where i = 1, 2, . . N 

longitudinal, lateral, and vertical velocity components 

true airspeed 

slipstream airspeed 

coordinate axes 

state vector 

components of x, where k = 1, 2, . . ., 8 
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X x>yx> z x 

x Y»yY’ z Y 

x Z’yz> z z 


y 

a 

Aa 

a a.’P 

a i 

6 

6 a> 6 e> 6 r 

6B,AB 

6 ik 

V 


VE’VT 


p 


a i a \ 


,o 

a i a i 


center -of-gravity offsets of X-axis accelerometer 
center -of-gravity offsets of Y-axis accelerometer 
center -of -gravity offsets of Z-axis accelerometer 

variables in performance index function 

parameter vector 

parameter change vector 

angles of attack and sideslip 

components of a, where i = 1, 2, . . 40 

tail -rotor blade angle 

control deflection vector 

aileron, elevator, and rudder control deflections 

relationships of blade angles for main engines 

Kronecker delta 

measurement noise vector 

main- engine and tail -rotor speeds 

components of rj, where i = 1, 2, . . ., 11 

atmospheric density 

correlation coefficient of and aj 

variance and standard deviation of 
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covariance of and rj- 




4>A^ 


roll, pitch, and yaw angles 


Stability and 

control derivatives: 

* 

r - 8Cx 
2V 

ac Y 

Cy p-^ 

2V 

8 C 7 

CZq ^ 

2V 

CX “a = ^ 

ac Y 

C Yr=^r 

®2V 

ac z 

Cz «a" 

9C; 

2V 

ac Y 

%‘TF ' 

aCz 

Cz «e = s5 e 

ac i 

C ‘r-^t 

S 2 V 

ac Y 

CY rS 

Tv 

8C n 

: C "p“-3b 

3 5v 

ac / 

C VW 

ac Y 

Cy 5 r = d 6 r 

3C n 

C "r=^b 

*27 

aC; 

Cu = -4 

* a^ 
2V 

r _ 8C m 

q ■ »«* 
a 2 V 

ac n 

dC L 

c \ = Mi 

C - 3C m 

Cm «a “ 9 <*a 

ac n 

Cn r> 

a 2V 

ac, 

c \-w t 

r - aCm 

Cm 4a ' 

Tv 

ac n 

Cl V^ 


c - aCfn 

Cm 6 e ‘ M e 

ac n 

Cn h = 95 r 


6 



Superscripts: 


-1 

inverse matrix operation 

M 

measured value 

o 

nominal evaluation 

T 

Subscripts: 

transpose matrix operation 

eq 

active variables in x 

J 

active variables in y 

t 

trim condition 


Dot over a symbol denotes time derivative. Arrow over symbol denotes vector. 

|R| denotes determinant of R. 

PROBLEM DESCRIPTION 

The stability and control parameters are the unknown coefficients in the differential 
equations of motion of the aircraft. The maximum likelihood estimation procedure, using 
the method of quasilinearization, estimates the stability and control parameters by mini- 
mizing the difference between the flight test measurements and the calculated solution of 
the differential equations of motion of the aircraft. 

The measured control deflections of the aircraft are used as inputs to the equations 
of motion, and the flight test measurements are assumed to be the true solution with mea- 
surement noise (Gaussian with zero mean). The initial conditions of the state are included 
as unknown parameters and the accelerometer equations supplement the equations of 
motion in the estimation procedure. 

Aircraft Mathematical Model 

The detailed nonlinear aircraft mathematical model is described in appendix B, 
where the equations of motion are represented in vector notation by 



For simplicity in describing the estimation procedure, let the equations of motion be 
written as 



(la) 


where the terms omitted in equation (la) are auxiliary relationships. The equation vari- 
ables (states) are 


x(a,t) = [u,v,w,p,q,r,0,$] T 


The parameter vector is 


a = 



a 


2 ’ 



( 2 ) 


(3) 


The control inputs are 


6 - pa^ej^r] — ® 


(4) 


Integration of the equations of motion yields the nominal solution x(a°,t), where 
a° is the nominal or current value of the parameter vector. These parameters are the 
aerodynamic coefficients (stability and control parameters) and the initial conditions of 
the states. The initial values of the coefficients are determined from a prior estimate 
(wind-tunnel data or analysis) and the initial states are determined from the flight test 
data. 

The accelerometer equations 



(5) 


are algebraic functions of the states and their derivatives. These equations need only to 
be evaluated and not integrated. 


Sensitivity Equations 


The sensitivity equations are derived by formally differentiating the equations of 
motion with respect to each parameter. (See ref. 1.) The sensitivity equations are 
(by using eq. (la)) 


8 / \ 

d / 3x \ _ V 9F [ dx k] 9F _ dx \ 3F 

dt\ 9 «i/ Li 9 aq \ 9 a i/ 9a i 

k=l 


(i = 1, 2, . . ., p’) 


( 6 ) 
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This system of equations represents p' sets of eight simultaneous differential equa- 
tions. Integration yields the sensitivity coefficients 9x^90^. All the initial conditions 
are zero except for 


9Xk 

9oq 


t=0 


9xjj 

3Xi(0) 


t=0 


~ 6 ik 


The accelerometer sensitivity coefficients 9a I/ / 9Q! i are functions of the sensitivity 
equations and coefficients previously defined. The sensitivity equations and acceler- 
ometer sensitivity coefficients are presented in detail in appendix C. 


Maximum Likelihood Estimation Equations 


The maximum likelihood estimation equations (using quasilinearization) are derived 
from the likelihood function in reference 1. The maximum likelihood parameter estima- 
tion procedure (the accelerometer equations being neglected) is diagramed in figure 1. 
The estimation procedure is initially formulated by using the complete mathematical 
model; then, by using variable dimensioning (appendix D), it is reduced to the specific 
flight test case being analyzed. 

The parameter-change equations are 


where 



N 

v — \ 

-1 

■"n 

V — ' 

A a = 

2 A T (ti) R-l A(ti) 


2 A T (ti) R- 1 fT(ti) 


i =1 


i=! 


A(ti) 


9y° 9y° 9y° \ 

9 a i 9 « 2 ’ ’ ' *’ 9a p . / 


^ (ti) = y M (ti) - y°(ti) 


( 7 ) 



(ti) 


x M (ti) 

_af(tO 


y°(ti) = 


x°(ti) 

a?(ti) 


Here y° 


denotes the variables in the performance index function. 
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The covariance matrix for the parameters is 


£ AT( ti ) R-l A(ti) 
i=i 

The covariance matrix for the measurement noise is 

N 

R°(N) 4 Estimate of R = — ^ 77 " (ti) rf^(ti) 

i=l 

and the performance index function to be minimized is 

J N (3°) = |r°(N)| (9) 

The flight test data are composed of the onboard instrument measurements of the 
aircraft behavior and are assumed to be the output of the aircraft mathematical model 
superimposed with instrument noise. These data contain many individual aircraft maneu- 
vers stored on one magnetic tape, each maneuver being easily accessible to the central 
memory of the computer. These data are used for comparison with the mathematical 

model output and for initialization of and control input to the equations of motion. The 
— M —M 

measurements y (tj) and 6 (ti) for i = 1, 2, . . ., N are known for all performance 
index variables and control deflections corresponding to the aircraft mathematical model. 

The steps in the maximum likelihood estimation procedure, corresponding to fig- 
ure 1 , are as follows: 

(1) Initialize the system parameters where a° denotes the nominal or current 
values of the parameters. 

(2) Integrate the equations of motion and the sensitivity equations to obtain the nomi- 
nal solution and the sensitivity coefficient matrix, respectively. 

(3) Calculate the comparisons of the flight test data and nominal solution for each 
data point time at tj where i = 1, 2, . . ., N and tj = 0 and tj^ = T. 

(4) Calculate the maximum likelihood estimation equations from the comparisons in 
step (3) and the sensitivity coefficient matrix in step (2). (Dashed lines in fig. 1 indicate 
accumulation of information over the flight test time period T.) 

(5) Calculate the performance index Jj^a 0 ). 

(6) Calculate the parameter changes Aa and the statistical information matrix 
R°(N). 
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(7) Update the nominal parameter values in step (1) to start the next iteration pro- 
cedure and repeat steps until convergence. Convergence of the estimation procedure is 
assumed when the change in the performance index is small enough to satisfy the criterion 
of the analyst. . 1 . : 


PROGRAM DESCRIPTION 
Labeled COMMON 

The following list contains the FORTRAN variables appearing in labeled COMMON 
and descriptions of each variable: 


COMMON 

label 

FORTRAN 

variable(s) 

Description 

/INTCOMM/ 

T 

t, updated in subroutine IGRATE1 


H 

Integration step size used in subroutine 
I GRATE 1 (H = DT) 


INT 

! Flow control parameter used in subroutine 
IGRATE1 


NEQ 

Number of integrations performed in sub- 
routine IGRATE1, NEQ s 249 


ISCHEME 

' Selects integration scheme in subroutine 
IGRATE1 


DERINT(2,249) 

Array of integrals and derivatives in sub- 
routine IGRATE1 

/INTINTR/ 

INTERN(5,249) 

Temporary storage array for subroutine 
I GRATE 1 

/REALTIM/ 

ADC(32) 

Analog-to-digital converter input array 


DAC(64) 

Digital -to-analog converter output array 


LDISI(108) 

Logical discrete input array 
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COMMON 

label 


/REALTIM/ 


/ALGOR/ 


FORTRAN 

variable(s) 


LDISO(196) 


> 

NOPER 


NHOLD 


NRESET 


NTERM 


NPRINT 

Not used with 


> subroutine 

NREAD 

PCCEXEC 


J 

NPAR 


IPAR 


INTP(30) 


IP 


INTV(8) 


IV 


INTA(ll) 



Description 

Logical discrete output array 


/Return addresses from subroutine RTMODE 


IA 

LAI 


Total number of parameters 

Maximum number of active parameters 

Array denoting active parameters 

Number of active parameters 

Array denoting active equation variables 

Number of active equation variables 

Array denoting active performance index 
variables 

Number of active performance index variables 

Number of active equation variables which are 
active performance index variables 


IA2 


IA-IA1 
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COMMON 

label 


/ALGOR/ 


/COMM1/ 


FORTRAN 

variable(s) 

PARAM(40) 

DPARAM(30) 

ALV(ll) > 

DVAR(8) 

DALG(ll) 

ALG(40) 

I AC (40) 


IEVEN 

WT(11,11) 

COM 

LI 

L2 

IRR 

IPL 

ISL 





j. 

TABLE(199) 


INTEG(99) 


Description 


Arrays of labels for printout 


Auxiliary storage array for AL 

Auxiliary storage array for INTEG(I) 
(I = 1, 2, . . 40) 

Used with LOGIC(ll) 

rJ _1 (n) 


Labels for printout 


Error condition set by subroutine NAMECRT 


Overlay level numbers used by subroutine 
PCCEXEC 


Array of floating-point numbers to be dis- 
played and/or changed 

Array of integer numbers to be displayed 
and/or changed 
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COMMON 

label 


/COMM1/ 


FORTRAN 

variable(s) 

LOGIC (20) 

NTAB 

NINT 

NLOG 

NADC 

NDAC 

NLDI 

NLDO 

NT 

AXI 

AYI 

AZI 

DRAD 

RADD 

PI 

IR(2) 

NRN 
NPTS . 


Description 

Logical array for selecting program options 


Dimensions used by subroutines DATABLX 
and INOUT 


Recording frequency for real-time disk file 

a X,I 

a Y,I 

a Z,I 

57.2958 

(DRAD)" 1 

3.141593 

Array of integers used by subroutine GETRAN 
Number of random number sequences 
Maximum number of data points 
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COMMON 

label 


/COMM1/ 


/COMM2/ 


FORTRAN 
. variable(s) 

ISKIP 

JSKIP / 

KSKIP 

MF 

TX 

TY 

TZ 

AMX . 

AMY 

AMZ 

VARCHNG 

ITYPE 

IVARBUF(5) 

DELX 

NOITSPS 

PTSINV 

MAXPAGE 


Description 


Skips initialization statements 


Output disk file 

Tx 

Ty 


T Z 

M X 

My 

M Z 

Flag set by subroutine DSPLAY 

Data format set by subroutine DSPLAY 

Buffer used by subroutine DSPLAY 

Flight data time interval used in program 
(INC • TT) 


Integration steps per flight data time 
DELX 


interval 


DT 


(N) 


-1 


CRT memory allotment 



COMMON 

label 


/COMM2/ 


FORTRAN 

variable(s) 

LABT ^ 

LABU 

LABV 

LABW 

LABP 

LABQ 

LABR 

LABTH \ 

LABPH 

LAB AX 

LABAY 

LABAZ 

LABDA 

LABDE 

LABDR 

FET25(17) 


LUN25 

NAM25 


Description 


Axis labels for plotting 


Array containing file environment table for 
flight data file (TAPE 2 5) 

Logical unit number for TAPE25 

Name of file TAPE25 used in subroutine 
CREATE F 
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COMMON 

label 

/COMM2/ 


/FLIGHT/ 


FORTRAN 
variable (s) 

ITRNMLT(7) 

IMULTA(7) 

INVSEN(7) 

IMULTB(7) 

INVWT(7) 

UM(201) 

VM(201) 

WM(201) 

PM(201) 

QM(201) 

RM(201) 

THM(201) 

PHM(201) 

AXM(201) 

AYM(201) 

AZM{201) 

FDA(201) 

FDE(201) 

FDR(201) 


Description 


Arrays of matrix operation and dimension 
information used by subroutine MASCNT 

U M (ti) 

V M (ti) 

w M (ti) 

p M (ti) 

q M (ti) 

r M (ti) 

0 M (ti) 

^ M (ti) 

*<t0 

a Y,l(ti) 

a z,i( t i) 

«5*(ti) 

6?(ti) 
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COMMON 

label 

FORTRAN 

variable(s) 

/PLIGHT/ 

FBT(201) 


FDB(201) 


FDELB(201) 


FPPER(201) 


FTPER(201) 


BE TAT 


ETA 


ETAT 


FCT(12) 


CTT 


CT 


Description 

^T^i) 

6B(ti) 

AB(ti) 

PE(ti) 

p t(U) 

Ve 

7 7 r p 

Array containing function for /3 T 

/Jrp i 

Ct t 

c T 


Display Arrays 

The real-time simulation program uses specified arrays (subroutine DSPLAY 
arrays) for displaying and/or changing the value of desired program variables. (See 
appendix A.) The desired program variables as defined in these specified arrays are 
assigned display addresses as shown in the following table: 


Subroutine DSPLAY 
arrays 

Display 

address 

Maximum number 
of elements 

TABLE (I) 

I 

199 

INTEG(I) 

I + 200 

99 

LOGIC (I) 

I + 300 

99 

ADC (I) (not used) 

I + 400 

99 

DAC(I) 

I + 500 

99 

LDISI(I) 

I + 600 

99 

LDISO(I) 

I + 700 

199 
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This type of addressing is called "IN TABLES" addressing. For program variables not in 
subroutine DSPLAY arrays, a type of addressing called "NO TABLES" addressing is used. 

Subroutine DSPLAY arrays are listed below with their associated FORTRAN vari- 
ables and descriptions (elements not mentioned are not used). The array elements are 
equivalenced to the FORTRAN variables, except where equality signs indicate. TABLE is 
a floating-point number array, and equivalence between FORTRAN variables is indicated 
by a semicolon. 


TABLE 

element(s) 


FORTRAN 

variable(s) 


Description 


1 

AL(l);UO 

u(0) 

2 

AL(2);CXO 

( c x)« 

3 

AL(3);CXAL 

CX -a 

4 

AL(4);CXQ 

C Xq 

5 

AL(5);THEO 

0(0) 

6 

AL(6);PHIO 

0(0) 

7 

AL(7);WO 

w(0) 

8 

AL(8);CZO 

' ( c z)« ; 

9 

AL(9);CZAL 

C Z a 

“a 

10 

AL(J.O); CZQ 

C Z q 

11 

AL(11);CZDE 

CZ «e 

12 

AL(12);QO 

q(0) 

13 

AL(13);CMO 

( Cm )o! 


a a,t’ 6 e,t 
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TABLE 

element(s) 

fortran 

variable(s) 

Description 

14 

AL(14);CMAL 

“a 

15 

AL(15);CMALD 

Cm,v 

“a 

16 

AL(16);CMQ 

Cm q 

17 

AL(17);CMDE 

Cm 6 

°e 

18 

AL(18);VO 

v(0) 

19 

AL(19);CYO 

( C y)o 6 
Pt> 6 a,t ; 

20 

AL(20);CYB 

C Yj3 

21 

AL(21);CYBD 

CY ^ 

22 

AL(22);CYP 

C ^P 

23 

AL(23);CYR 

Cy r 

24 

AL(24);CYDR 

c y 5 

°r 

25 

AL(25);PO 

p(o) 

26 

AL(26);CLO 

(Cl) 

[ Va,t 

27 

AL(27);CLB 

c lp 

28 

AL(28);CLBD 

Cj- 
6 (3 

29 

AL(29);CLP 

Cl 

‘'P 

30 

AL(30);CLR 

C *r 
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TABLE 

element(s) 

FORTRAN 
variable (s) 


Description 

31 

AL(31);CLDR 

c h 

°r 


32 

AL(32);CLDA 

C lb 

°a 


33 

AL(33);RO 

r(0) 


34 

AL(34);CNO 



35 

AL(35);CNB 

C ' 1 ;3 


36 

AL(36);CNBD 

S 


37 

AL(37);CNP 

C "P 


38 

AL(38);CNR 

c n r 


39 

AL(39);CNDR 

Cn 6 r 


40 

AL(40);CNDA 

C n 6a 


50 

UMAX 



51 

VMAX 

! 


52 

WMAX 



53 

PMAX 



54 

QMAX 

> CalComp, CRT, 

and DAC scaling information 

55 

RMAX 



56 

THMAX 



57 

PHMAX 



58 

AXMAX 

y 
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TABLE 

element(s) 

59 

60 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 . 

75 

76 

77 

78 


FORTRAN 
• variable(s) 

AYMAX 

AZMAX 

DAMAX 

DEMAX 

DRMAX 

UMULT 

VMULT 

WMULT 

PMULT 

QMULT 

RMULT 

THMULT 

PHMULT 

AXMULT 

AYMULT 

AZMIJLT 

DAMULT 

DEMULT 

DRMULT 

/ 








J 


Description 


CalComp, CRT, arid DAC scaling information 


Flight test data multipliers 


22 



TABLE 

element(s) 

FORTRAN 

variable(s) 

Description 

79 

UBIAS 


80 

VBIAS 


81 

WBIAS 


82 

PBIAS 

- 

83 

QBIAS 


84 

RBIAS 

> Flight test data biases 

85 

THBIAS 


86 

PHBIAS 


87 

AXBIAS 


88 

AYBIAS 


89 

AZBIAS 


90 

DABIAS 

5 a,t 

91 

DEBIAS 

^e,t 

92 

DRBIAS 

5 r,t 

98 

RUN 

Run number 

99 

PASS 

Iteration number 



N 

r - ^ rn 

100 

AJP 

(t-i) ^J^i) 
i=l 

101-104 

DRM(4) 

Array of desired random means 



TABLE 

element(s) 

FORTRAN 

variable(s) 

Description 

117 

ADC 

ix 

118 

AIY 

Iy 

119 

AIZ 

Iz 

120 

ABCZ 

!XZ 

121 

WEIGHT 

Weight of aircraft (mg) 

122 

GRAV 

g 

123 

RHO 

P 

124 

S 

S 

125 

B 

b 

126 

CBAR 

c 

127 

DEAMPL 

Amplitude of 6 e for test case 

128 

DEFREQ 

Frequency of 6 e for test case 

129 

ALPHAT 

a a,t 

130 

DT 

Integration step size (H = DT) 

131 

TT 

Time interval for flight test data tape 

132 

TS 

Tape starting time for putting in flight data 

133 

TIMF 

Final time for CRT display 

134 

=T 

t 

135 

=U 

u 
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TABLE 

element(s) 

FORTRAN 

variable(s) 

Description 

136 

=V 

V 

137 

=w 

w 

138 

=p 

P 

139 

=Q 

q 

140 

=R 

r 

141 

=THE 

e 

142 

=PHI 

<P 

143 

=PSI 

* 

144 

=UDOT 

u 

145 

=VDOT 

V 

146 

=WDOT 

w 

147 

=PDOT 

p 

148 

=QDOT 

q 

149 

=RDOT 

r 

150 

=THEDOT 

e 

151 

=PHIDOT 


152 

=PSIDOT 

* 

154 

! XX 

X x 


i 
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TABLE 

element(s) 

FORTRAN 

variable(s) 

Description 

155 

YX 

y x 

156 

zx 

z x 

157 

XY 

X Y 

158 

YY 

yY 

159 

ZY 

Z Y 

160 

XZ 

X Z 

161 

YZ 

y z 

162 

ZZ 

Z z 

169 

CMCON 

Constant used with LOGIC (6) 

170 

DALMLT 

Parameter step size multiplier 

171 

UCRTBI 

u bias for CalComp, CRT, and DA C 



presentations 

172 

AIW 

% 

173 

CTO 

C T,o 

174 

CTTO 

c Tt 

A o 

175 

CTB 

c t b 

176 

CTBT 

c % 

177 

D 

d e 

178 

CAPDT 

Dfp 
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table 

element(s) 

FORTRAN 

variable(s) 

Description 

179 

ELBAR 

i 

180 

-ELTP 

^TP 

181 

RB 

r b 

182 

PPER 

P E 

183 

TPER 

P T 

184 

BTBIAS ^ 


185 

PPERBI 

> ' Flight test data biases 

186 

TPERBI 



187 

DETl 

. Inverse determinant of parameter covariance 
matrix 

188 

DET2 

|rj<n)| 

191-194 

DRSD(4) 

Array of desired standard deviations of 
random numbers 

INTEG is an integer number array. 


INTEG 

element(s) 

FORTRAN 

variable(s) 

Description 

1-40 


Array denoting activeness of parameters 

51-61 

INTY(ll) 

Array denoting activeness of performance 
index variables 

62 

NOPTS 

N 

63 

INC 

Sample rate for flight data tape 
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INTEG 

element(s) 

FORTRAN 
variable (s) 

Description 

64 

IPRINT 

Selects options in overlay level (2,0) 

65 

NPLOT 

Density of plotting symbols for CalComp 
plots 

66 

IREAD 

Selects options in overlay level (3,0) 

67 

KSCAN 

Scan rate, used by subroutine SCANNER 

81-88 

INTX(8) 

Array denoting activeness of equation 
variables 


LOGIC is a logical array and the descriptions are for a true (.T.) value. 


LOGIC 

element(s) 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


Description 

Calculate Rj (N) s 

Diagonalize Rj(N) 

Set P E = FPPER(l) 

Set P T = FTPER(l) 

Calculate C-p 0 T from FCT(12) and /3 T 
/3rp 

Set Czg = CMCON • C^g 

Set longitudinal states equal to flight test data 

Set o a> t= tan-1 S® 

Calculate parameters to trim x to zero 


Automate LOGIC (9) 
Automate two -pass updating 
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Y 

1 

DAC is an output array for the time history recorder. 


DAC 

element(s) 

FORTRAN 

variable(s) 

Description 

1 

=(U-UCRTBI)/UMAX > 

!! 

2 

=V/VMAX 


3 

=W/WMAX 


4 

=P/PMAX 

| 

> Normalized time history recordings 

5 

=Q/QMAX 


6 

=R/RMAX 


7 

=THE/THMAX 


8 

=PHI/PHMAX 

J 


LDISI is a logical discrete input array where the descriptions are for true (.T.) 

values of the switches. 


LDISI 

element(s) 

FORTRAN 
variable (s) 

Description 

1-16 


Data entry keyboard 

17 


OPERATE (OPER) mode switch 

18 


HOLD mode switch 

19 


RESET mode switch 

20 


TERMINATE (TERM) mode switch 

21 


CHANGE mode switch 

22 


SCAN mode switch 
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/ 


LDISI 

element(s) 


FORTRAN 

variable(s) 


Description 



23 


RELEASE mode switch, releases CHANGE or 
SCAN mode switch 

27 


ERASE mode switch 

29 


IDLE mode switch 

30 


PRINT mode switch 

31 


READ mode switch 

32 


RELEASE mode switch, releases ERASE, 
IDLE, PRINT, or READ mode switch 

33 

FSS(l) 

Accelerometer variables for CRT display and 
CalComp plot 

34 

FSS(2) 

CalComp plot of flight data only 

35 

FSS(3) 

Fill flight data arrays with pseudo data 

36 

FSS(4) 

Initialize variable dimensioning 

37 

FSS(5) 

Iteration printout on MF file 

39 

FSS(7) 

Lateral variables for CRT display and 
CalComp plot 

40 

FSS(8) 

Control variables for CRT display and 
CalComp plot 

41 

FSS(9) 

Pseudo flight controls 

42 

FSS(10) 

Skip update of a 0 

43 

FSS(ll) 

Activates typewriter for subroutine DS PLAY 
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LDISI 

FORTRAN 


element(s) 

variable(s) 

Description 

44 

FSS(12) 

Retain x(0) in ALG array 

45 

FSS(13) 

Exit CRT or CalComp loop 

46 

FSS(14) 

Replot the CRT or CalComp plot 

47 

FSS(15) 

Enter CRT loop 

48 

FSS(16) 

"IN TABLES" addressing (false (.F.) for 
"NO TABLES" addressing) 

LDISO is 

a logical discrete output array used to turn the white indicator lights (WL) 

on (.T.) and off (.F.). The diagnostics are 

given for .T. value. 

LDISO 

FORTRAN 


element 

variable(s) 

Diagnostic 

61 

WL(1) 

Default value used for IPRINT 

62 

WL(2) 

CalComp plot completed 

63 

WL(3) 

Error in CalComp plot loop 

66 

WL(6) 

Default value used for IREAD 

67 

WL(7) 

Attempted to read flight data beyond end of 



file 

70 

WL(10) 

N > NPTS 

71 

WL(ll) 

IP > IPAR 

72 

WL(12) 

Inactive equation variable in performance 
index 

73 

WL(13) 

|0| > 1.5 , 


LDISO 

element(s) 

FORTRAN 

variable(s) 

Diagnostic 

74 

WL(14) 

u < 5 

75 

WL(15) 

|v| > u 

76 

WL(16) 

|Qf a | > 1-5 

77 

WL(17) 

Invalid flight data location used 

78 

WL(18) 

DET1 = 0 

79 

WL(19) 

DET2 = 0 

80 

WL(20) 

In CRT loop 


Subroutine Descriptions 

The following are the subroutines used and brief statements describing their 
functions. 


Subroutine 

ASCALE 

ATERM 

AXES 

CALPLT 

CLRPLOT 

CLRTABL 

CREATE F 

CRTPLOT 


Description 

Determines origin and scale factor of time axis for CalComp plots 

Does final processing and halts execution 

Draws and annotates axes for CalComp plots 

Moves plotter pen to new location or signals end of CalComp job 

Clears the CRT plot of the calculated variables 

Clears existing plot parameter tables for CRT variables 

Creates disk file for flight data tape 

Establishes plot parameters for CRT variables and/or generates 
annotated plotting grids 
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Subroutine 


Description 


CYCLE 

DATABLX 

DAYTIM 

DSPLAY 

ENDPLOT 

ERASE 

GETRAN 

GRID 

HALT 

I GRATE 1 

INOUT 

LDRSEC 

LEROY 

LINE 

MASCNT 

NAMECRT 

NM218 

NOTATE 


Sets up return address for subroutine RECYCLE 
Specifies the variable arrays for subroutine DSPLAY 
Provides date and time of day 

Activates data entry keyboard and digital decimal display unit 

Marks an end to the CRT plot 

Erases data on real-time disk file 

Generates Gaussian random numbers 

Draws grid for CalComp plots 

Signals completion of real-time portion of program 

Integrates variables in DERINT(2,J) and stores results in 
DERINT(1,J) (J = 1, 2, . . NEQ) 

Sets up arrays for input/output conversion 

Provides for PRINT and READ mode entries into subroutine 
RTMODE 

Controls CalComp plotting with liquid ink pen 

CalComp routine to draw a continuous line and/or a symbol 

Performs matrix algebra operations 

Identifies and initializes usage of the CRT console 

Initializes usage of typewriter 

Draws alphanumeric characters for CalComp plots 



Subroutine 


Description 


NUMBER 

OPERATE 

PCCEXEC 

PLAYBAK 

PRINTER 

PSEUDO 

READOUT 

READY 

RECIN 

RECORD 

RECYCLE 

RITE CRT 

RTMODE 

SCANNER 

TYPE VAR 

UNLODE 

XDSPLAY 


Draws floating point numbers for CalComp plots — - - ^ ... 

Causes readout of DAC and LDISO arrays and readin of LDISI 
array 

Controls overlay loading for RTS jobs 

Plays back data recorded on real-time disk file, one frame at a 
time 

Causes routing of MF disk file to printer 

Saves plotting information in order to use CalComp postprocessor 

Specifies quantities to be recorded and frequency of recording for 
the real-time disk file 

Signals that the program is ready for real-time operation 

Inputs flight test data 

Records quantities as specified in subroutine READOUT 

Signals end of a cycle when in the OPERATE mode and returns to 
address specified by subroutine CYCLE 

Issues plotting requests to CRT system 

Entry into the mode control routine 

Increments display address 

Types out data displayed on DDDU 

Erases data stored on CRT file 

Initializes data entry keyboard and DDDU, and routes program 
listing 
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OVERLAY Program Descriptions, Flow Charts, and Listings 

The program uses RTS central memory overlay capabilities controlled by subrou- 
tine PCCEXEC where level (4,0) is the real-time operational level. A brief description, 
flow chart, and FORTRAN source listing are presented for each overlay level (excluding 
subroutines previously described). 

OVERLAY(XC142FL,0,0) .- The main overlay, level (0,0), is always resident in cen- 
tral memory. It includes the labeled COMMON, the initial call to subroutine PCCEXEC, 
and other system initialization calls. 
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Ov£RlAY(XC142FL.O.O > A 0001 

PROGRAM XC142A( INPUT .OUTPUT) A-0002- 

LOGICAL LOISI .LDISO. LOGIC. VARCHNG a 0003 

COMMON /InTCOMM/ T.H. INT.NEO. I SCHEME. DERI NT (2. 249) A OOOA 

COMMON /INTINTR/ INTERN ( 5 .249 ) A 0005 

COMMON /REALTIM/ ADC ( 32 > i DAC (64 > .Lei S l ( 1 08 > ,LD 1 SO ( 1 96 ) ,NOPER .NHOLD A 0006 
1 .NRESET.NTERM.NPRINT.NREAD A 0007 

COMMON /ALGOR/ NPAR » I PAR . INTP ( 30 > . IP . I NTV (8 > . I V. INTA ( 1 1 ) . I A , l A1 . I A A 0008 
I2» PARAM (40) .DPARAM { 30 ) .ALV( 1 1 ) »DVAR(8)» DALG ( 1 I ) .ALG (40) . I AC (40) . IE A 0009 
2VEN.WT( 11 . n ) .COM, n ,L2 A 0010 

COMMON /COMM 1 / I RR « I PL « I SL « T ABLE ( I 99 ) . I NTEG ( 99 ) . LOG I C ( 20 ) » NT AB . N I N A 0011 
1 T. NLOG.NAOC.NO AC. NLOl ,NLOO. NT, AX 1 , AY I . AZ l ,ORAD .RAOD.P l . IR(2) .NRN.N A 0012 
2PTS. I SK IP, JSK IP, K SKIP, MF.TX.TY.TZ. AMX. AMY, AMZ, VARCHNG, I TYPE. I VAR0U A 0013 
3F(5) ,DELX,NOITSPS.PTSINV A 0014 

COMMON /C0MM2/ MAxPAGE » LAST , LABU.LAB V»L ABW.LABP »LABQ»LA8R .LABTH .LA A 0015 
1BPH,lABAX,LABAY.LABAZ,LABDA,LABDE,LABDR,FET25( 17) , LUN25, NAM25 . I TRN A 0016 
2MLT ( 7 ) » I MULTA ( 7 ) » I NVSEN ( 7 ) , I MULTB ( 7 ) , I NVWT ( 7 ) A 0017 

COMMON /FLIGHT/ UM(201 ) ,VM(201 ),WM(201 ) ,PM(201 ) ,QM(20I ) ,RM(201) ,TH A 0018 
1M (201 ) ,PHM(20l ) , AXM ( 20 1 ) , A YM { 20 1 ) , AZM (201 ) » FDA ( 20 I ) .FOE (201 ) ,F0R(2 A 0019 
201 ) ,FBT <201 ) ,FDB( 201 ) ,F0ELB(201 ) .FPPER <20l ) ,FTPER(201 ) .BETAT »ET A,E A 0020 
3TAT ,FCT (I2),CTT,CT A 0021 

CALL PSEUDO A 0022 

CALL NAMECRT (6LCRTTPE. IRR > A 0023 

XC 1 42FL=7LXC 1 42FL A 0024 

CALL PCCEXEC ( XC 1 42FL . MF , 90034S . LD I S I . I PL . I SL ) A 0025 

CALL DSPLAY A 0026 

CALL CLRTABL A 0027 

90034 FORMAT ( 54H XC-142A Cl 123 1 3927T S. C. MAYHEW ROOM 21310 A 0028 

END A 0029 
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OVERLAY (XCi42FL,l > 0) .- The initialization overlay, level (1,0), is automatically 
loaded by the initial call to subroutine PCCEXEC. Upon completion of level (1,0), 
overlay (4,0) is automatically loaded. 
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OVERLAY (XC 1 42FL « 1 .0) B 0001 

PROGRAM I NIT B 0002 

LOG I CAL LD 1 S I »LDI SO .LOGIC • VARCHNG B 0003 

DIMENSION AL ( 40 ) « 1NTX(8). INTY(ll) B 0004 

COMMON /INTCOMM/ T . H « I NT « NEQ . I SCHEME . DER I NT ( 2 • 249 ) B 0005 

COMMON /INTINTR/ I NTERN (5*249) B 0006 

COMMON /REALTlM/ ADC(32> »DAC(64 ) .LD1 SI ( 108 ) *LD1S0( 196 ) .NOPER.NHOLD B 0007 
1 .NRESET.NTERM.NPRINT.NREAD B OOOB 

common /Algor/ NPAR. IPAR. !NTP(30 > « IP. INT.V(8> • I V. INTA{ 1 1 ) , IA. I Al . IA B 0009 
1 2.PARAM (40 ) .DPARAM ( 30 ) .ALV< 1 1 ) «DVAR ( 8) .DALG (11) «ALG(40 > * I AC (40) . IE B 0010 
2VEN.WT* 1 I , | I ) .COM. LI .L2 B 0011 

COMMON /COMMl/ IRR. IPL. I SL. TABLE ( 1 99 ) . lNTEG(99 ) .LOGIC (20 ) .NTAB.NIN B 0012 
1T.NL0G.NADC.NDAC.NLDI .NLDO.NT.AXI .AY I .AZI .DRAD.RADD.PI . IR(2) .NRN.N B 0013 
2PTS. ISKIP, JSKIP.KSKIP. ME. TX . TY . T Z . AMX . AMY. AMZ . VARCHNG. I TYPE . IVARBU B 0014 
3F(5) .OELX.NOITSPS.PTSINV B 0015 

COMMON /C0MM2/ MAXPAG£.LABT.LABU.LABV.LABW.LA8P,LABQ.LABR.LABTH.LA B 0016 
1BPH.lABAx,LABAY.LABAZ.LA3DA,LABDE.LABDR.FET25( 17) .LUN25.NAM25. I TRN B 0017 
2MLT (7 ) . IMULTA (7) . INVSEN (7) . I MULTB ( 7 > . INVWT (7 ) 8 001 B 

COMMON /FLIGHT/ UM ( 20 1 ) . VM ( 20 1 > • WM ( 20 1 ) . PM ( 2 0 1 ) . QM ( 20 1 ) . RM ( 20 1 ) • TH B 0019 
1M (201 ) . PHM ( 20 1 ) , AXM (201 ) . AYM (201 ) . AZM (201 ) .FDA (201 ) *FDE(201 ) »FDR(2 B 0020 
201 > < F8T ( 20 1 ). FOB (201 >.FDELB(20! ).FPPE R(20l ).FTPER(20l ) , BETAT .ETA ,E B 0021 
3TAT « fCT (12)«CTT«CT B 0022 

DATA PARAM/2HUO . 3HCX0 . 4HCX AL . 3HCX0 » 4HTHEO » 4HPH 10. 2HW0. 3HCZO .4HCZAL B 0023 

1 .3HCZQ.4HCZDE.2HQ0.3HCM0.4HCMAL.5HCMALD.3HCMQ.4HCMDE.2HV0.3HCY0.3H 8 0024 
2C YB . 4HC YBD . 3HCYP . 3HCYR . 4HC YDR, 2HP0 ♦ 3HCL0 . 3HCLB .4HCLBD, 3HCLP . 3HCLR. B 0025 
34HCLDR . 4HCLDA . 2HR0 . 3HCN0 . 3HCNB . 4HCNBD . 3HCNP . 3HCNR . 4HCNDR . 4HCNDA/ B 0026 

DATA ALV/1HU. 1HV. 1HW. 1HP. 1 HQ . 1 HR . 3HT HE . 3HPH 1 . 3HAX I . 3HAY I . 3HAZ I / B 0027 
DATA FCT/O. . .008, .025. .046. .072, .097, . 125, .151 .. 174. ,193. .208. .22/ B 0028 
EQUIVALENCE (DERINT( 1 . 1 ) .U) . (OERINT ( 1 .2 ) ,V> . ( DER I NT (1,3). W). (DEB 0029 

1RINT ( 1 ,4 ) ,P ) . ( DER I NT ( 1 « 5 ) . Q ) « ( DER I NT ( l , 6 ) . R ) • ( DERI NT ( 1 .7 ) • THE ) • B 0030 

2 (DERINT( 1 ,8) .PHI ) , (AL(l).UO), (AL(2).CX0). ( AL ( 3 ) »CXAL ) . (ALC7). B 0031 

3WO ) . ( AL (8 ) , CZO ) , ( AL (9 ) .CZAL) « (AL(12).00), ( AL (14). CMAL ) . (AL(16B 0032 

4 ) , CMQ ) , (AL(18). VO ) , ( AL ( 20 ) , CYB ) . ,( AL ( 22 ) . C YP ) . ( AL ( 23 ) . CYR ) , (AL B 0033 

5(25). PO), ( AL ( 27 ) ,CLS ) , ( AL ( 29 ) , CLP ) » < AL ( 30 ) ♦ CLR ) . (AL(33).R0). ( B 0034 

6AL (35 ) . CNB ) , (AL (37) ,CNP) . ( AL ( 38 ) .CNR ) , ( TABLE ( 1 ). Al ( 1 )) . (TABLEf B 0035 

750). UMAX). (TABLE (51 ) »VMAX ) , ( TABLE ( 52 ). WMAX ) , ( TABLE ( 53 ) .PMAX ) . * B 0036 

STABLE ( 54 ) ,QMAX) . ( TABLE ( 55 ). RMAX >, ( TABLE ( 56 >. THMAX ) . (TABLE (57 ) .P B 0037 

9HMAX), (TABLE (58 ) .AXMAX > , ( T ABLE ( 59 > . A YM AX ) . ( TABLE ( 60 ). AZMAX ) . (T B 0038 

SABLE (62 ) .DAMAX ) , (TABLE(63) .DEMAX) . ( TABLE ( 64 ) .DRMAX ) , (TABLE(99)» B 0039 

SPASS), (TABLE ( I 17 ) . AIX) . ( TABLE ( 1 1 B ) . A I Y ) . ( TABLE ( 1 1 9 ) . A I Z ) . (TA8L B 0040 

*E( 121 )'. WEIGHT ) , (TABLE ( 122 ) ,GRAV ).( TABLE ( 123 ) .RHO ) . (TABLE ( I 24 ) .S B 0041 
S), (TABLE ( 1 25 ) .8 ) » ( TABLE ( 1 26 ), CBAR ) . (TABLE ( 1 28 )♦ DEFREQ ) » (TABLE* B 0042 

*130), DT)'. (TABLE ( 131 ), TT ) , ( TABLE ( 1 33 > »T IMF ) . ( TABLE ( 1 54 ) .XX ) . (TAB 0043 
SBLE ( 1 55 ) »YX ) . ( T ABLE ( 1 56 ) »ZX ) » ( TABLE ( 1 57 ), XY ) ♦ ( TABLE ( I 58 ) . YY) . ( B 0044 
STABLE* 159) ,ZY) , (TABLE ( 160 ) . XZ ) . ( T ABLE ( I 6 1 > . YZ ) . ( TABLE ( 162 ) .ZZ ) B 0045 

EQUIVALENCE (TABLE( 170 ) .DALMLT) , ( TABLE ( 1 7 1 ) .UCRTB I ) . (TABLE ( 173) . B 0046 
1 CTO ) . (TABLE* 177) ,D) . ( I NTEG (5 1 > « I NT Y ( 1 > > . ( I NTEG ( 62 > .NOPTS ) . ( INT B 0047 
2EG(63). INC) » ( INTEG(64 > . IPRINT) » ( I NTEG ( 65 ) .NPLOT ) . ( 1NTEG (66 ) » IRE B 0048 


3AD ) , ( INTEG (67) .KSCAN) , ( I NTEG ( 8 1 ) . I NTX ( 1 > ) 

B 

0049 

ISCHpME*3 

B 

0050 

NPAR=40 

B 

0051 

I PAR=30 

B 

0052 

I EVEN*- 1 

B 

0053 

COM* 1 H. 

C 

0054 

L 1 =8h LOGIC* 

C 

0055 

L2*2h>* 

C 

0056 

I PL®4 

B 

0057 

I SL = 0 

B 

0058 

NT AB* 1 99 

B 

0059 

DO 10 1*1 .NT AB 

B 

0060 
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10 TABLE! I )=0. 

DO 20 1=65.78 
20 TABLE! I 1 = 1. 

N 1 NT =99 

DO 30 I =1 «NI NT 
30 I NTEG ! I ) =0 

NLOG=20 

DO 40 1=1 .NLOG 
40 LOGIC!I)=.E. 

LOGIC! 1 )=.T. 

LOG I C < 2 ) = .T.» 

NADC=0 
NDAC=8 
NLDI=48 
DO 50 1=1 . NLDI 
50 LD 1ST ! I ) = »F « 

NLDO* 1 80 
DO 60 I=l.NLOO 
60 LDISO! I )=.EALSE. 

CALL DATABLX (TABLE .NT AB . I NTEG. N 
1DISI .NLDI .LDISO.NLDO) 

CALL INOUT IADC.NADC.DAC.NDAC.LD 
CALL NM218 ! 5LOSCAR ) 

NT = l 

CALL READOUT !6 .NT « U« V. W.P .0 .R > 

CALL READOUT ( 6 .NT . T . THE. PHI . AX I 

DR AD =57. 2958 

RADD=1 •/'DRAD 

Pt=3. 141593 

IR!1 )*13 

I R ! 2 ) *2357 

NRN=4 

NPTS-201 

ISKIP*JSKIP*KSK1P=0 
TV = 0. 

MAxPAGE*10 
LABTxIOHTIME I SEC) 

LABU= 1 OHU (FT/SEC) 

, LABV*10HV (FT/SEC) 

LABW=10HW (FT/SEC) 

LABP= 1 OHP ( RAD/SEC ) 

LABQsl OHO (RAD/SEC ) 

LABRcl OHR (RAD/SEC ) 

LABTH*lOHTHETA (RAD) 

LABPH*10HPHI (RAD) 

LABAx= 1 OHAXI (G) 

LABAY=1 OHAY I (G) 

LABAZMOHAZI (G) 

LABDA* 1 OHOA (RAD) 

LABDE-10HDE (RAD) 

LABOR* 10HDR (RAD) 

NAM25*6LT APE25 
DO 70 1*1.7 
I TRNMLT ( I ) =0 
IMULTA ( I ) =0 
I NULTB ( I ) *0 
INVSEN1 I ) *0 


B 0061 
B 0062 
B 0063 
B 0064 
B 0065 
B 0066 
B 0067 
B 0068 
B 0069 
B 0070 
B 0071 
B 0072 
B 0073 
B 0074 
B 0075 
B 0076 
B 0077 
B 0070 
B 0079 

NT .LOGIC .NLOG. ADC. NADC.DAC.NDAC.L B 0080 

B 008 1 

SI .NLDI .LDISO.NLDO) B 0082 

8 0083 
B 0084 
B 0085 

AYI.AZI) B 0086 

B 0087 
B 0088 
B 0089 
8 0090 
B 0091 
B 0092 
B 0093 
B 0094 
B 0095 
B 0096 
B 0097 
B 0098 
B 0099 
B 0100 
B 0101 
B 0102 
B 0103 
B 0104 
B 0105 
B 0106 
B 0107 
B 0108 
B 0109 
B 0110 
B 01 1 1 
B 0112 
B 0113 
B 01 14 
B 0115 
B 0116 
B 0117 
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70 


I NVWT ( I ) *0 
I TRNMLT < 1 ) = 23 

I TRNMLT <5 > = I TRNMLT (6 >* I MULT A(6)*l NVWT ( 5 ) * 1 1 

I TRNMLT C 7 > * I MULTA < 5 ) = I MULT A < 7 > = I NVSEN (5 ) = I MULTB ( 5 > = I MULTB ( 6 ) = IP AR 
I MULTB ( 7 ) = I PAR 
I MULTA < 1 ) a I MULTB < 1 ) <*20 
IMULTA f 4 ) * I MULTB ( 4 ) * 1 
INVSFN { 1 )*INVWT(1 ) = 10 
I NVSEN ( 4 ) a I NVWT < 4 1*2 
UMAX*20. 

VMAX=40 « 

WMAX«10. 

PMAX*OMAX=RMAX».2 
THMAX*PHMAXbAXMAX=AYMAX*«4 
AZMAX=2. 

OAMAX»DEMAX*DRMAX = . 1 
PASS«-1. 

A I X* 1 30000, 

A I Y*= 1 28000. 

A 12=270000* 

WEIGHT»36050. 

GRAV=32« 1 7 
RHO=. 002242 
S-534. 

B=67«5 
CBAR*6.07 
DEFRE0*1» 

DT=. 03125 
TT». 03125 
TIMFal 0. 

XX»10. 

YX=2. 

ZX*5. 

XY=7. 

YY*9. 

ZY*4, 

XZ=2. 

YZ»5. 

ZZ*7. 

0ALMLT«1 . 

UCRT0I =200* 

CTOb.03 
0=15.5 
00 80 1=1 »8 
INTX(I)=1 


B 0118 
B 0119 
B 0120 
B 0121 
B 0122 
B 0123 
B 0124 
B 0125 
B 0126 
B 0127 
B 0128 
B 0129 
B 0130 
B 0131 
B 0132 
B 0133 
B 0134 
B 0135 
B 0136 
B 0137 
B 0138 
B 0139 
B 0140 
8 0141 
B 0142 
B 0143 
B 0144 
B 0145 
B 0146 
B 0147 
B 0148 
B 0149 
B 0150 
B 0151 
B 0152 
B 0153 
B 0154 
B 0155 
B 0156 
B 0157 
B 0158 
B 0159 
B 0160 
B 0161 
B 0162 
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INTY { I 1 = l 

B 

0163 

NOPTS«ei 

B 

0164 

INC-4 

B 

0165 

I PRINT- 1 

B 

0166 

NPLOT-2 

B 

0167 

I READ* 1 

B 

0168 

KSCAN-2 

B 

0169 

ETA- 1232. *.75/60. 

B 

0170 

U0-120.*I .689 

B 

0171 

TMP»s*U0**2 

B 

0172 

CXO— 8.*ETA**2*D**4*CT0/TMP 

B 

0173 

CXA|_=-0.29 

B 

0174 

WO— 5 . 

B 

0175 

CZO-— 2 . *WE I GHT / (RHO*TMP ) 

B 

0176 

C2AL--4.30 

B 

0177 

Q0--.15 

B 

0178 

CMAL--1.49 

B 

0179 

CMQ= —3 1.14 

B 

0180 

VO-10. 

B 

0181 

CYB--1 .70 

B 

0182 

CYP-0.50 

B 

0183 

CYR-0.40 

B 

0184 

P0--.15 

B 

0185 

CLB— 0.18 

B 

0186 

CLP-— 0.72 

B 

0187 

CLR-0.20 

B 

0188 

R0-.15 

B 

0189 

CNB-0.06 

B 

0190 

CNP-O.05 

B 

0191 

CNR— 0.70 

B 

0192 

I NTEG < 2 > = INTEG (8 ) = I NTEG (9) - I NTEG C 1 4 > - I NTEG ( 1 6 ) « I NT EG ( 20 ) » 1 

B 

0193 

I NTEG < 22 ) - I NTEG ( 23 ) = I NTEG ( 27 > - I NTEG ( 29 ) « I NTEG < 30 ) - I NTEG ( 35 ) » 1 

B 

0194 

I NTEG < 37 ) = I NTEG ( 38 1 = 1 

B 

0195 

DO 90 1=1 » NPAR 

B 

0196 

ALG< I ) =AL ( I ) 

B 

0197 

1 AC ( I ) = INTEG < I ) 

B 

0198 

ALG ( 2 ) = — » 1 

B 

0199 

ALG (8 ) =— 1 • 

B 

0200 

ALGC9J--3. 

B 

0201 

ALG ( 1 4 ) *— 1 . 

B 

0202 

ALG (16)— 22 

B 

0203 

ALG (20 ) =0 . 

B 

0204 

ALG (22 ) - .3 

B 

0205 

ALG ( 23 ) = .2 

B 

0206 

ALG127 ) . 1 3 

B 

0207 

ALG < 29 ) -—.4 

B 

0208 

ALG (30) -.15 

B 

0209 

ALG ( 35 ) = .03 

B 

0210 

ALG ( 37 ) = .02 

B 

021 1 

ALG ( 38 ) -— .4 

B 

0212 

RETURN 

B 

0213 

END 

B 

0214 


41 



OVERLAY (XC142FL,2,0) .- The print overlay, level (2,0), is loaded when the PRINT 
mode is selected on the program control console. The integer IPRINT, preset in level 
(4,6), selects one of the five options in the PRINT overlay. The CalComp plot option is 
presented in appendix E. The primary overlay, level (4,0), is automatically loaded upon 
completion of level (2,0). 
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OVERlAY(XC1 42FL.2.0) C 

PROGRAM XCPRINT C 

LOGICAL FSS< 16) .LDI SI .LDISO.LOGI C. VARCHnG.WL (39) C 

DIMENSION X(203>. Y1 (203) . Y2(203>. Y3(203). Y4(203)t Zl(203>. Z2 ( C 
1203), Z3 ( 203 ) , Z4(203)« DATE(2)« DRM(4>, DRSD(4). RAN(4.201)* RN(4 C 
2), RSD ( 4 ) C 

COMMON /INTCOMM/ T, H, I NT, NEQ, I SCHEME, DERI NT (2, 249) C 

COMMON /INTINTR/ 1 NTERN ( 5 , 249 ) C 

COMMON /REALTlM/ ADC (32) ,DAC(64) ,LDISI ( 1 08 > ,LDISO( 1 96 ) ,NOPER.NHOLD C 
1 .NRESET.NTERM.nPRINT.NREAD C 

COMMON /'ALGOR/' NPAR. IPAR. INTP130 ) , IP. INTV (8) . I V. INTA ( 1 1 > . I A, I A1 , IA C 
12.PARAM (40) , DPARAM (30 ) , ALV ( 1 1 ) ,DVAR(8) , DALG ( 1 1 ) «ALG(40) , I AC (40) . IE C 
2VEN.WT (11,11) , COM. LI ,L2 C 

common /CO MM1/ IRR. IPL, I SL . TABLE ( 1 99 ) , I NTEG ( 99 ) . LOG I C ( 20 ) . NTAB. NI N C 
1 T,NLOG,NADC,NDAC.NLDI ,NLOO,NT,AXI ,AYI ,AZI .DRAD.RADD.PI . IR(2,),NRN.N C 
2PTS. I SKIP. JSKIP.KSKIP.MF.TX.TY.TZ.AMX.AMY.AMZ.VARCHNG. I TYPE. IVARBU C 
3F(5) ,DELX»NOITSPS,PTSINV C 

COMMON /C0MM2/ MAXPAGE .LAST ,LABU ,L ABV »L ABW , L ABP ,L ABO , LABR ,LABTH ,LA C 
1 BPH.lABAX.LABAY.LABAZ,LABDA,LABDE»LA8DR.FET25( 17) .LUN25.NAM25, I TRN, C 
2MLT ( 7 ) « I MULTA ( 7 ) , I NVSEN ( 7 ) , I MULTB ( 7 ) , I NVWT ( 7 > C 

COMMON /FLIGHT/ UM ( 201 ) . VM ( 20 l ) , WM ( 20 1 ) ,PM (20 1 ) , QM < 201 ) »RM ( 20 1 ) » TH C 
1 M(20l > ,PHM(20l ) , A XM ( 2 0 1 ) » A YM ( 20 1 ) , AZM (201 ) .FDA (201 ) .FOE (201 ) .FDR (2 C 
201 ) , FBT ( 20 1 ) , FDB ( 20 1 ).FDELB(201 ),FPPER(20l ).FTPER(201 ),.BETAT.ETA.E C 


3TAT.FCT ( 12) .CTT.CT 

EQUIVALENCE (DERI NT ( 1 . 1 ) ,U ) . ( DERI NT < 1 .2 > • V > » ( DER I NT (1 , 3 > , W ) , (DE 

1 R I NT (1 »4),P). (DERINT( 1.5) ,Q) . ( DER I NT ( 1 »6 ) «R ) » (DER I NT ( 1 .7 ) . THE ) ... 

2 (DERINT(1 , 8 ) .PHI ), (LD I SO ( 6 1 ) * WL ( 1 ) > , ( FSS ( 1 ) .LDI S 1 ( 33 )). (TABLE ( 

350), UMAX), (TABLE (51 ) .VMAX ) , ( TABLE < 52 ), WMAX ) . (TABLE ( 53 ) ,PMAX) . ( 

4TABLE < 54 ) .QMAX) , (TABLE (55 ) , RMAX ) « (TABLE(56) .THMAX) . - (TABLE(57).P 
5HMAX). (TABLE (58 ) .AXMAX ) . ( T ABLE < 59 ) . A YMAX ) . ( TABLE ( 60 ) .AZMAX ).( T 

6 ABLE (62 1 .DAMAX) » (TABLE (63 ) .DEMAX ) . ( TABLE ( 64 ) .DRMAX) . (TABLE (98 ) » 

7RUN), (TABLE(99) .PASS) . (TABLE ( 1 01 ) »0RM { 1 >) » (TABLE ( 1 7 1 ) . UCRTBI > , 

8 (TABLE ( 191 ) ,DRSD( 1 ) ) , ( INTEG (62 ) .NOPTS ) , ( INTEG(64 ) ♦ IPR1NT ) , ( I NTE 

9G (65 ) . NPLOT ) 

IPRlNT=l TO PRINT ICS 

TO REWIND. MF FILE . ■ 

TO OBTAIN RANDOM NUMBERS 

FOR CALCOMP PLOT (SET VALUE BEFORE EXITING CRT LOOP AND 
PRINT - USES FSS( 1 ,2.7.8, 13. 14) ) 

TO RETURN (ROUTES OUTPUT) 

INDICATES DEFAULT VALUE USED FOR IPRINT 
INDICATES CALCOMP PLOT COMPLETED 
INDICATES ERROR IN CALCOMP PLOT LOOP 


C 

C I PR I NT = 2 

C I PR I NT = 3 

C I PR I nT “4 

C ENTERING 

C IPRINT=5 

C WL ( 1 ) 

C WL ( 2 ) 

C WL ( 3 ) 

L* I PR I NT 
1 PR I nT = 5 
WL ( 1 )=.F. 

IF ( (L.GE.l ).AND. (L.LE.5) ) GO TO 10 
WL ( 1 )*.T. 

RETURN 

10 GO TO (20.30.40,120.340). L 

20 IEVEN=-1 

RUN»RUN+1 . 


PASS=-1 • 


WRITE 

(MF, 

■ 350 ) 

WR 

ITE 

(MF, 

,360) 

WR 

ite 

(MF, 

,370 ) 

WR 

ite 

( MF, 

,380 ) 

WR 

ite 

( MF, 

,390 ) 

WR 

ITE 

(MF, 

,400 ) 

WR 

ite 

(MF, 

>410 ) 

WR 

ite 

( MF, 

,420 ) 

WR 

ITE 

(MF, 

,430 ) 

WR 

ITE 

( MF, 

>440 ) 


RUN. ( (L1.1.L2 .LOG IC(1>)«I S 1«11 ) 
DVAR ( 1 > , ( ( COM .DVAR ( I ) > . 1=2. IV) 
DALG ( 1 ) , ( (COM.DALG ( I > ) . 1=2. I A > 

( ( PARAM ( I ) .TABLE ( I ) . I NTEG (!)>,!= 
( T ABLE ( I ) « 1=65.92 > 

(TABLE ( I ). 1=1 17, 133) 

( TABLE ( I ). 1=154. 162 ) 

(TABLE! I ) .1=169,186) 

TX.TY « TZ . AMX. AMY , AMZ 
INTEG (62). I NTEG (63) 


1 .NPAR) 


RETURN 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
001 1 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 
0061 
0062 

0063 

0064 
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30 


40 

50 

60 

TO 

80 

90 
1 00 

I 10 
120 

C 

C 

130 

140 

150 


REWInO MF 
RETURN 

CAUL GETRAN ( IR« 1 .2.DMY. Y1 .Y2> 

DO 50 I “1 «NRN 
DO 50 J-t .NPTS 

CALL GETRAN ( I RV2.2 . RAN (I.J).Yl .Y2 > 

CONTINUE 

DO 60 ! “1 *NRN 

RN( I )=0* 

RSD (11=0# 

DO 60 J = 1 «NOPTS 
RN( I )*RN( I ) +RAN ( I «J) 

RSD ( I >*RSO< I 1+RANU . J)**2 
DO 70 1=1 tNRN 
RNtl) =RN t 1 ) *PTS 1 NV 

RSD ( I >*SQRT<RSD< I )*PTS1NV-RN< I )**2> 
RSD ( t 1 *DRSD ( 1 ) /RSD ( I ) 

RN ( I ) =DRM ( 1 l-RSD ( 1 ) *RN ( 1 ) 

00 80 1 *1 »NRN 
DO 80 J*1 « NOPT S 

RAN( I • J)*RAN< I « J)*RSD< I >+RN< I > 

CONTINUE 

DO 90 I =1 «NRN 

RN( I )=0» 

RSD< 1 )»0. 

DO 90 J“1 «NOPTS 

RN ( I )=RNCI ) +RAN ( I .J) 

RSD< I 1 *RSD < 1 >+RAN(I « J 1**2 
DO 100 1*1* NRN 
RN< I )=RN( I )*PTSINV 

RSD < I ) * SORT ( RSD ( I )*PTS 1 NV-RN < 1 1**2) 
CONTINUE 

DO 110 1 = 1* NOPTS 
UM( I )*UM< I )+RAN(l , I > 

WM ( I )=WM( I ) +RAN <2*1 ) 

QM(I )=QM( I)+RAN<3» t 1 

THM( I > = THM( I ) +RAN ( 4 * I J 

WRITE ( MF 1 450 ) RUN « NOPTS » < < ORM < I > « RN 

RETURN 

WL(3)=.F. 

IF (JSKIP.GT.O) GO TO 130 
CALL CALPLT (0.«0»«-3) 

CALL LEROY 

Ke-i FOR SYMBOL EVERY DATA POINT 

K*-2 FOR SYMBOL EVERY OTHER 

Kb-nplOT 

JSKIP=JSKIP+1 

BI=SF4=0. 

IF (,N0T.FSS (1 ) ) GO TO 140 

SF1=aZMAX 

SF2=AYMAX 

SF3*AXMAX 

GO TO 1 TO 

IF ( .NOT.FSSIT) > GO TO 150 

SFlsPHMAX 

SF2=v*AX 

SF3*RMAX 

SF4*PMAX 

GO TO 170 

IF (,N0T.FSS(8) 1 GO TO 160 

SFl*DRMAX 

SF2*oEMAX 

SF3=dAMAX 

GO TO 170 


C 0065 
C 0066 
C 0067 
C 0068 
C 0069 
C 0070 
C 0071 
C 0072 
C 0073 
C 0074 
C 0075 
C 0076 
C 0077 
C 0078 
C 0079 
C 0080 
C 0081 
C 0082 
C 0083 
C 0084 
C 0085 
C 0086 
C 0087 
C 0088 
C 0089 
C 0090 
C 0091 
C 0092 
C 0093 
C 0094 
C 0095 
C 0096 
C 0097 
C 0098 

v C 0099 

C 0100 
C 0101 

( I ) »DRSD( I ) .RSD( I ) ) , 1*1 .NRN) C 0102 

C 0103 
C 0104 
C 0105 
C 0106 
C 0107 
C 0108 

PO I NT «ETC • C 0109 

C 0110 
C 0111 
C 0112 
C 0113 
C 0114 
C 0115 
C 0116 
C 0117 
C 0118 
C 0119 
C 0120 
C 0121 
C 0122 
C 0123 
C 0124 
C 0125 
C 0126 
C 0127 
C 0128 
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160 SF1=THMAX 
SF2*OMAX 
SF3=WMAX 
SF4=UMAX 

bi*ucrtbi 

170 ZT *81 +SF4 
Z8*B1 -SF4 

IF <FSS<2).0«.FSS<8> ) GO TO 220 

1*0 

J=-l 

ISO CALL PLAYBAK (21 5S 1 
J*J+! 

IF {M0D(J«N01TSPS).GT.01 GO TO 180 
1*14-1 

IF (I.GT.NOPTS) GO TO 180 
IF ( .NOT.FSSI 1 ) ) GO TO 190 
Y1 ( 1 ) *AZ I 
Y2 ( I )=AYI 
Y3< I > =AX I 
Y4 ( I )*0* 

GO TO a 10 

190 IF (.NOT.FSSI7) ) GO TO 200 
Y 1 ( 1 ) *PH I 
Y2 ( t )'V 
Y3 (I ) *R 
Y4 { I ) *P 
GO TO 210 

200 Y 1 ( 1 )*THE 
Y2 < I > *Q 
Y3< I )*W 
Y4(I1*U 

2 1 0 IF ( Y1 I t 1 *GT «SF1 ) Y t ( I 1 *SF 1 
IF (Y2( I 1.GT.SF2) Y2 < I ) *SF2 
IF (Y3< I 1.GT.SF3) Y3 ( I ) *SF3 
IF (Y4( I J.GT.ZT) Y4 ( I )*ZT 
IF (Yl < I l.LT.-SFl ) Y 1 ( 1 1 =— SF 1 
IF (Y2( I 1.LT.-SF2) Y2(I)=-SF2 
IF (Y3< I 1.LT.-SF3) Y3<I)=-SF3 
IF fY4< I >.LT.ZB> Y4 ( I ) *ZB 
GO TO 180 

215 IF (I.EQ.NOPTS) GO TO 220 
WL <3) = «T, 

RETURN 

220 DO 270 I * 1 ♦ NOPTS 
X< I »=OELX*( I-! ) 

IF < .NOT.FSSI 1 ) ) GO TO 230 
Z1 < I )=AZM< 1 ) 

Z2 ( I i = AYM ( I ) 

Z3(n=AXM(I) 

Z4 (I ) *0 • 

GO TO 260 

230 IF l ,NOT,FSS(7 ) ) GO TO 240 
Z1 ( 1 )*PHM< I ) 

Z2 < I )*VM(I ) 

Z3 ( I )*RM( I > 

Z4 ( I )*PM< I ) 

GO TO 260 

240 IF < .NOT «FSS <8)1 GO TO 250 
Z1 ( 1 ) =FDR ( 1 ) 

Z2 ( I >=FDE(I ) 

Z3 ( I ) *FDA ( I 1 
Z4 { I )*0* 

GO TO 260 


C 0129 
C 0130 
C 0131 
C 0132 
C 0133 
C 0134 
C 0135 
C 0136 
C 0137 
C 0138 
C 0139 
C 0140 
C 0141 
C 0142 
C 0143 
C 0144 
C 0145 
C 0146 
C 0147 
C 0148 
C 0149 
C 0150 
C 0151 
C 0152 
C 0153 
C 0154 
C 0155 
C 0156 
C 0157 
C 0158 
C 0159 
C 0160 
C 0161 
C 0162 
C 0163 
C 0164 
C 0165 
C 0166 
C 0167 
C 0168 
C 0169 
C 0170 
C 0171 
C 0172 
C 0173 
C 0174 
C 0175 
C 0176 
C 0177 
C 0178 
C 0179 
C 0180 
C 0181 
C 0182 
C 0183 
C 0184 
C 0185 
C 0186 
C 0187 
C 0188 
C 0189 
C 0190 
C 0191 
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250 Z1 ( I >»THM( I ) 

Z2 ( ! )»QM( I ) 

Z3(I)«WM(!) 

Z4 ( I )*UM(1 ) 

260 IF <Z! ( I l.GT.SFl > Z1 ( I >*SF1 
IF (Z2( I ).GT.SF2) Z2(I)=SF2 
IF (Z3< I ).GT.SF3) Z3(I)=SF3 
IF (Z4( I ).GT.ZT> Z4 ( I } *ZT 
IF (Z1 < I ) .LT.-SFJ > Z1(I)=-SF1 
IF (Z2C I ).LT.-SF2) Z2(I)«-SF2 
IF (Z3( I ).LT.-SF3) Z3(I)*-SF3 
IF (Z4( I l.LT.ZB) Z4(I)=ZB 
270 CONTINUE 

CALL DAYTIM (DATE I 

CALL NOTATE (- I t5< 0 .0 . . 07, DATE ( 1 ) , 90. . 1 0 ) 

CALL NOTATE ( - 1 . 5 .2.4 « . 07. DATE (2 ) .90. • 1 0 } 

CALL NOTATE (- 1 .5 <4 .8. • 07 « 7HRUN N0-.90..7) 

CALL NUMBER (-1 .5,6.5, .07.RUN.90. »-l > 

CALL NOTATE ( - 1 .5 .7.5 . . 07, 5H I TER= .90. » 5 ) 

CALL NUMBER ( - 1 .5 . 8. 7. . 07 . PASS . 90 • * - 1 > 

CALL ASCALE (X.5. .NOPTS. I , 1 0. ) 

LPT=nOPTS+1 

UPT*nOPTS+2 

XM=X(LPT> 

XS*X < JPT ) 

CALL GRID (0. .0. . .5. .5. 10.4 ) 

CALL AXES <0. .0. ,0. .5. .XM.XS.l •» 1 0. .LAST, . 14.-10) 

ym*-sfi 
YS »SF 1 

Y1 <LPT)*Z1 (LPT)-YM 
Y1 (JPT)»ZI (JPT)=YS 
IF ( .NOT.FSSI 1 ) > GO TO 280 

CALL AXES ( O. , 0 . . 90. .2. .YM.YS. .5 * 10.. 12 NAZI (G UN I TS 1 . . 1 4 . 1 2 > 
CALL LINE (X.Zt .NOPTS. 1 .K.3. .07) 

IF ( ,NOT.FSS(2 ) ) CALL LINE < X . Y1 .NOPTS. 1 . 0 . 0 ♦ .07 1 
CALL CALPLT (0..2.75.-3) 

CALL GRID (0..0.. .5, .5.10.4) 

CALL AXES (0. .0. ,0. .5. .XM.XS.l .. 1 0. .LAST. . 14.-10 ) 

YM=— SF2 
YS=SF2 

Y2 (LPT ) «Z2 ( LPT ) = YM 
Y2 < JPT ) =Z2 ( JPT ) *YS 

CALL AXES (0. ,0. .90. .2. . YM. YS..5. 1 0. . 12HAYI (G UN I TS ) « • 1 4 , 1 2 ) 
CALL LINE (X.Z2.N0PTS. 1 .<.3. .07) 

IF ( ,NOT.FSS(2 ) ) CALL LINE ( X . Y2 . NOPTS . 1 . 0 , 0 . .07 ) 

CALL CALPLT (0..2.75.-3) 

CALL GRIO (0..0...5..5.10.4) 

CALL AXES (0. .0. .0. .5. .XM.XS. 1 •* 1 0. .LABT. • 14.-10 ) 

YM=-SF3 

YS=Sf3 

Y3 (LPT ) «Z3 (LPT ) =YM 
Y3 ( JPT ) «Z3 ( JPT ) «YS 

CALL AXES (O. .0..90. .2. .YM.YS..5.10.. 12HAXI (G UNI TS ) . . 14 . 1 2 ) 
CALL LINE (X.Z3. NOPTS. 1 ,K,3. .07) 

IF ( ,NOT.FSS(2 ) ) CALL LINE (X. Y3 . NOPTS . 1 . 0 . 0 .. 07 ) 

CALL CALPLT ( 1 2. ,-5.5. -3 ) 

GO TO 310 


C 0192 
C 0193 
C 0194 
C 0195 
C 0196 
C 0197 
C 0198 
C 0199 
C 0200 
C 0201 
C 0202 
C 0203 
C 0204 
C 0205 
C 0206 
C 0207 
C 0208 
C 0209 
C 0210 
C 0211 
C 0212 
C 0213 
C 0214 
C 0215 
C 0216 
C 0217 
C 0218 
C 0219 
C 0220 
C 0221 
C 0222 
C 0223 
C 0224 
C 0225 
C 0226 
C 0227 
C 0228 
C 0229 
C 0230 
C 0231 
C 0232 
C 0233 
C 0234 
C 0235 
C 0236 
C 0237 
C 0238 
C 0239 
C 0240 
C 0241 
C 0242 
C 0243 
C 0244 
C 0245 
C 0246 
C 0247 
C 0248 
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280 


IF ( .NOT.FSS(7) ) GO TO 290 

CALL AXES ( 0 . .0 . . 90* <2. « YM, YS. .5 » 10* *LABPH« *14*10) 

CALL LINE CX.Zl .NOPTS. 1 .K.3. .07) 

IF ( .NOT .FSS (2 ) ) CALL LINE ( X. VI .NOPTS. 1 . 0. 0. . 07 ) 

CALL CALPLT (0..2.75.-3) 

CALL GRID (0. .0. « .5. .5. 10.4 ) 

CALL AXES (0. .0. ,0. .3. .XM.XS.l •« 1 0. .LABT. • 14,-10 ) 

YM=— SF2 
YS=SF2 

Y2 (LPT ) =Z2 ( LPT ) =YM 
Y2 ( JPT > *Z2 ( JPT ) ■' YS 

CALL AXES (0. .0. .90. .2. ,YM. YS. .5, 1 0. ,LABV«. 14, 10 ) 
VMIN=YM*.3048 
VMX=YS*.304B 

CALL AXES (3. .0* .90. .2. .VMIN.VMX. .5. 1 0. ,9HV (MXSEC ) . . 1 4 . -9 ) 
CALL LINE (X.Z2. NOPTS. 1 .K.3. .07) 

IF ( .NOT .FSS 1 2 ) ) CALL LINE ( X . Y2 .NOPTS . 1 . 0 . 0. .07 ) 

CALL CALPLT (0. .2. 75.^3 ) 

CALL GRID (0..0...5..5.10.4) 

CALL AXES (0..0. .0. ,5. .XM.XS.l .. 10. .LAST. • 14,-10) 

YM=-SF3 
YS=SF3 

Y3 (LPT ) *Z3 ( LPT ) *YM 
Y3 ( JPT ) =Z3 ( JPT ) =YS 

CALL AXES (0. .0..90. .2. .YM.YS. .5. 10. «LABR. . 14, 10 ) 

CALL LINE (X.Z3. NOPTS. 1 .K.3. .07) 

IF ( .NOT .FSS(2) ) CALL LINE { X. Y3 .NOPTS. 1 . 0 .0. • 07 ) 

CALL CALPLT (O., 2.75. -3) 

CALL GRID (0. ,0. « .5. *5. 10,4 ) 

CALL AXES (0..0..0..S..XM.XS.1..10..LABT..14.-10) 

YM*-SF4 
YS=SF4 

Y4 (LPT ) *Z4 (LPT ) *YM 
Y4 ( JPT ) =Z4 ( JPT ) * YS 

CALL AXES (0 . ,0 . .90. *2 » » YM » YS « .5 . 1 0. »LABP ,.14.10) 

CALL LINE ( X » Z4 , NOPT S.1.K.3..07) 

IF ( .NOT .FSS(2) ) CALL LINE ( X. Y4 .NOPTS . 1 . 0, 0. . 07 ) 

CALL CALPLT ( 1 2. . -8 .25. -3 ) 

GO TO 310 
290 IF ( .NOT .FSS (8 ) ) GO TO 300 

CALL AXES (0. ,0. .90. ,2. ,YM. YS. .5. 10. .LABDR, • 14 . 10) 

CALL LINE (X.Z1 .NOPTS. 1 .K.3. .07) 

CALL CALPLT (0..2.75.-3) 

CALL GRID (0..0.. .5. .5,10.4) 

CALL AXES (0. ,0. .0. ,5. .XM.XS. 1 ., 1 0. .LABT. . 14,-10 ) 

YM=— SF2 
YS»SF2 

Y2 (LPT ) =Z2 ( LPT ) =YM 
Y2 ( JPT ) *Z2 ( JPT ) =YS 

CALL AXES (0. .0..90. ,2. , YM.YS. .5, 10. .LABDE, . 14 .10 ) 

CALL LINE (X. Z2. NOPTS. 1 .K.3..07) 

CALL CALPLT (0..2.75.-3) 

CALL GRID (0. ,0. . .3, .5, 10.4 ) 

CALL AXES (0. .0. .0. .5. .XM.XS.l .. 1 0. .LABT, • 14.-10 ) 

YM*— SF3 
YS=SF3 

Y3 (LPT ) *Z3 (LPT ) =YM 
Y3 ( JPT ) *Z3 ( JPT ) *YS 

CALL AXES (0. ,0. .90. ,2. ,YM. YS. .5, 10. .LABDA, . 14 , 10 > 

CALL LINE (X.Z3. NOPTS. 1 .K.3. .07 > 

CALL CALPLT ( 1 2. . -5.5 . -3 ) 

GO TO 310 


C 0249 
C 0250 
C 0251 
C 0252 
C 0253 
C 0254 
C 0255 
C 0256 
C 0257 
C 0238 
C 0259 
C 0260 
C 0261 
C 0262 
C 0263 
C 0264 
C 0265 
C 0266 
C 0267 
C 0268 
C 0269 
C 0270 
C 0271 
C 0272 
C 0273 
C 0274 
C 0275 
C 0276 
C 0277 
C 0278 
C 0279 
C 0280 
C 0281 
C 0282 
C 0283 
C 0284 
C 0285 
C 0286 
C 0287 
C 0288 
C 0289 
C 0290 
C 0291 
C 0292 
C 0293 
C 0294 
C 0295 
C 0296 
C 0297 
C 0298 
C 0299 
C 0300 
C 0301 
C 0302 
C 0303 
C 0304 
C 0305 
C 0306 
C 0307 
C 0308 
C 0309 
C 0310 
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300 


3 1 0 
320 


330 

340 

350 

360 

370 

380 

390 


CONTINUE C 0311 

CALL AXES (0. .0. .90. .2. .YM. YS« .5 . 10. .LABTH, . 14. 10 ) C 0312 

CALL LINE (X.Zl .NOPTS. 1 .K.3. .07) C 0313 

IF ( • NO T.FSSI2) ) CALL LINE < X . Y 1 .NOPTS . 1 . 0 . 0 . • 07 ) C 031* 

CALL CALPLT I0..2.75.-3) C 0315 

CALL GRID (0. .0.. .5. .5. 10.4 ) C 0316 

CALL AXES (0. .0. .0* .5. .XM.XS. 1 • « 1 0. .LABT. . 14.-10 ) C 0317 

YM=- SF2 C 0318 

YS*Sf2 C 0319 

Y2(LPT)«Z2(LPT>*YM C 0320 

Y2(JPT)*Z2( JPTI^YS C 0321 

CALL AXES (0..0..90..2..YM.YS..5.10..LABQ..14.10) C 0322 

CALL LINE (X.Z2. NOPTS. 1 .K.3..07) C 0323 

IF (.NOT.FSS<2>> CALL LINE ( X . Y2 .NOPTS . 1 , 0 , 0 .. 07 ) C 0324 

CALL CALPLT (0..2.75.-3) C 0325 

CALL GRID <0.,0...5..5.10,4> C 0326 

CALL AXES (0. .0. .0. ,5. .XM.XS. 1 ., 10. .LABT, . 14.-10 ) C 0327 

YM=-SF3 C 0328 

YS=SF3 C 0329 

Y3(LPT)-Z3(LPT)=YM C 0330 

Y3 ( JPT ) =Z3 ( JPT ) = YS C 0331 

CALL AXES (0. .0. .90. .2. «YM. YS..5 • I 0. .LABW. • 1 4. 10 ) C 0332 

WMIN*YM*.3048 C 0333 

WMX=YS*.3048 C 0334 

CALL AXES <5. .0. .90. .2. .WMIN.WMX. .5. 10. .9HW ( M/SEC > . . 1 4 , -9 ) C 0335 

CALL LINE ( X. Z3 .NOPTS. 1.K.3..07) C 0336 

IF ( .NOT «FSS (2 ) ) CALL LINE (X . Y3 . NOPTS . 1 . 0 . 0 .. 07 ) C 0337 

CALL CALPLT (0..2.75.-3J C 0338 

CALL GRID <0. .0. . .5. .5. 10,4 ) C 0339 

CALL AXES (0. .0. .0. ,5. .XM.XS. 1 •• 1 0. .LABT. . 14,-10 ) C 0340 

YM=ZR C 0341 

YS=SF4 C 0342 

Y4(LPT)=Z4(LPT)=YM C 0343 

Y4 I JPT ) =Z4 ( JPT ) »YS C 0344 

CALL AXES (0. .0. .90. .2. . YM. YS. .5 . 1 0. .LABU. . 14. 10 ) C 0345 

UMIN»YM*.3048 C 0346 

UMX=YS#.3048 C 0347 

CALL AXES (5. .0. .90. .2. .UMIN.UMX, .5. 1 0. ,9HU < M/SEC ) . . 1 4 .-9 > C 0348 

CALL LINE <X.Z4, NOPTS. 1 .K.3..07) C 0349 

IF ( .NOT .FSS ( 2 1 ) CALL LINE (X . Y4 . NOPTS , 1 . 0 , 0 .. 07 > C 0350 

CALL CALPLT ( 1 2. . -8.25 • -3 > C 0351 

WL(2)«.T. C 0352 

CALL DSPLAY C 0353 

CALL OPERATE C 0354 

IF ( .NOT .FSS (141) GO TO 330 C 0355 

WL<2)*.F. C 0356 

CALL OPERATE C 0357 

GO TO 130 C 0358 

IF ( .NOT .FSS (131) GO TO 320 C 0359 

WL < 2 ) * • F • C 0360 

RETURN C 0361 

FORMAT (/V/6H RUN=.F4.0//(6(A8. 12. A2.L1 ) > ) C 0362 

FORMAT (X39H ACTIVE EQUATION VARIABLES ARE . . . . A3. 7 ( A 1 , A3 ) ) C 0363 
FORMAT I/39H ACTIVE ALGORITHM VARIABLES ARE ... .A3. 1 0 ( Ai , A3 ) ) C 0364 

FORMAT (/5I25H PARAM VALUE ACT > / ( 5 ( 2X A5. 1 XE1 3.6 . 1 XL3 ) ) ) C 0365 

FORMAT </9H UMULT *.E13.6.9H VMULT =.EI3.6.9H WMULT *.E13.6.9H C 0366 

1 PMULT “.E13.6.9H QMULT =.E13.6.9H RMULT =.E13.6/9H THMULT * . El 3 C 0367 

2.6. 9H PHMULT =,E13.6.9H AXMULT =»E13.6.9H A YMULT* . El 3. 6 . 9H AZMUL C 0368 
3T=.El3.6.9H OAMULT = ,E 1 3 .6/9H DEMULT*= , E 1 3. 6 , 9H DRMULT = » E 1 3 • 6 » 9H C 0369 

4 UBIAS -.E13.6.9H VBIAS *.E13.6.9H WB1AS *,EI3.6.9H PBIAS =,E13 C 0370 
5.6/9H QB IAS -.E13.6.9H RBIAS *,E13.6.9H THB I AS“ .El 3.6 . 9H PHBIA C 0371 
6S*.El3.6.9H AXB1AS=»E13.6.9H AYBI AS= . E 1 3.6/9H A ZB IAS-.E13.6. 9H C 0372 

7 DAB IAS“»E13.6. 9H DEB I AS* » E 1 3. 6 . 9H DRB I As= . E 1 3 . 6 ) C 0373 
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400 FORMAT </9H AIX ®.E13.6«9H A!y =.E13.6.9H AIZ =.E13.6.9H C 0374 

1 AtXZ =.E13.6.9H WE I GHT= « E 1 3. 6 . 9H GRAV =.E13.6/9H RHO ®,E13 C 0375 

2.6.9 H S S.E13.6.9H 8 =.E13.6.9H CBAR =.EI3.6.9H DEAMP C 0376 

3L® *E 1 3. 6.9H OEFREQ® ,E 1 3.6/9H ALPHAT = . E 1 3. 6 . 9H DT ®.EI3.6.9H C 0377 

4 TT ®.E13.6.9h TS =.E13.6.9H TIMF =.E13.6> C 0378 

410 FORMAT </9H XX ®,E13.6»9H Yx =.E13.6.9H ZX =.E13.6.9H C 0379 

1 XY *.E13.6.9H YY =.E13.6.9H ZY =.E13.6/9h XZ ®.E13 C 0380 

2.6.9H YZ = »E I 3.6 . 9H ZZ =.E13.6> C 0381 

420 FORMAT (/9H CMCON »,E13.6.9H D ALML T= * E 1 3 • 6 « 9H UCRTB I = , El 3 .6. 9H C 0382 

1 AIW =.E13.6.9H CTO =.E13.6.9H CTTO ®.E13.6/9H CTB ®.E13 C 0383 

2.6.9H CTBT =,E13.6.9H D =*E13.6.9H CAPDT =.E13.6.9H ELBAR C 0384 

3 ® »E 1 3. 6 «9H ELTP ®,E13.6/9H RB ®»E13.6.9H PPER ®.E13.6.9H C 0385 

4 TPER =.E1 3.6.9H BTB I AS® . El 3.6 , 9H PPERB I » . E 1 3 . 6 . 9H TPERBI = ,E!3 C 0386 

5.6) C 0387 

430 FORMAT (/9H TX =.E13.6.9H TY =.E13.6.9H TZ *.E13.6.9H C 0388 

1 AMX ® »E l 3.6 »9H AMY ®.E13.6.9H AMZ *.E13.6) C 0389 

440 FORMAT </8H NOPTS® .14* 6H INC®. 14) C 0390 

450 FORMAT (✓✓/6H RyN«.F4.0«8H NOPTS® . I 4///5X. 6HDRM ( I) . 9X. 5HRN (I) . 9X C 0391 

1 »7MORSD< I ) »9X.6HRSD< I )//<4 (2XE13.6 ) > ) C 0392 


ENO 


C 0393 


OVERLAY (XC142FL,3,0) .~ The read overlay, level (3,0), is loaded when the READ 
mode is selected on the program control console. The primary overlay, level (4,0), is 
automatically loaded upon completion of level (3,0). The integer IREAD preset in level 
(4,0) selects one of the six options in the READ overlay, one option being to input flight 
test data. 

The flight test data input option uses a tape that was made by altering the original 
time records of all runs so that monotonically increasing time serves as the tape index 
key. Another aspect of the tape is that the time interval between consecutive points is 
constant for a given run, but does vary from run to run. The analyst can readily deter- 
mine TS, TT, and NOPTS (tape starting time, tape time interval, and number of points, 
respectively) from the tape printout, but must consider frequencies involved to determine 
INC (sample rate of flight data tape). A major factor in determining INC is that DELX 
(program time interval of flight data, where DELX = TT . INC) must be an integral 
multiple of DT (program integration step size). 

This option also provides the analyst with the opportunity to use biases and multi- 
pliers to eliminate certain apparent anomalies in the flight test data. An example of this 
usage is the putting in of control input trim conditions as biases. 

To input flight test data, the determined values for TS, TT, NOPTS, INC, and any 
multipliers and/or biases must be entered into the computer. Then IREAD is changed to 
6 and the READ switch is depressed. Input of flight test data is completed when the READ 
mode light comes on; if WL(7) is also on, an error in entering these values is apparent. 
The bottom RELEASE switch is depressed to complete the return to level (4,0). 

The flow chart of the flight test data input option follows the flow chart of overlay 
level (3,0). 
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FLIGHT TEST DATA INPUT OPTION 


Entry when IRE AX) = 6 


WL(7) = .F. 


(DO 130) j Yes 


Rewind LUN25 
UM(I) = 1000000 
I - 1, 2, . . NPTS 


(D0136) * 

JaJ+l 


(D0137) 4. 

^ Call RECIN 


END FILE X, Yes 
. LUN25 


(D0143)>^X 

NoX S , artin6 


time 

reached^. 


— <^MOD(J,INC) = (T> 

(D0145) Ye °| 

UM(I) = UMULT ■ (UM<I) - UBIAS) 

(see lines D0145 D0162) 
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c 

c 

c 

c 

c 

c 

c 

c 


overlay < xc 1 42fl • 3 ♦ o j 
program xcread 

logical LDISI .LDISO.LOGIC.VARCHNG.WL (39) ,FSS( 16) 

DIMENSION FWA25( 1025). AL(40). INTX<8). INTY(II) 

COMMON / 1 NT COMM/ T « H . I NT .NEQ . I SCHEME .DER I NT < 2 , 249 ) 

COMMON /INTINTR/ INTERN (5. 24 9) 

COMMON /REALT IM/ ADC (32 > .DAC ( 64 ) . LDI S I ( 1 08 ) «LD t SO ( 1 96 ) .NOPER • NHOLD 
1 .NRESET.NTERM.NPRINT.NREAD 

COMMON /ALGOR/ NPAR » I PAR. INTP(30 ) . IP. INTV (8 ) .IV.INTA(11).IA.IA1.IA 
12.PARAM (40) .DPARAM ( 30 ) . ALV ( 1 1 ) .DVAR(S) «DALG ( 1 I ) ,ALG (40 ) . I AC (40) . IE 
2VEN.WT ( 1 1.11) .COM. LI .L2 

COMMON /COMM 1 / I RR . I PL . I SL ♦ T ABLE ( 1 99 ) . I NTEG ( 99 ) .LOG I C ( 20 ) .NTAB.NIN 
1T.NL0G.NADC.NDAC.NL0 I .NLDO.NT.AXI ,AYI .AZI «DRAD .R ADD. PI . IR (2) . NRN.N 
2PTS. I SKIP, JSKIP.KSKIP.MF, TX. TY. TZ , AMX . AMY. AMZ, VARCHNG. I TYPE. I VARBU 
3F(5),0ELX,NOITSPS.PTSINV 

COMMON /C0MM2/ MAXPAGE.LABT.LABU.LABV.LABW.LABP.LABQ.LABR.LABTH »LA 
1 BPH , L AB AX . L ABA Y « L A8AZ , L ABDA « LABDE .LABDR , FET25 ( 17 ) . LUN25 . NAM25 . I TRN 
2MLT ( 7 ) . IMULTA (7 ) « !NVSEN(7) . I MULTB ( 7) . I NVWT ( 7 > 

COMMON /FLIGHT/ UM ( 20 1 ) , VM(20l ) . WM (201 ) .PM (201 ) , QM (201 ) .RM (201 ) .TH 
1 M ( 20 1 ) .PHM(201 ) , AXM( 201 ) . AYM (201 ) , AZM (201 > .FDA (201 ) .FDE (20 1 ) .FDR (2 
201 ) .FBT (201 ) .FOB (201 ) , FDELB ( 201 > .FPPER (201 ) .FTPER ( 201 ) .BE TAT .ETA.E 
3TAT »ECT (12) .CTT.CT 

EQUIVALENCE (LD I SO ( 61 ) « WL ( 1 > ) » (LD I S I ( 33 ) . FSS ( 1 ) ) , (AL(l).UO). (AL 
1(5), THEO ) » ( AL ( 6 ) .PH I 0 > » ( AL ( 7 ) » WO ) . (AL(12).Q0)» <AL(18).V0). (AL 

2(25). PO). ( AL ( 33 ) ,R0 ) » ( TABLE ( l > . AL ( 1 > > • ( TABLE ( 65 ) .UMULT ) . (TABLE 

3 ( 66 ) . VMULT ) , ( TABLE ( 67 ) .WMULT) . ( TABLE ( 68 ). PMULT ) . ( TABLE ( 69 ) .OMUL 

4T ) , ( TABLE ( 70 ) .RMULT ) , ( TABLE ( 7 1 > .THMULT ) . ( TABLE ( 72 > . PHMULT ) , (TA 

50LE (73 ) .AXMULT ) . ( TABLE ( 74 ). AYMULT ) . (TABLE ( 75 ) .AZMULT ) , ( T ABLE ( 76 

6).DAMULT), (TABLE (77 ) .DEMULT ) , ( TABLE ( 78 ) .DRMULT ) . ( TABLE ( 79 ) .UBI A 

7S), (TABLE ( 80 ) , VB IAS), ( TABLE (8 1 ) . WB I AS ) . ( TABLE ( 82 ) . PB I AS ) . (TABL 
BE ( 83 ) . OBI AS ) , ( TABLE ( 84 ) . RB I AS ) . (TABLE (85 ) . THB I AS ) . ( TABLE ( 86 > .PH 
9BIAS). (TABLE ( 87 ) « AXB IAS). ( TABLE (88 ) . AYBI AS ) . ( TABLE ( 89 ) . AZB I AS ) » 

* ( TABLE (90 ) . DAB IAS), ( T ABLE ( 91 ) .DEB I AS ) , ( TABLE ( 92 > . DRB I AS ) . (TABL 
*E(120).A1XZ>. (TABLE ( 123) .RHO) . { TABLE ( 1 24 ). S ) , ( TABLE ( 1 30 ). DT ) , ( 

STABLE (131). TT ) , ( TABLE ( 1 32 ). TS ) . ( TABLE ( 1 33 ) » T I MF ) . ( TABLE ( 1 7 1 ) .UC 

SRTBI). ( T ABLE ( 1 75 ) , C TB ) . { TABLE ( 1 76 ). CTBT ) , ( TABLE ( 1 78 ). CAPDT ) . (T 
SABLE ( 179 ) .ELBAR) , ( TABLE ( 1 80 ). ELTP ) . (TABLE ( I 81 ) , RB ) , ( TABLE (184), 

S8TB1AS). (TABLE ( 185 ) .PPERBI ) . (TABLE ( 1 86 ) . TPERB I > . ( I NTEG ( 51 ) . I NTY 

S(l)), ( INTEG(62 ) .NOPTS) . ( I NTEG ( 63 ) . INC ) . ( I NTEG ( 66 ) . I RE AD ) . ( I NTE 

SG(81 ) . I NTX ( 1 ) ) 

t READ 3 1 TO STORE ALG IN AL.ETC. 

STORE AL IN ALG. ETC. (USES FSS ( 1 2 > ) 

GET NOMINAL CONDITIONS 
GET LONGITUDINAL CONDITIONS 
GET LATERAL CONDITIONS 
FILL FLIGHT DATA ARRAYS 
INDICATES DEFAULT VALUE USED FOR IREAD 

INDICATES ATTEMPT TO READ FLIGHT DATA BEYOND END FILE 


TO 

TO 

TO 

TO 

TO 


10 

20 

30 

40 


lREAO*2 
I REA0*3 
IREA0=4 
IREAD*5 
IREA0*6 
WL (6) 

WL (7) 

L* IREAD 
IREAO-2 
WL (6 )*.F. 

IF ( (L.GE.l ).AND. (L.LE.6) ) GO TO 10 
WL ( 6 ) * • T • 

RETURN 

GO TO (20,40.60.80,110.140), L 
DO 30 1*1, NPAR 
AL ( I >*ALG ( I ) 

I NTEG ( I )*IAC( I ) 

RETURN 

DO 50 1=1 .NPAR 
ALG( I )=TABLE( I ) 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 
0009 
001 0 
001 1 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 
0061 
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50 


I AC ( I ) = I NTEG < I > 

IF (FSS ( 12 ) ) RETURN 
ALG < 1 )=UM< 1 ) 

ALG< 18 > =VM < 1 ) 

ALG < 7 ) *WM { 1 ) 

ALG ( 25 1 *PM ( 1 ) 

ALG.< 1 2 > =OM < 1 > 

ALG <33 ) *RM< 1 ) 

ALG ( 5 ) = THM ( 1 ) 

ALG<6 )=PHM( 1 > 

RETURN 
60 A 1X2*8000. 

RHO*. 0021 86 
S=534.4 
0T= • 05 
TT*.05 
TIMFtB. 

CTB*.B423 

CTBT«.573 

CAPDT*8. 

ELBAR*20.54 
ELTP*32.09 
RB*1 .6 
INC”? 

DO 70 1=154.162 
70 TABLE < I) =0. 

RETURN 

80 DO 90 I =1 . 1 I 

IF (I.LT.9) I NTX < 1 ) =0 
90 I NTY t I ) = 1 

I NTX < 1 ) = I NTX (3) = INTX<5)»INTX(7) = 1 

I NTY ( 2 ) = I NTY ( 4 ) * I NTY (6) = INTY(8)=I NTY { 1 0 > = 0 

D! =U0 

D2 = V»0 

D3=QO 

04= the o 

00 I oo I=1,NPAR 
AL ( I 1=0. 

INTEG< I 1 = 0 

IF (I.LE.17) 1 NTEG { I ) = 1 

100 CONTINUE 

I NTEG ( 6 ) =0 

UO*Df 

W0=D2 

00*03 

THEO=04 

RETURN 

1 10 00 120 1*1.11 

IF (I.LT.9) I NTX ( I ) *0 
120 I NTY ( I ) *0 

I NTX ( 2 ) * I NTX ( 4 ) = I NTX <6) = INTX<8) = 1 
I NTY ( 2 ) = I NTY < 4 ) * I NTY <6) = INTY<8) = I NTY < 1 0 ) = 1 

01 *UO 
02= VO 
D3*PO 
04*R0 
D5*PHlO 

DO 130 1*1 .NPAR 

AL < I ) *0 • 

I NTEG ( I ) "1 

IF (I.LE.17) INTEG ( I ) = 0 


0 0062 
0 0063 
0 0064 
D 0065 
D 0066 
0 0067 
D 0068 
D 0069 
0 0070 
D 0071 
D 0072 
D 0073 
D 0074 
D 0075 
0 0076 
O 0077 
0 0078 
D 0079 
D 0080 
0 0081 
D 0082 
D 0083 
0 0084 
0 0085 
O 0086 
D 0087 
0 0088 
0 0089 
D 0090 
D 0091 
0 0092 
D 0093 
D 0094 
D 0095 
D 0096 
D 0097 
D 0098 
D 0099 
D 0100 
D 0101 
0 0102 
0 0103 
D 0104 
D 0105 
D 0106 
0 0107 
O 0108 
O 0109 
D 0110 
D 0111 
0 0112 
D 0113 
0 0114 
O 01 15 
D 0116 
0 0117 
0 0118 
D 01 19 
0 0120 
0 0121 
D 0122 
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130 


140 

150 

160 

170 

180 

190 

200 

210 


CONTINUE D 

I NTEG ( 6 ) * 1 D 

UO=Ol D 

V0*D2 , . .... . ... .. . _ J ._. ■ D 

PO*D3 ’ D 

R0=D4 D 

PH I 0=D5 D 

RETURN D 

WL ( 7 ) = . F . D 

IF (KSKIP.GT.O) GO TO 150 D 

KSKIP=1 D 

CALL CREATEF ( LUN25 . FWA25 . 1 025.FET25 .NAM25 > D 

REWIND LUN25 D 

! DO 160 1*1 * NPTS D 

UM(I)=t.E6 D 

1*0 D 

J = -l D 

1=1+1 0 

J=J+1 0 

CALL RECIN (LUN25 » 1 « I COUNT » X »UM (I).VM(l). WM ( I ) *PM ( I ) »QM ( I ) « RM (I ) « T D 
1HM( n • PHM ( I ) » AXM ( I ) , AYM( I ) • AZM ( I > .FDA < 1 ) , FOE ( I ) »FDR ( I ) »FBT ( I ) »FDB ( D 
2 1 ) .FOELB ( I ) .FPPER ( I ) .FTPER ( I ) ) D 

IF (FNDFILE LUN25) 200.210 D 

WL (7 ) * • T* D 

RETURN D 

IF (X.LT. (TS-.001 ) ) GO TO 190 D 

IF (MOD< J, INC) «NE.O ) GO TO l 80 D 

UM( I )=UMULT*(UM( I )-UBI AS) D 

VM( I )=VMULT*(VM( I)-VBI AS ) D 

WM< I )=WMULT*(WM( I >-WB! AS) D 

PM< I )»PMULT*(PM( I )-PBI AS) D 

QM( I )=QMULT*(QM( I )-QBIAS) D 

RM ( I )*RMULT* (RM ( I)-RBI AS) D 

THM (!) *THMULT* (THM < I ) -THB IAS) D 

PHM( I )=PHMULT*(PHM( l )-PHBI AS) D 

AXM ( 1 ) *AXMULT* ( AXM ( 1 )-AXB J AS ) D 

A YM ( r ) = AYMULT* (AYM( I )— AYBTAS) D 

AZM ( 1 )* AZMULT* (AZM( T )— AZBIAS) D 

FDA( I )=DAMULT*(FDA( I )-DA8IAS) D 

FOE ( I ) =DEMULT* (FDE ( I )— DEB IAS) D 

FDR ( 1 ) =DRMULT* (FDR ( I ) -DRB IAS) D 

FBT ( I ) =FBT ( I )-BT8IAS D 

IF(FDELB(I ).LT.O. ) FDELB ( I ) =0. D 

FPPER ( I ) =FPPER ( I ) -PPERB I D 

FTPER ( I ) “FTPER ( I ) -TPERB I D 

IF (I.LT.NOPTS) GO TO 170 D 

I *UM ( 1 ) D 

UCRTB I * I D 

UO=UM ( 1 ) D 

VO*VM( 1 ) D 

WO=WM( 1 ) D 

PO*PM( 1 ) D 

QO=OM( 1 ) D 

RO=RM ( 1 ) D 

TMEO=THM ( 1 ) D 

PHIO*PHM(l) D 

RETURN D 

END D 
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OVERLAY (XC142FL,4 > 0) .- Overlay level (4,0) nominally operates in real time and 
is automatically loaded upon exiting levels (1,0), (2,0), and (3,0). Level (4,0) can cause 
the loading of levels (2,0) or (3,0) by selecting PRINT or READ mode, respectively, on the 
program control console. Level (4,0) contains the maximum likelihood estimation pro- 
cedure (fig. 1) and CRT display loop. 

The CRT display loop has been developed to present time history comparisons 
between flight data and calculated results. The performance index variables and control 
deflections are divided into four separate displays (each selectable from the program 
control console). The displays consist of multiple grids and annotated axes with 

(1) A (.) symbol to represent flight data points 

(2) Continuous vectors between calculated data points. 

Displays available are selected as shown in the following table: 


Display 

FSS(l) 

FSS(7) 

FSS(8) 

Longitudinal (u, w, q, 9) 

.F. 

.F. 

.F. 

Lateral (p, r, v, (p) 

.F. 

.T. 

.F. 

Accelerometer (a x j, a Y p a z j) 

.T. 

.F. 

.F. 

Controls (6 a , 6 e , 5 r (flight data only)) 

.F. 

.F. 

.T. 


The CRT display loop is entered by any of the following methods: 

(1) Depressing FSS(15) and then releasing when in the RESET mode 

(2) Depressing FSS(15) and then releasing when an operational error occurs 

(3) Automatically upon completion of any iteration 

(4) Automatically upon completion of pseudo data fill. 

All four methods result in the program awaiting operator action, but method (4) first 
clears the existing display and plots the selected display. The operator action awaited is 
that of either requesting another display or exiting the CRT display loop. Requesting 
another display is accomplished by setting the appropriate FSS and depressing FSS(14) 
and then releasing. Exiting the CRT display loop is accomplished by depressing FSS(13). 

The flow chart for the CRT display loop follows the flow chart for overlay level (4,0). 
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CRT DISPLAY LOOP 
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0VERlAY(XC142FL,4.0) 

PROGRAM XCMAIN 
REAL LIX.MIY.NIZ 

LOGICAL FSS< 16) »LDISI .LDISO.LOGIC.VARCHNG.WL (39) 

DIMENSION AL ( 40 ) » 00ELA(30). DDX (11)* DPXTWDX ( 30 ) . DX ( 1 1 ) , INTX<8) 
1« INTY(II). PX ( 1 1 « 30 ) » PXD ( 6 *40 ) • PXTDX(30>. PXTPX ( 30 . 30 ) . PXTWT<3 
20.11). SO(lt.ll). XB AR (11).' DRM ( 4 ) • DRSD (4 ) . G(8.8). PFC8.40). DEL 
3A (40 ) 

COMMON /INTCOMM/ T.H, I NT. NEQ. I SCHEME. DERI NT (2. 249) 

COMMON /INTINTR/ INTERN (5.249) 

COMMON /REALTIM/ ADC (32 > .DAC (64 ) .LD I S I ( 1 08 ) .LD I SO ( 1 96 ) . NOPER « NHOLD 

1 .nreset.nterm.nprint.nread 

COMMON /ALGOR/ NPAR. I PAR. INTP (30 ) , IP . INTV (8 ) . I V, INTA ( 1 1 ) . I A . I A1 • I A 
12. PARAM ( 40 ) .OPARAM ( 30 ) . ALV( 1 1 ) .DVAR <8 ) « DALG ( 1 1 ) .ALG(40 ) * I AC (40) . IE 
2VEN.WT (11.11) .COM. LI «L2 

COMMON /COMM 1 / IRR. I PL . ISL.TABLE ( 199 ). 1NTEGC99 ) .LOGIC (20 ) .NTAB.NIN 
IT. NLOG . NADC . NDAC . NLD I . NLDO « NT . AX I . AY I .AZ I . DRAD .RADD.P I . IR (2) .NRN.N 
2PTS. I SKIP. JSKIP.KSKIP.MF.TX.TY.TZ. AMX. AMY . AMZ . VARCHNG . I TYPE. IVAR8U 
3F<5) .OELX.NOITSPS.PTSINV 

COMMON /COMM2/ MAXPAGE . LAST • LABU.L AB V «LABW . LABP ,LABQ ,LABR .LABTH ,LA 
1 BPH.LABAX «LABA Y .LABAZ . LABDA . LABDE . L ABDR . FET25 ( 1 7 ) , LUN25. NAM25 . I TRN 
2MLT(7) . IMULTA (7) . INVSEN(7 ) . I MULTB ( 7) . INVWT(7) 

COMMON /FLIGHT/ UM (201 ) .VM(201 > ,WM(201 ) .PM (201 ) .QM (201 ) .RM(201 ) .TH 
1M (201 ) « PHM ( 20 1 ) . AXM( 201 ) . AYM(201 ) «AZM (201 ) .FDA (201 ) .FDE(201 ) .FDR (2 
201 l.'FBTieoi ). FOB (201 ).F0£LB(201 ).FPPER(201 ).FTPER(20l ) .BETAT ,ETA ,E 
3TAT.FCT( 12) .CTT.CT 

EQUIVALENCE (PF( 1 ,3) .F13) , (PF ( 1 « 4 >. F 1 4 >,’ (PF ( 2 . 1 9 >. F2 1 9 > « (PF(2.2 
10),F220). ( PF (2.21 ).F221 ). (PF ( 2 . 22 ) .F222 ) , (PF (2 . 23 ) .F223 ) . (PF(2 
2 . 24 ) , F224 ) . ( PF (3.9) «F39 ) . (PF( 3. 10) .F31 O') , ( PF ( 3 . 11 ) .F3 1 1 > . (PF(4 
3 . 26 ) ,F426 ) . (PF (4.34 ) .F434 ) , ( PF ( 4 . 35 ) .F435 ) • (PF ( 4 . 36 ) .F436 ) . (PF 
4(4,37). F437 ) . (PF (4 . 38 ) . F4 38 ) . ( PF ( 4 . 39 ). F439 ) . (PF (6. 34 ) .F634 ) . ( 
5G (2,1) ,G21 ) , ( G ( 2 .2 ) «G22 ) , (G (2 .3 ) . G23 ) , (G ( 2 « 4 ) , G24 > . (G(2.6).G26 
6). (G(2«7).G27). (G ( 2 .8 ) . G28 ) , ( G ( 3. 1 ) . G3 1 ) . ( G ( 3 , 2 ) . G32 ) . (G (3 .3 ) 

7. G33). ( G ( 3 » 4 ) , G34 ) , ( G ( 3 . 5 ) .G35 ) . ( G ( 3 . 7 ) , G37 ) . ( G ( 3 . 8 ) . G38 > , (LD 
8IS0(6l 1 »WL( 1 ) ) , ( FSS ( 1 ) . LD I S I (33 ) ) » ( DER I NT ( l . 1 ) . U > . ( DER I NT ( 2. 1 ) . 
9UD0T ) . (DERINT( 1 .2) .V), ( DER I NT ( 2 . 2 ) . VDOT ) , (DERI NT ( 1 , 3 ) . W > . (DERI 
*NT (2,3) » WDOT ) , ( DER I NT ( 1 , 4 ) ,P ) . (DER 1 NT ( 2 .4 > . PDOT > . ( DER I NT ( 1 ,5 ) . Q 
*). ( OER INTC2.5). QDO T ) , ( DER I NT ( 1 , 6 ) , R > . ( DER INT ( 2 . 6 ) . ROOT ) . (DERIN 
*T (1.7). THE ) . ( DER INT (2.7) , THEDOT ) » ( DER I NT ( 1 . 8 ). PH I ) , (DERINT(2.8) 
S.PHIOOT). ( DER I NT (1,9). PS 1 ) « ( DER 1 NT ( 2 . 9 ) .PS I DOT > 

EQUIVALENCE (AL(l).UO). (AL(2).CX0). ( AL ( 3 ) . CX AL ) . (AL(4).CX0). (A 
1 L (5 ) » THEO ) » (AL (6 ) .PHIO ) » (AL(7).W0). (AL(S).CZO). ( AL (9 ) . CZAL ) . ( 
2AL (10). CZQ ) ♦ (AL(11 >.CZDE). (AL(12).00). (AL (13). CMO ) , (AL(14)»CMA 
3L ) , (AL ( 15) .CMALD) . ( AL ( 1 6 ) . CMQ ) * ( AL (17). CMDE > . ( AL ( 1 8 ) . VO ) . ( AL ( 
419). CYO). ( AL(20 > ,CY8 ) . ( AL ( 2 1 ) . C YBD ) » ( AL ( 22 ) . CYP ) , ( AL ( 23 ) . CYR ) . 
5 (AL (24 ) .CYOR) . ( AL ( 25 ( ) »P0 ) . (AL ( 26 ) . CLO ) . ( AL (27 ) . CLB ) . (AL(28).C 
6LBD) « ( AL ( 29 ) . CLP ) « ( AL ( 30 ) .CLR ) » ( AL (31). CLDR ) . ( AL ( 32 ) . CLDA ) . (A 
7L ( 33 ) *R0 ) . ( AL ( 34 ) « CNO ) » ( AL (35 > «CNB ) . ( AL ( 36 ) .CNBD ) . (AL(37).CNP) 

8. ( AL ( 38 ) . CNR ) . ( AL ( 39 ) » CNDR ) . ( AL ( 40 ) • CNDA ) . (TABLE ( 1 ). AL ( 1 >) . (T 
9 ABLE ( 50 ) .UMAX ) , ( TABLE ( 51 ). VMAX ) . ( TABLE ( 52 >. WMAX ) . ( TABLE ( 53 ). PMA 
*X)» ( T ABLE ( 54 ) . QM AX ) , ( TABLE (55 ) ,RMAX ) . ( TABLE ( 56 ). THMAX ) » (TABLE! 
*57) . PHM AX) , (TABLE (58) , AXMAX ) . (TABLE (59 ) , AYMAX) . ( TABLE ( 60 ) .AZMAX 
*)» ( TABLE ( 62 ) .DAMAX) , ( TABLE ( 63 ) .DEMAX ) . ( TABLE ( 64 ) .DRMAX ) . (TABLE 
*(65) .UMULT) , ( TABLE ( 66 ) .VMULT). ( TABLE ( 67 ), WMULT ) , ( TABLE ( 68 ) .PMUL 
*T ) » (TABLE ! 69 ) .QMULT). ( TABLE ( 70 ). RMULT ) . ( T ABLE ( 7 1 ) . THMULT ) , (TAB 
®LE ( 72 ) ♦ PHMULT ) , ( TABLE ( 73 ). AXMULT ) , ( TABLE ( 74 ), AYMULT ) . ( T ABLE ( 75 ) 
*»AZMULT)» ( T ABLE ( 76 ) .DAMULT ) « (TABLE (77 ), DEMULT ) , ( TABLE ( 78 ). DRMUL 
*T), ( T ABLE ( 79 ) , U8 IAS) 
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EQUIVALENCE ( T ABLE ( 80 ) ♦ VB I AS )•» ( TABLE (8 1 ) . WB I AS ) . ( TABLE ( 82 ) .PB IAS , 

1). (TABLE ( 83 ) • QB I AS ) • < TABLE (84 ) . RB I AS ) . ( TABLE ( 85 > . THB I AS ) . (TABL 
2E ( 86 ) « PHB I A S ) • < TABLE ( 87 > . AXBI AS > . ( TABLE ( 88 ) . AYB I AS ) . <TABLE(89)« 

3AZBIAS). (TABLE (90 ) « DAB I AS ) . (TABLE ( 91 > . DEB I AS > . ( TABLE (92 > .DRB I AS 

4 ) « (TABLE (98 ).RUN). ( TABLE < 99 ). PASS ) . ( TABLE ( 1 00 ) . A JP ) . ( TABLE ( 1 0 1 

5) «ORM(l)>. (TABLEU 171.AIX) . (TABLE ( 1 1 8 ) . A I Y > . ( TABLE ( 1 1 9 ) . A I Z > . C 

6TABLE ( 1 20 ) * A I XZ ) « ( TABLE ( 1 2 1 1 . WE I GHT ) . < TABLE ( 1 22 ). GRAV ) . (TABLE ( 1 

723). RHO). (TABLE( 124 ) .S) . < TABLE ( 1 25 ) .B ) . ( TABLE ( 1 26 >• C8AR ) . (TABL 

8E( 127) .DEAMPL ) . ( TABLE ( 1 28 ). DEFREQ ) . ( TABLE ( 1 29 ). ALPHAT ) . ( TABLE (1 

930). OT). (TABLE! 131 ) .TT) . ( TABLE ( 1 32 >.TS > . ( TABLE ( 1 33 ) . T I MF ) , (TAB 

SLE( 154 ) .XX) ♦ (TABLE! 155 ) .YX) . (T ABLE < 1 56 ). ZX >.( TABLE (1 57 ). XY > . (T 
SABLE ( 1 58 ) . YY ) . ( TABLE ( 1 59 ). ZY ) . (TABLE ( 160 ). XZ ) . ( TABLE (1 6 1 > . YZ > . 

S (TABLE ( 162) .ZZ) . ( TABLE ( 1 69 ) .CMCON ) . (TABLE < 1 70 ). DALMLT ) . (TABLEU 

S71 ) .uCRTBI ) . (TA8LE( l 72 ) . A I W > . ( TABLE <1 73 ). CTO > . ( TABLE ( 1 74 >, CTTO ) 

S « (TABLE ( 1 75 ) « CTB ) . (TABLE ( 1 76 ). CTBT ) . ( TABLE ( 1 77 ). D ) , (TABLE.U78) 
S.CAPQT). ( TABLE ( 1 79 ) .ELBAR), (T ABLE < 1 80 ) » ELTP ) . (TABLE ( 1 8 1 ) .RB ) . ( 

STABLE < 1 82 ) * PPER ) « ( TABLE ( 1 83 )» TPER ) « ( TABLE ( 1 84 > . BTB I AS ) . (TABLEU 

S85 ) .PPERB I ) . (TABLE ( 186) .TPERBI ) . ( TABLE (1 87 ). DET 1 ) . (TABLE ( 1 88 ) . D, 

SET2 ) 

EQUIVALENCE (TABLE (191) .DRSD ( 1 ) ) ♦ ( I NTEG (51 ) * I NTY (1 > ).» (INTEG(62). 

INOPTS). ( INTEG(63>. INC) • ( I NTEG ( 64 ) . I PR I NT ) . ( INTEG ( 65 ) .NPLOT ) . ( 1 

2NTEG ( 66 > . I READ > * ( I NTEG (67 ) .KSCAN ) . ( INTEG (81 ) . INTX (1 > ) 

ARRAYS DIRECTLY INVOLVED IN ALGORITHM 


AL ( 1 ) 

I NTEG ( I ) 
INTX ( I ) 
INTY ( 1 ) 
INTP ( I ) 
INTV ( I ) 
INTA ( 1 ) 
G(I.J) 
PFU.J) 

PX( I . J) 
PXD(I.J) 

DX ( I ) 

DDX ( I ) 

XBAR ( I ) 
SD(I.J) 

WT{ I , J) 
PXTWT ( I .J) 
PXTDX ( I . J ) 
PXTPX ( I .J) 
DPXTWDX ( I ) 
DDEL A ( I ) 
DELA ( I ) 
WL(IO) 

WL U 1 ) 

WL <1 2 ) 

WL (1 3 ) 

WL (1 4 ) 

WL ( 1 5 ) 

WL (1 6 ) 

WL (1 7 ) 

WL ( 1 8 ) 

WL ( 1 9 ) 

WL (20) 

FSS( 1 )=«T» 
FSS(2)=«T. 


INDICATES 
INDICATES 
INDICATES 
INDICATES 
INDICATES 
INDICATES 
INDICATES 
INDICATES 
INDICATES 
INDICATES 
INDICATES 
FOR 
FOR 


CONTAINS VALUES FOR PARAMETERS ( ALPHAS > 

=1 FOR ACTIVE PARAMETER. =0 FOR INACTIVE PARAMETER 

=1 FOR ACTIVE STATE. =0 FOR INACTIVE STATES 

=1 FOR ACTIVE. =0 FOR INACTIVE ALG. VARIABLE 

ACTIVE PARAMETERS 

ACTIVE STATE VARIABLES 

ACTIVE ALGORITHM VARIABLES 

SENSITIVITY EQUATION MATRIX 

EXPLICIT PARTIALS IN SENSITIVITY EQUATIONS 
SENSITIVITY COEFFICIENTS 
DERIVATIVES OF PX 

DIFFERENCE OF MEASURED AND CALCULATED STATE VARIABLES 

PACKED DX ARRAY 

mean OF MEASUREMENT NOISE 

STANDARD DEVIATION MATRIX 

WEIGHT MATRIX 

INTERMEDIATE CALCULATION 

ACCUMULATED DPXTWDX FROM TIME=0 TO END OF ITERATION 
COVARIANCE MATRIX OF PARAMETERS 
RIGHT HAND SIDE OF PARAMETER CHANGE EQUATION 
PACKED DELA ARRAY 

DELTA ALPHAS FOR PARAMETERS (UPDATES) 

NOPTS .GT. NPTS 
IP .GT. IPAR 

ERROR IN INITIALIZATION OF. STATES 


ABS (THE ) .GT. 
U .LT. 5. FPS 
ABS ( V ) .GT. U 
ABS(ALF) .GT. 
ATTEMPTING TO 


1.5 RADIANS 


1 .5 
USE 


COVARIANCE MATRIX 


RADIANS 

UNFILLED 

SINGULAR 


FLIGHT DATA 
(DETl =0. ) 


WEIGHT MATRIX SINGULAR (DET2=0.) 
IN CRT LOOP 
ACCELERATIONS ON CRT 
CALCOMP OF FLIGHT DATA ONLY 


FSS ( 3 ) *.T. FOR PSEUDO-DATA-FILL 

FSS ( 4 ) s.T. TEMPORARILY FOR VARIABLE DIMENSIONS 
FSS ( 5)*.T. FOR OUTPUT LISTING 
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10 

90003 


30 


FSS ( 7>*.T. 
FSS ( 8)*.T. 
FSS ( 9>«VT* 
FSS <jO) = .T« 
FSSIl 1 )*.T. 
FSS ( 1 2 ) * .T. 
FSS<l3)*.T. 
FSS< 14) = .T. 
FSS(l5)*.T. 
FSS < 1 6 ) = .T • 


FOR LATERAL STATES ON CRT 
FOR CONTROLS ON CRT 

FOR PSEUDO FLIGHT FUNCTIONS ' ' ' ' " 

TO SKIP UPDATE OF ALPHAS 
TO ENABLE TYPEWRITER 

TO RETAIN STATE I.C. S DURING STORING OF AL IN ALG 
TO EXIT CRT LOOP 
TEMP. FOR REPLOT 
TO ENTER CRT LOOP 
FOR TABLES ADDRESSING 
FSS( 1 .7.8 >=.F. FOR LONGITUDINAL STATES ON CRT 
LOG I C ( 1 ) = »T. TO CALCULATE WEIGHT UPDATES 
LOG I C 1 2 ) *.T. FOR DIAGONALIZED WEIGHT MATRIX 
LOGIC ( 3 )*.T. FOR MAIN PROP = FPPER ( 1 ) 

LOG I C < 4 ) = . T • FOR TAIL * FTPER ( 1 I 
FOR CTT »F ( BET AT I 
FOR CZOE=CMDE*CMCON 

FOR LONGITUDINAL STATES=FL I GHT DATA 
FOR ALFT = AT AN ( WO « UO ) 

FOR TRIM CONDITIONS 
FOR AUTOMATIC TRIM 


LOGIC < 5 ) * • T • 
LOG IC<6)*.T. 
LOGIC ( 7 > * .T ♦ 
LOGIC ( 8 ) * • T* 
LOGIC <9 ) = .T. 
LOGIC ( 10)*. T. 


LOGIC < 2 1 >*«T« FOR AUTOMATIC 2 PASS SYSTEM 
I PR I NT = 1 TO PRINT ICS SET 

I PR I NT *2 TO REWIND MF FILE SET 

I PRINT *3 TO OBTAIN RANDOM NUMBERS SET 

IPRINT*4 FOR CALCOMP PLOT (SET VALUE BEFORE 
ENTERING PRINT - USES FSSIl. 2.7. 8. 

I PR I NT =5 TO RETURN (ROUTES OUTPUT) 

I READ 3 1 TO STORE ALG IN Al.ETC. 

I READ*2 TO STORE AL IN ALG. ETC. 

IREAO*3 TO GET NOMINAL CONDITIONS 
IREAD*4 TO GET LONGITUDINAL COND I T 
IREAO=5 TO GET LATERAL CONDITIONS 
IREAD-6 TO FILL FLIGHT DATA ARRAYS 
CALL LDRSEC 

IF .( ISKIP.GT.O) GO TO 10 
I SK I P*1 

CALL CYCLE (90006S) 

ASSIGN 90001 TO NOPER 

90002 TO NHOLO 

90003 TO NRESET 
9000 A TO NTERM 
9001 A TO NPRINT 
90015 TO NREAD 


VALUE BEFORE PRINT 
VALUE BEFORE PRINT 
VALUE 8EF0RE PRINT 
EXITING CRT LOOP AND 
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IONS SET 
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BEFORE 

READ 

E 

0151 

SET 

VALUE 

BEFORE 

read 
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SET 

VALUE 

BEFORE 

read 
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40 


ASSIGN 
ASSIGN 
ASSIGN 
ASSIGN 
ASSIGN 

CALL XDSPLAY 
CALL REAOY 
CONTINUE 

IF (NOPTS.LE.NPTS) 
WL 1 1 0 ) = • T . 

GO TO 740 

IF ( ,N0T.FSS(4 ) ) GO 
DO 40 1*1 . I PAP 
DPARAM( I ) * 1 OH 
I NTP ( I ) =0 
IP*0 


(LDISI .LDISO.VARCHNG, I TYPE . I VARBUF .FSS ( 1 6 ) ) 


GO TO 30 


TO 110 
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0174 


DO 50 1*1 «NPAR 

IF ( I NT EG (I).EQ.O) GO TO 50 

IP.IP+1 

I NTP ( I P ) * I 

DPARAM( IP)*PARAM( I ) 


E 0175 
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60 



50 


60 
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1 00 
1 10 


120 


1 30 


140 

150 

160 
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190 


CONTINUE 

IF (IP.LE.IPAR) GO TO 60 
WL< 1 1 )*.T. 

GO TO 740 
DO 70 1*1.8 

IF ( 1 NTY ( I ) *EQ. 0 ) GO TO 70 
IF { I NTX ( I ) «EQ. 1 ) GO TO 70 
WL(12)*.T. 

GO TO 740 
CONTINUE 
I A 1 *0 
I V*0 

DO 80 I =1 .8 
DVAR ( I ) * 1 OH 
I NTV ( 1 1 *0 

IF ( INTY( I) .EQ. 1 ) IA1-IA1+1 

IF ( I NTX (!) . EQ.O ) GO TO 80 
I VMV+1 

DVAR( IV >*ALV( I ) 

I NTV ( IV)*l 
CONTINUE 
!A = 0 

DO 90 I ”1 . 1 1 
DALG ( I )*10H 
INTA (I 1 *0 

IF ClNTYm.EQ.O) GO TO 90 
IA*IA+1 

DALG ( I A ) *AL V ( I ) 

I NT A ( 1 A ) = I 
CONTINUE 
I A2* I A- I A 1 

I TRNMLT (2 ) a I TRNMLT (41*1 MULT A (31* I NVWT ( 2 > * I N VWT ( 3 > * l A 

I TRNMLT (3 )* I MULT A (2 ) = l NVSEN < 2 ) = l NVSEN ( 3 > = I MULTB ( 2 ) = I MULTB (3 ) * IP 

IF (. NOT. LOGIC! 1 ) ) GO TO 110 

DO 100 1*1. I A 

DO 100 J=1,IA 

WT(I , J)=0. 

IF (l.EQ.J) WT ( I « J ) * 1 . 

CONTINUE 
NEQ*9+ I P* I V 
IF (FSS < 3 ) ) NEQ=9 
DO 1?0 l«h|A 
XBAR ( I ) *0. 

DO 120 J=1 . IA 
SD ( I . J 1 “O. 

DO 130 1*1. IP 
PXTDX ( I )*0. 

DO 130 J=1 . IP 
PXTPX <I«J)*0. 

<«9 

DO 230 1*1 . IV 
I I*!nTV<I ) 

DO 230 J*1 . IP 
JJ*INTP(J) 

K*K+1 

DERINT( 1 .K)*0. 

GO TO (140.150.160.170.180.190.200.210). II 

IF ( JJ.EO.l ) 220,230 

IF (JJ.EQ.18) 220.230 

IF (JJ.EQ.7) 220,230 

IF (JJ.E0.25) 220,230 

IF (JJ.E0.12) 220.230 

IF (JJ.E0.33) 220.230 
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200 IF (JJ.EO.5) 220.230 

210 IF (JJ.EQ.6) 220.230 

220 DERINT< 1 .K>=1 • * ' “ • ' 

230 CONTINUE 

DO 240 J=1.NPAR 
OELA (J ) =0 • 

DO 240 1*1.8 
IF (I.LT.7) PXD(I.J)=0. 

IF (J.LT.9) G ( I . J ) =0 . 

240 PF ( I . J 1 *0 • 

G(8«4)=l. 

INT = 0 
H*DT 

VARCHNG-. FALSE. 

T *0. 

U*UO 

v*vo 

w=wo 

P=PO 

0*00 

R*RO 

the*theo 

phi=phio 

psi*o. 

C YCL*2 . *P I /DEFREO 
ALFT*ALPHAT*RA0O 

IF <LOGIC(8n ALFT=ATAN2 (WO .UO) 
AJ*0. 

I TS=— 1 

NOEL*0 

DELX=TT*INC 

NOITSPS=DELX/'DT 

TMAXX*DELX*(NOPTS-l > 

PTSINV-1 ./NOPTS 

GINV*1./GRAV 

AM AS INV*GRAV/WE I GHT 

AIYlNVal ,/AIY 

BTWO».S*B 

CBAP?*.5*CBAR 

RHOS2=RHO*S*.5 

A I =RHOS2*AMASINV 

A2=A 1 *CBAR2 

A3-A1 *BTWO 

A4=RHOS2*B 

A5=A4*B TWO 

CONI N* 1 »/ ( A I X*A I Z-A I XZ**2 ) 

A 6* A 1 Z*CONI N 
A7*(AIZ-AIX)»AIYINV 
A8*RHOS2*CBAR*A I Y I NV 
A9*A8*CBAR2 

VDTCON* 1 » / ( 1 •— A3*CYBD ) 

PF21 =A3*VDTC0N 
A9CMAL0*A9#CMALD 
81 *A I XZ*CONI N 
B2* A I Y— A I Z 
B3*AIX-AIY 
B4=A IXZ*AI YINV 
B5*AJX#C0NtN 

PF4 1 * A5* ( A6*CLBD+B 1 *CNBD ) 

PF61 =A5* (B1 *CLBD+B5*CNBD ) 

842=2. *B4 
Cl W*C0S ( A I W ) 
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SIW=SlN(A IW> 

CLDAP=CLDA*C I W-CNDA*S I W 
CNOAP=CLDA*S I W+CNDA*C I W 
C0N=2 • »RH0*D**4 
CONI *2.*C0N 
C0N2=EL8AR*CTB*C0N 
CON3=Ct W*CONl 
CON4=SIW*CON! 

CON5=SIW*CON2 
CON6=CIW*CON2 
C0N7«RB*C0N1 
CON8=RHO*C AP0T**4 
C0N9=ELTP*C0N8 
DCON»8.*D*D/Pl 
90006 CONTINUE 
I TS= I TS+ 1 
90002 CONTINUE 
90008 CONTINUE 
K*9 

DO 260 1=1 . IV 
DO 260 J=1 .IP 
K*K+1 

260 PX < I « J ) *DER I NT < 1 « K ) 

IP < .N0T.FSSI9) > GO TO 270 

BETAT^DAeDB-DEsDELB^DR^PPER^TPERsO . 

IF (T.LE.CYCL) DE=DEAMPL*S I N <DEFREQ*T ) 

GO TO 280 
270 FT = T 

IS=FT/DELX 
SS=FT/DELX-IS 
IS1=IS+1 
I S2= IS+2 

IF ( IS2.GT.NOPTS) !S2=NOPTS 
SSI =1 .-SS 

DA=SS1 *FDA ( IS1 )+SS*FOA ( ! S2 ) 

DE=SS1*FDE( I SI )+SS*FOE< IS2) 

DR=SS1*FDR( IS1 >+SS*FDR< IS2) 

BETAT“SS1 «FBT< IS1 >+SS*FBT< I S2 ) 

DB=SS1»FDB( I SI )+SS*FDB( I S2 ) 

DEI_B = SS1*FDELB ( I SI >+SS*FDELB C I S2 ) 

PPER=FPPER ( 1 ) 

IF (.NOT.LOGICO) 1 PPER = SS 1 *FPPER < I S 1 >+SS#FPPER < I S2 ) 
TPFR=FTPER( 1 ) 

IF ( .NOT .LOG I C < 4 ) ) TPER = SS1*FTPER( IS! > + SS*FTPER< I S2 ) 
IF ( .NOT. LOGIC (7 ) ) GO TO 280 
U=SS 1 *UM ( I SI J + SS*UM< IS2) 

WaSSl *WM ( I SI >+SS*WM ( IS2 ) 

Q*SS 1 *QM ( I S 1 >+SS*OM( IS2) 

THE=SSl*THM< 1SI ) + SS*THM< 1 S2 ) 

280 CONTINUE 

IF ( ABS (THE l.LT.l .3 ) GO TO 290 
WL<1?)=.T. 

GO TO 740 
290 CONTINUE 

IF (u.GT.5.) GO TO 300 
WL < 1 4 ) = • T . 

GO TO 740 
300 CONTINUE 

IF (ABS(V).LT.U) GO TO 310 
WL(13)=»T. 

GO TO 740 
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31 0 


320 


330 
34 0 


UINVbI./TJ 

ALF=ATAN2 (W.U)-ALFT 

IF (ABS1ALF > .LT. 1.3) GO TO 320 

WL(16)*.T, 

GO TO 740 

U2=U*U 

V2=V*V 

W2=y*W 

VR2*u2+V2+W2 
VR*SORT < VR2 ) 

VR I NV* 1 ./VR 

BET = AS I N < V*VR I NV ) 

ETA* 1 232.4PPER/60 . 

ETAT*40.*TPER 
CT=CTO+CTB*DELB 
IF ( .NOT .LOG IC ( 5 ) ) GO TO 330 
FBETsABS <8ETAT*DRAD > 

IF (FBET.GT.22. ) FBET-22. 

IS*F8ET*.5 

SS*FBET*.5-IS 

TDUM* ( 1 SS ) *FCT ( IS + 1 >+SS*FCT( 1S + 2) 

IF (BETAT.LT. 0. ) TDUM*— TDUM 

CTT=TDUM+CTTO 

GO TO 340 

CONTINUE 

ctt=ctbt*betat+ctto 

CONTINUE 

TEMPnET A**2 

TEMPI *CT*TEMP 

TEMP2*DB*TEMP 

TEMP3»ETAT**2*CTT 

TX=CON3*TEMPI 

T7=-CON4*TEMPl -CONB*TEMP3 

AMX=CON5*TEMP2 

AMY*-CON7*TEMPl -CON9*TEMP3 

AMZ=CON6*TEMP2 

VSS2=OCON* TEMPI 

VSS*SORT ( VSS2 ) 

VSR*VSS+VR 
VSRVR* VSR*VR I NV 
VSR2*VSR**2 
STHE*S I N (THE ) 

CTHE*COS (THE ) 

SPHI »S I N (PH I ) 

CPH!«C0S(PH! ) 

gcthe*grav*cthe 

GSTHF*GRAV*STHE 

TTHEsTAN(THE) 

TXM<=AMASINV*TX 

TYM*AMASINV*TY 

TZM=AMASINV*TZ 

AMY I Y* A I Y I NV*AMY 

A1U*A1*U 

A 1 V= A 1 *V 

A I A 1 

A4U=A4*U 

A4V= A4*V 

A4W*A4»W 

A8U*A8»U 

A8V* A8*V 

A8W=A8*W 
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A 1 VR»A 1 *VR 

A2VRoA2#VR 

A3VRa A3*VR 

A4VRaA4*VR 

ASVR=AS*VR 

A9vRaA9*VR 

A1 VR2*A 1 *VR2 

A4VR?*A4*VR2 

A4VSR2=A4*VSR2 

A8VR2*A8*VR2 

B2VRaBTW0*VRINV 

PF22aPF21*VR 

PF42aA6*A5VR 

PF43=A4VSR2*0A 

PF44aBl*A5VR 

PF62*B3*A5VR 

CBAR2VR*CBAR2*VR I NV 

P2=P*P 

Q2aQ*Q 

R2=R*R 

PO«=P*Q 

PR=P*R 

QR=Q*R 

Q1=Q*AIXZ 

CXI =cxo+alf*cxal 

CX2=CXQ*Q*CBAR2VR 

CY1 =CYO+BET*CY0+DR*CYDR 

CY2«B2VR*(P*CYP+R*CYR) 

czi =czo+alf*czal+de*czde 

CZ2=c2Q*0*C8AR2VR 
FXM= A 1 VR2* i CXI +CX2 ) 

FYM=A 1 VR2* ( CY1 +CY2 ) 

FZMr A 1 VR2* ( CZI +CZ2 ) 
UDP=R*V-0*W-GSTHE 

vdp=p*w-r*u+gcthe*sphi 

WDP=Q*U-P*V+GCTHE*CPH! 

UDOToUDP+FXM+TXM 

VDOT =VDTCON*< VDP+FYM+TYM ) 

WDOT=WDP+FZM+TZM 

IF { I NTX ( 1 ) »EQ« 0 ) UDOT *0 • 

IF ( INTX(2 ) .EQ.O ) VD0T=0. 

IF ( INTX<3> .EQ.O) WD0T=0. 

IF ( .NOT. LOGIC (9) ) GO TO 350 

CALL HALT 

CXI =CX1 -CXO 

CY1 =CYl-CYO 

CZI =CZ! -CZO 

TMP1 a I ./A I VR2 

CXO=CXO-UDOT*TMPl 

CYO=CYO-VOOT*TMPI /VDTCON 

CZO=cZO-WDOT*TmP1 

UDOT = 0 • 

VDOT = 0 . 

WDOTaO. 

CXI =CXI +CXO 
CY1*CYI+CY0 
CZI=CZ1+CZ0 
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350 CONTINUE 

A LFO=WOOT»UINV = 

BETD= VD0T*VR I NV 

BDB2vR=BETD*B2VR 

CL1 =CLO+BET*CLB+DR*CLDR 

CL2 = ( BETD*CLBD+P*CLP+R*CLR >*B2VR 

CL3=OA*CLDAP 

CM1 =CMO+ALF*CMAL*DE*CMDE 

CM2 = (ALF0*CMALD+0*CMQ)*CBAR2VR 

CN1 =CNO+BET*CNB+DR*CNDR 

CN2= (BETD*CNBD+P*CNP+R*CNR )*B2VR 

CN3=0A*CNDAP 

LIX=A4VR2*(CL1+CL2)+A4VSR2*CL3 
M I Y* A8VR2* ( CM 1 +CM2 ) 

N I Z« A4VR2* ( CN1 +CN2 ) -A4 VSR2*CN3 
F4P-Q* ( B2*R+P*A I XZ ) +L I X+AMX 
F6P=Q* ( B3*P-R*A I XZ ) +NI Z+AMZ 
PDOT=A6*F4P+Bl *F6P 
QD0T=A7*PR+B4* (R2-P2 ) +M I Y+AM Y I Y 
ROOT =B 1 *F4P+B5*F6P 
IF ( INTXI4 ) *EQ.O) POOT = 0 • 

IF ( INTX(5) .EQ.O) QDOT = 0 • 

IF ( INTX(6) .EQ.O) ROOT=0. 

IF ( .NOT. LOGIC 19) ) GO TO 360 
LOG I C < 9 ) = .F . 

CL 1 =CL I — CLO 
CMl*CMI-CMO 
CN1 *CN1 -CNO 
TMP4*1 •/’A4VR2 
CLO«CLO-F4P*TMP4 
CMO=CMO-QDOT/’A8VR2 
CNO=CNO-F6P*TMP4 
POOT = 0 . 

QDOT = 0. 

ROOT = 0 • 

CL.1-CL1+CLO 
CM 1 =CM 1 +CMO 

CNI«CNl+CNO ^ 

CALL READY 
360 CONTINUE 

AXCG*G I NV* ( UDOT-UDP ) 

AYCGeG I NV* ( VDOT-VOP ) 

AZCG«GINV*(WDOT-WDP) 

AX I * AXCG+GI NV* ( — ( R2 + Q2 ) *XX+ < PQ-ROOT )* YX+ (PR+QDOT )*ZX ) 

AY I = AYCG+G I NV* { (PQ+RDOT ) *X Y- (P2+R2 > * YY+ ( QR-PDOT ) *ZY ) 

AZ I = AZCG+G I NV* ( ( PR-QDOT ) *XZ + < QR+PDOT ) *YZ- (P2+Q2 )*ZZ ) 

THEDOT«Q*CPH I -R*SPH I 

PS I DOT * <Q*SPH!+R*CPHI I/CTHE 

PH I DOT *P+PS IDOT*STHE 

IF ( INTX<7) .EQ.O) THEDOT=0. 

IF ( INTX(8>«EQ.O) PHIDOT=0. 

IF ( .NOT .FSS ( 3 ) ) GO TO 370 

IF (MOD(ITS»NOlTSPS).NE.O) GO TO 530 

IF (.NOT.LOISI (IB)) NOEL*NOEL+ 1 

UM (NOEL )=U 

VM<NOEL)*V 

WM(NOEL)«W 

PM (NOEL )*P 

QMINOEL )*Q 

RM<NOEL)»R 

THM(NOEL)*THE 

PHMINOEL )=PHI 
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AXM(nOEL)=AXI 
AYM(nOEL)=AYI 
AZM(NOEL)*AZI 
FDA(nOEL)*DA 
FDE(NOEL)=DE 
FDR ( NOEL ) *OR 
FBT ( NO EL ) “BET AT 
FOB (NOEL >*DB 
FDELB < NOEL > =DELB 
FPPER ( NOEL 1 =PPER 
FTPER < NOEL ) =TPER 
GO TO 330 
370 CONTINUE 

PF( 1 ,2> *A1 VR2 

PF< 1 ,3>*A1VR2*ALF 

PF( 1 ,4 )*A2VR*Q 

PF ( 2 , 1 9 ) = VDTCON* A 1 VR2 

PF < 2 , 20 ) =F2 1 9*BET 

PF (2*21 )«PF21*VOOT 

PF(2,22 >=PF22*P 

PF(2.23 >=PF22*R 

PF (2,24 ) S F2 1 9*DR 

PF(3,6)=A1VR2 

PF (3,9) ®F 1 3 

PF (3,10) *F! 4 

PF (3,11 )=A1 VR2*DE 

PF (4,19) “F2 1 9*PF4 1 

PF ( 4 ,20 )=F220*PF41 

PF (4,21 >«F221*PF41 

PF ( 4 , 22 ) «F222*PF4 1 

PF ( 4 , 23 ) *F223*PF4 1 

PF ( 4 , 24 > =F224*PF4 1 

PF(4,26)«A6*A4VR2 

PF(4,27 >=F426*BET 

PF (4 ,28 )*PF42*BETD 

PF(4,29>»PF42*P 

PF(4,30 )*PF42*R 

PF (4,31 )*F426*DR 

PF (4 , 32 )»PF43*(A6*CIW-Bl*SIW) 

PF ( 4 , 34 )=BI *A4VR2 

PF(4,35 )«F434*BET 

PF (4 , 36 ) *PF44*BETD 

PF (4,37 1 *PF 44 *P 

PF(4,38 >=PF44#R 

PF ( 4 . 39 ) «F4 34*DR 

PF (4,40 )=-PF43*(Bl*CIW+A6*SIW) 

PF(5,8)*A1 VR2*A9CMALD 

PF ( 5 , 9 ) =F39*A9CMALD 

PF ( 5 » 10 )«F310*A9CMALD 

PF (5,1 1 ) S F3 1 1*A9CMALD 

PF (5 , I3)=A8VR2 

PF ( 5 , 14 >=A8VR2*ALF 

PF ( 5 , I 5 ) *A9VR*ALFD 

PF (5,16 )=A9VR*Q 

PF (5 . 17>«A8VR2*DE 

PF (6 , 1 9 )*F219#PF61 

PF ( 6 , 20 ) *F220*PF6 1 

PF (6,21 ) “F22 1 *PF6 1 

PF ( 6 , 22 > =F222*PF6 1 

PF ( 6 , 23 ) =F2 23*PF6 1 

PF ( 6 , 24 ) =F224*PF6 I 


E 0546 
E 0347 
E 0348 
E 0549 
E 0350 
E 0551 
E 0532 
E 0553 
E 0554 
E 0555 
E 0536 
E 0557 
E 0558 
E 0559 
E 0560 
E 0561 
E 0562 
E 0563 
E 0564 
E 0563 
E 0566 
E 0567 
E 0568 
E 0369 
E 0570 
E 0371 
E 0572 
E 0373 
E 0574 
E 0575 
E 0576 
E 0577 
E 0578 
E 0579 
E 0580 
E 0381 
E 0582 
E 0583 
E 0584 
E 0585 
E 0586 
E 0587 
E 0588 
E 0589 
E 0590 
E 0591 
E 0592 
E 0593 
E 0594 
E 0593 
E 0596 
E 0597 
E 0598 
E 0599 
E 0600 
E 0601 
E 0602 
E 0603 
E 0604 
E 0603 
E 0606 


67 



PF (6.26) “F4 34 
PF(6.27)«F435 
PF (6.28 >=F436 
PF(6.29 )=F437 
PF (6 . 30 )«F438 
PF (6.31 >=F439 

PF (6. 32 >«PP43*<B1*C I W-85*S I W) 

PF (6.34) ®B5*A4 VR2 

PF ( 6 . 35 ) =F634*8ET 

PF ( 6 • 36 J »PF62*BET0 

PF(6.37)=PF62*P 

PF ( 6 « 38 ) “PF 62*R 

PF(6.39)*F634*0R 

PF (6.40 )»-PF43*(8S*CIW+81*SIWJ 

G1P=2.*CX1+CX2 

G2P® 2»*CYI+CV2-BET*CYB 

G3Pj»2»*C2) +CZ2 

G4P=2« * (CL 1 +CL3*VSRVR )+CL2-BET*CLB— BD82VR*CLBD 
G5P*2«»CMI+CM2 

G6P*2»*(CN1 -CN3*VSRVR)+CN2-BET*CNB-BDB2VR*CNB0 
G< 1 .1 )«A1U*G1P— At W*CXAl_ 

G ( 1 « 2 ) * A 1 V*G 1 P+R 

G ( t .3)*AlW*GlP+AtU*CXAL-Q 

G(1 . *5 ) a A2VR*CXQ— W 

Sit »6)=V 

G ( ) .7>=-GCTHE 

G ( 2 . ] )»VOTCON*( A)U*G2P-R) 

G ( 2 ♦ 2 ) * VDTCON* ( A 1 V»G2P+A 1 VR*CYB ) 

G < 2. 3 ) ■ VDTCON* ( A 1 W»G2P+P ) 

G (2.4 > * VDTCON* ( A3VR*CYP+W ) 

G ( 2 . 6 >= VDTCON* { A3VR*CYR-U ) 

G ( 2 . 7 ) a — GSTHE*SPH I *VDTC0N 

G ( 2 . 8 ) *GCTHE*CPH t *VDTC0N 

G ( 3. 1 ) “ At U*G3P— A 1 W*CZAL+Q 

G(3.2)=At V*G3P-P 

G ( 3 . 3 > = A 1 W*G3P+A 1 U*CZAL 

G (3.4 ) =— V 

G(3.3 >«A2VR*CZQ+U 

G ( 3 . 7 ) a -GSTHE*CPH I 

G ( 3. 8 ) =-GCTH£*SPH I 

G4 1 a A4U*G4P 

G42a a4V*G4P+A4VR*CLB 

G43=A4W*G4P 

G44=0I +A5VR*CLP 

G45aP*B2+P*AtXZ 

G46=0*B2+A5VR*CLR 

G61 =A4U*G6P 

G62=AAV*G6P+A4VR*CNB 

G63aA4w*G6P 

G64=Q*B3+A5VR*CNP 

G65*P*B3-R*AIXZ 

G66= a5VR*CNR— Q I 

G(4, 1 ) a A6*G4 1 +B1 *G6 1 +PF4 1 *G2 1 
G ( 4 . 2 ) a A6*G42+B1 *G62+PF4 1 *G22 
G ( 4 , 3 ) a A6*G43+B t *G63+PF4 1 *G23 
G (4.4 ) =A6*G44+B1 *G64+PF4 1 *G24 
G (4 .8 ) = A6*G45+B1 *G65 
G ( 4 .6 ) a A6*G46+Bt *G66+PF41 *G26 
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G(4,7)*PF41*G27 
G ( 4 . 8 ) = PF4 1 *G28 

G ( 5. 1 )=A8U*G5P-A8W*CMAL+A9CMALD* (G31 

G ( 5. 2 ) * A8V*G5P+A9CMAlD*G32 

G(5.3 )=A8W*G5P+A9CMALD*G33+A8U*CMAL 

G(5.4 >=R*A7-B42*P+A9CMALD*G34 

G ( 5 . 5 > = A9VR*CMQ+ A9CMALD*G35 

G(5.6)=P*A7+B42*R 

G (5.7 ) = A9CMALD*G37 

G(5.8)*A9CMAL0*G38 

G (6. 1 )aBl +G41 +S5#G6 1 +PF61 *G21 

G (6.2 ) *BI *G42+B5*G62+PF6 I *G22 

G ( 6. 3 ) *B1 *G43+B5#G63+PF61 »G23 

G ( 6. 4 ) =B1 *G44+B5*G64+PF6 1 *G24 

G(6.5)»B1*G45+B5*G65 

G(6.6)=B1*G46+B5#G66+PF61*G26 

6(6.7) =PF61 *G27 

G (6.8 ) =PF61 *G28 

G(7.5)=CPHI 

G ( 7. 6 ) *-SPH I 

G ( 7. 8 )*-PS I DOT#CTHE 

G ( 8. 5 ) =SPHI *TTHE 

G (8. 6 ) = CPH I *TTHE 

G ( 8, 7 ) »PS IOOT/CTHE 

G ( 8. 8 ) =THEOOT*TTHE 

L*9 

DO 390 1*1 , I V 

11= INTV ( I ) 

DO 390 J=1 . IP 
L=L+1 

JJ=1NTP( J) 

GPX=0. 

DO 380 K*1 , IV 
KK=tNTV(<) 

GPX = GPX+G( I I .KK )*PX (KiJ) 

DERInT ( 2 .L ) *GPX+PF ( I I.JJ) 

IF I1I.LT.7) PXD (II. J J ) =DER I NT ( 2 .L > 
CONTINUE 

IF ( I A 1 »EQ • I A ) GO TO 450 
G(1 .1 )*0. 

G ( 1 . 2 ) *-R 
G( 1.3)*Q 

G ( 1 .4 ) = YX*0+ZX*R 
G ( 1 .5)=-2.*XX*Q+YX*P+W 
G ( 1 .6)=-2.*XX*R+ZX*P-V 
G ( 1 .7>=GCTHE 
G (2. 1 ) *R 
G (2.2 )*0« 

G ( 2. 3 ) ■— P 

G(2.4)=-W+XY*Q-2.*YY*P 

G (2.5 ) = XY*P+ZY*R 

G ( 2 . 6 ) *U-2. *YY*R+ZY*Q 

G ( 2 . 7 ) =GSTHE*SPH I 

G ( 2. 8 ) = — GCTHE*CPH I 

G( 3. 1 )=-Q 

G(3.2)=P 

G(3.3)»0. 

G(3.4 )=V7XZ*R-2.*ZZ*P 
G ( 3 . 8 ) = -U+ YZ»R-2 . *ZZ*0 
G ( 3 « 6 ) * XZ*P+ YZ*Q 
G ( 3. 7 ) =GSTHE*CPHI 
G ( 3. 8 ) =GCTHE*SPH I 
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L = I V 
M» I A 1 

DO 44 O 1 = 1 • IA2 
L=L+ 1 
M = M+1 

11 = 1 NT A (M ) —8 
DO 44 O J=1 » IP 
JJ=INTPCJ> 

gpx=o« 

DO 400 K*1 . TV 
KK=INTV(t<) 

400 GPX=GPX+G< I I «KK )*PX(K« J> 

GO TO (410«420«430) « II 

410 PX(U.J)=(GPX+PXD(1 < JJ)+ZX*PXD(5« JJ )- YX*PXD { 6 • JJ) >*GINV 
GO TO 440 

420 PX(l_« J) = (GPX+PXD(2« JJ ) -ZY*PXD (4 • JJ )+XY*PXD (6 * JJ) >*GINV 

GO TO 440 

430 PX(L« J)* (GPX+PXD (3« JJ ) + YZ4PXD (4 « JJ )— XZ*PXO (5 * JJ ) )*GtNV 

440 CONTINUE 

G(1.4)=G(2*5)=G(3.6>=0. 

450 CONTINUE 

IF (IA1.E0.IV) GO TO 480 
J=0 
L = 0 

DO 470 1=1 t 1 1 

IF ( I »LE*8 ) J=J+INTX(I) 

IF ( INTYI IJ «EQ.O) GO TO 470 
L*L+ 1 

IF (I.GT.8) J=J+1 
DO 460 K=1 . IP 
460 PX(l_«K)=PX( J.K) 

470 CONTINUE 
480 CONTINUE 

IF <MOD< ITS.NOITSPSJ.NE.O) GO TO 530 
IF (.NOT.LDISI (18) ) N0EL=N0EL+1 
IF (UM(NOEL).LT.l .E5) GO TO 490 
WL ( 1 7)=»T» 

GO TO 740 
490 CONTINUE 

DX < 1 )=UM (NOEL )— U 
DX ( 2 ) =VM (NOEL )— V 
DX ( 3 ) *WM ( NOEL ) — W 
DX ( 4 ) =PM ( NOEL ) -P 
DX ( 5 ) *QM (NOEL )— Q 
DX ( 6 ) =RM ( NOEL ) -R 
DX ( 7 ) =THM (NOEL ) —THE 
DX ( 8 ) =PHM ( NOEL ) -PH I 
DX (9 ) =AXM (NOEL )— AX I 
DX ( 1 0 ) = AYM (NOEL )-AY I 
DX ( 1 I ) =AZM (NOEL )-AZ I 
DO 500 1 = 1 . IA 
I I = I NT A ( I ) 

DDX (!) =DX (II) 

A J= AJ+DDX ( I )**2 

500 XBAR ( I )=XBAR( I )+DDX( I )*PTSINV 
DO 510 1 = 1 * IA 

DO 51 0 J=1 * IA 

510 SO ( I t J ) =SD ( I « J )+DDX ( I )*DDX( J >*PTSINV 

C PXTWT= (PX)T(WT) 

CALL MASCNT ( I TRNMLT »PX« WT«PXTWT ) 
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c DPXTWDX= (PXTWTKDOX) 

CALL MASCNT ( I MULTA .PXTWT . DDX. DPXTWDX > 

DO 520 I ■ 1 • I P 

PXTDX ( I ) *PXTDX ( I ) +DPXTWDX ( I ) 

DO 520 J=I . IP 
DO 520 K*1 • I A 

520 PXTPX * I . J)=PXTPX( I « J I+PXTWT ( I iK )*PX (K« J) 

530 CONTINUE 

TABLE * 1 34 )=T 
TABLE <1 35 )=U 
TABLE* 1 36 )=V 
TABLE* 137) 

TABLE C 1 38 ) «P 
T ABLF * I 39 ) *Q 
TABLE* 140 )*R 
TABLE* 141 >*THE 
TABLE* 142) =PH I 
TABLE* 143) *PS I 
TABLE* 144 >»UDOT 
TABLE* 145) = VDOT 
TABLE* 146)»W00T 
TABLE* 1 47 ) =PD0T 
TABLE* 1 48 ) =QD0T 
TABLE *1 49 )*RDOT 
TABLE* 150)- THEDOT 
TABLE* 151 )*PHID0T 
TABLE* 152)=PSID0T 

C DACS FOR TIME HISTORY RECORDERS 

DAC * 1 ) = (U-UCRTB I ) /UMAX 
DAC ( 2 ) =V/VMAX 
DAC * 3 ) * W/WMAX 
DAC(4)»P/PMAX 
DAC (5 )=Q/QMAX 
0AC*6)=R/RMAX 
DAC < 7 ) =THE/THMAX 
DA C ( 8 ) = PH I /PHM A X 

C RITECRT PLOTS IN REAL TIME IF T »LE» TMAXX 

IF (. NOT. FSS<3). AND. (T.LE. TMAXX) ) CALL RITECRT 
1 AGE ) 

IF (LD I S I * 22 ) ) CALL SCANNER *KSCAN) 

CALL OSPLAY 

IF (LOISI (17) ) GO TO 90050 
IF (VARCHNG.AND.FSS* 1 1 ) ) CALL TYPEVAR 
IF (FSS * 1 1 ) .AND.LDIS I * 14 ) ) CALL TYPEVAR 
IF (FSS *15)) GO TO 770 
90050 CALL RTMODE 
90001 CONTINUE 

CALL RECORD 

IF (NOEL.GE.NOPTS) GO TO 560 
90005 CONTINUE 

CALL IGRATE1 
I NT * 1 

CALL RECYCLE 

560 IF ( FSS ( 3 ) ) GO TO 770 x 
CALL HALT 

IF (FSS (10)) GO TO 690 

IF ( .NOT.LOGtC( 1 1 ) ) GO TO 570 

IEVEN»!EVEN+1 

IF ( MOD ( I EVEN. 2 ) *EQ« 0 ) GO TO 690 
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570 CONTINUE 

■■ PASS=PASS+1 • 

DO 580 J*1 » IP 
DO 580 I=J« IP 

580 PXTPX ( I • J ) =PXTPX < J « I ) 

C PXTPX® (PXTPX) INV 

CALL MASCNT { INVSEN.PXTPX < DET1 . AA 1 ) 

IF (DET1.NE.0. ) GO TO 590 
WL(18)=.T. 

GO TO 750 
590 AUP=AJ 

C DDELA“ (PXTPX) (PXTDX ) 

CALL MASCNT { I MULTB .PXTPX « PXTDX • DDELA ) 

DO 600 I®1 « IP 
I I = I NTP ( I ) 

600 DELA (II) “DDELA ( I ) *DALMLT 

C FSS (5 ) ON TO WRITE AL AND DELA ON TAPE 50 

IF { .NOT «FSS(5) ) GO TO 670 
WRITE (MF.960) RUN, PASS* DET1 • DET2 . A JP 
WRITE (MF.970) DALG ( 1 ) . ( ( COM , DALG ( l ) ) , I =2 , I A > 

WRITE ( MF * 980 ) 

WRITE ( MF ,990 ) (XBAR ( I ) , I » 1 , I A ) 

WRITE (MF.1000) 

DO 61 0 1 = 1 , I A 

WRITE ( MF ,990 ) (SD( I , J ) , J= 1 , IA ) 

610 CONTINUE 

WRITE ( MF ,1010) 

DO 620 I»1 , IA 

WRITE ( MF .990 ) (WT ( I , J ) , J= 1 • I A ) 

620 CONTINUE 

WRITE ( MF , 1 020 ) ( (PARAM ( I ) , AL ( I > .DELA ( 1 ) ) . 1 = 1 .NPAR ) 

I ERR=0 

DO 640 1 = 1 . IP 

IF (PXTPX (I , I ) «GT. 1 *E— 20 ) GO TO 630 
I ERR = 1 

WRITE (MF. 1 030 ) RUN. PASS. I . I .PXTPX { I .1 ) 

PXTPX ( I » I )=1 .E-20 
630 CONTINUE 

640 DDELA ( I ) “SORT (PXTPX ( I. I ) ) 

IF (IERR.EQ.1) CALL PRINTER 
DO 650 I = 1 , IP 
DO 650 J*1 , IP 

PXTPX ( I .J) “PXTPX ( I , J )/( DDELA ( I )*DDELA ( J) ) 

IF (I.EQ.J) PXTPX (I . I > “DDELA ( I ) 

650 CONTINUE 

WRITE (MF.1040) 

DO 660 K = 1 . IP, 8 
KK=K+ 7 

IF (KK.GT.IP) KK»IP 

WRITE (MF.1050) (DPARAM (KKK ) »KKK=K ,KK ) 

DO 660 1 = 1 . IP 

WRITE (MF.I060) (DPARAM ( I ). (PXTPX ( I . J ). J=K ,KK ) ) 

660 CONTINUE 

WRITE (MF.1070) 

670 CONTINUE 

DO 680 1=1 .NPAR 
680 AL ( I ) = AL ( I ) +DELA ( I ) 

IF (LOGIC (6)) CZDE«CMDE*CMCON 
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690 

700 

710 

720 

C 

730 

740 

750 

760 

770 

780 

790 


800 


810 


CONTINUE 

IF ( .NOT. LOGIC ( 1 ) ) GO TO 730 
DO 700 1«1<1A 
DO 7o0 J = I » I A 
WT( I . J ) = S0 ( I .J) 

WT(J, I ) =SD ( I . J) 

IF ( .NOT. LOGIC (2 )) GO TO 720 
DO 710 1*1 • 1A 
DO 710 J=1 . I A 
IF (I.NE.J) WT(I.J)=0. 

CONTINUE 
CONTINUE 
WT = (WT)INV 

CALL MASCNT ( INVWT . WT . DET2 « AA2 ) 
IF ( DET2.NE .0 • ) GO TO 730 
WL ( 19) =.T. 

GO TO 750 
CONTINUE 

IF (LOG I C ( 1 0 ) ) LOGI C (9 ) = «T. 

GO TO 780 
CALL HALT 
CALL DSPLAY 
CALL OPERATE 

IF < .NOT .FSS (151) GO TO 750 
DO 760 1=10.19 
WL (I >=.F. 

GO TO 780 
CALL HALT 
CONTINUE 
WL (20) = . T. 

CALL ENDPLOT 

IF (.NOT.FSSO) ) GO TO 920 

CALL UNLODE 

CALL CLRTABL 

NPASS=1 

NFREQ=1 

N»ITF=0 

NTYPF=0 

TGAIN=TIMF*.5 

TOFF=-l • 

UOFF ■ - U CRTS I /UMAX 

IF ( ,NOT.FSS( 1 ) ) GO TO 810 
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CALL CRTPLOT ( 1 . 3 .NFREO .NT YPE .NR I TE. T . TGA I N .TOFF , L ABT , AX I . AXMAX .0. E 
l.LABAX) E 

CALL CRTPLOT ( 2 . 3 .NFREQ .NT YPE .NR I TE. T .TGA I N . TOFF .L ABT . AY I . AYMAX « 0. E 
l.LABAY) E 

CALL CRTPLOT ( 3 . 3 .NFREQ .NTYPE .NR I TE . T . TGA I N. TOFF ,L ABT • AZ I . AZMAX .0. E 
l.LABAZ) E 

GO TO 840 E 

IF ( .NOT .FSS(7) ) GO TO 820 E 

CALL CRTPLOT ( 1 .4 .NFREQ .NTYPE .NR I TE . T.TGA I N ♦ TOFF .LAST .P. PMAx . 0« .LA E 
IBP) E 

CALL CRTPLOT ( 2 .4 .NFREQ . NTYPE. NR ITE . T . TGA I N . TOFF ,L ABT .R.RMAx . 0. .LA E 
1 BR ) E 

CALL CRTPLOT ( 3 .4 .NFREQ . NTYPE .NR I TE . T . TGA I N . TOFF ,L ABT . V . VMAx .0. .LA E 

1 BV ) E 

CALL CRTPLOT ( 4 . 4 .NFREQ . NTYPE. NR I TE . T . TGA I N. TOFF ,L ABT .PH I .PHMAX .0* E 

I.LABPH) E 

GO TO 840 E 


0904 

0905 

0906 

0907 

0908 

0909 

0910 
091 1 

0912 

0913 

0914 

0915 

0916 

0917 

0918 

0919 

0920 

0921 

0922 

0923 

0924 

0925 

0926 

0927 

0928 

0929 

0930 

0931 

0932 

0933 

0934 

0935 

0936 

0937 

0938 

0939 

0940 

0941 

0942 

0943 

0944 

0945 

0946 

0947 

0948 

0949 

0950 

0951 

0952 

0953 

0954 

0955 

0956 

0957 

0958 

0959 

0960 

0961 

0962 


73 



820 
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IF < .NOT.FSSI8) > GO TO 830 
CALL CRTPLOT <1 «3 .NFREQ , NTYPE »N 
1LABDA) 

CALL CRTPLOT < 2 . 3 .NFREQ .NT YPE .N 
1LABDF ) 

CALL CRTPLOT ( 3 » 3, NFREQ .NTYPE « N 

1 LABOR) 

GO TO 840 
CONTINUE 

CALL CRTPLOT ( 1 .4. NFREQ. NTYPE. Is 
lLASU) 

CALL CRTPLOT < 2 »4 .NFREQ .NTYPE »N 
1BW) 

CALL CRTPLOT < 3 « 4 .NFREQ .NTYPE .K 

2 8Q) 

CALL CRTPLOT < 4 .4 .NFREQ .NT YPE ,N 
1 .LABTH) 

CONTINUE 

IF (nPASS.EQ.O) GO TO 900 
t=~delx 

IF (LO I S I < 1 9 ) ) NOEL a NOP T S 
DO 890 1=1, NOEL 
T=T+DELX 

IF ( .NOT.FSSI 1 ) i GO TO 850 
AX I = AXM ( I ) 

AY I =AYM ( I ) 

AZl=AZM< I ) 

GO TO 880 
CONTINUE 

IF (,NOT.FSS(7) ) GO TO 860 
P=PM ( I ) 

R=RM ( I ) 

V=VMC I ) 

PHI.pHM < I ) 

GO TO 880 

IF (.NOT.FSSVS) >.: GO TO 87C 
DA»FOA ( I ) * 

oe=foe(I) '■ " 

OR = FoR< * ) 

GO TO 880 
U=UM ( I ) : 

W=WM ( I ) 

QaQM(I) 

THE.THM ( I ) 

CALL RITECRT ( .T . , .F. , MAXP AGE ) 

CONTINUE 

CALL ENOPLOT 

IF (FSS ( 8 ) ) GO TO 910 

CALL CLRTABL 

NPASS=0 

NFREO*N01TSPS 

NR ITEM 

ntype»- l 

GO TO 800 

CALL PLAYBAK (91 OS ) 

CALL RITECRT (»T«,#F« < MAXP AGE ) 
GO TO 900 
CALL ENDPLOT 
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CALL DSPLAY 

CALL OPERATE 

IF (FSS ( 1 4 ) > GO TO 790 

IF < .NOT «FSS < 1 3 ) ) GO TO 920 

CALL CLRPLOT 

IF ( IPRINT.NE.4 ) CALL ERASE 


WL <20)=.F. 

GO TO 10 
90004 CONTINUE 

CALL CALPLT (0..0..999) 

CALL ATERM 

90014 CONTINUE 

90015 CONTINUE 
RETURN 

960 FORMAT </6H RUN=.F4.0«7H PASS* .F4 ■ 0 *7H DET 1 « .E 1 4.6.7H DET2=.E1 


970 
980 
990 
1 000 
1010 
1020 


14.6.6H AJP=«E14.6) 

FORMAT (/39H ACTIVE ALGORITHM VARIABLES ARE • . • . A3. 1 0 ( A 1 . A3 > ) 

FORMAT (/19H PACKED MEAN ARRAY) 

FORMAT (11E12.4) 

FORMAT I/32H PACKED NOISE COVARIANCE MATRIX) 

FORMAT I/22H PACKED WEIGHT MATRIX) 

FORMAT </2X.3(5HPARAM.6X»5HVALUE. 1 OX . 6HCHANGE • 5X )/ ( 3 (2XA5 « 2XE1 3.6« 


1030 

1 040 
1050 
1060 
1070 


12XE13.6))) 

FORMAT (/6H RUN=.F6.0,7H PASS= <F6< 0 «8H PXTPX ( • I 2 . 1 H. 

14.6) 

FORMAT </38H MODIFIED PARAMETER COVARIANCE MATRIX) 
FORMAT </8<lOXA5)) 

FORMAT (2XA5«8(E15.6 ) ) 

FORMAT < 1 H I » 

END 


I 2.2H > = .El 


E 1021 
E 1022 
E 1023 
E 1024 
E 1025 
E 1026 
E 1027 
E 1028 
E 1029 
E 1030 
E 1031 
E 1032 
E 1033 
E 1034 
E 1035 
E 1036 
E 1037 
E 1038 
E 1039 
E 1040 
E 1041 
E 1042 
E 1043 
E 1044 
E 1045 
E 1046 
E 1047 
E 1048 
E 1049 
E 1050 


PROGRAM USAGE 


The use of the program is demonstrated by showing the setup of the longitudinal 
equations of motion, run procedure, and output listings for the test case. The test case 
consists of pseudo flight data, which are generated by integrating the longitudinal equations 
of motion for fixed parameter values (called the true values) and then adding random 
number sequences (measurement noise) to the variables. The parameter values are then 
offset to become the starting point of the estimation program. The integration scheme 
(from subroutine IGRATE1) used is second-order Adams-Bashforth, a 1-pass integration 
scheme. (See ref. 1.) 


Test Case Setup 

The longitudinal equations of motion are 


u = -qw - g sin 9 + | £ V 2 s(c x , 0 ) 


W = qu + g cos 9 + | £ V 2 s(c z>0 + C Zc ^a!a + C^S 


J 5 e e 


( 10 ) 

(ID 
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where 


. , H& v2se ( c 


m,o 


-m 




c H£.+ 
^ m q 2V 



( 12 ) 


e = q 


(13) 


®e “ 


I 0.1 sin 2.5t 

|o 


V = \ju 2 + 


= tan-1 w 
a u 


(0 §t £ ir/1.25) 
(t > jt/1.25) 


The longitudinal equations are generated from the equations of motion in appendix B 
by the variable -dimensioning arrays described in appendix D. 

For the active equation variables u, w, q, and 6, the input array INTX is 
INTX = (1, 0, 1, 0, l, 0, 1, 0) 

and hence 

INTV = (1, 3, 5, 7, 0, . . 0) (IV = 4) 

For the active performance index variables u, w, q, and 6, the input array 
INTY is 


INTY = (1, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0) 

and hence 

INTA = (1, 3, 5, 7, 0, . . 0) (IA = IA1 = 4) 
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For the active parameters (see TABLE array in section "Display Arrays"), all the 
INTEG40OD values are 0 except for I = 2, 8, 9, 11, 13, 14, 16, 17, and hence 

INTP = (2, 8, 9, 11, 13, 14, 16, 17, 0,. . ., 0) (IP = 8) 

By putting in these arrays and the FORTRAN variables listed in the subroutine 
DSPLAY arrays, the test case is set up. 

Test Case Run Procedure 

The program deck or the data cell control cards are read into the computer after 
selecting RESET mode and both RELEASE switches, and FSS 3, 4, 9, and 16 true. The 
dynamic check case is then run before the test case or normal use of the parameter esti- 
mation program. 

The step-by-step procedure for running the test case from the control console is 
described as follows: 


Steps 

(1) Change constants, ol and activeness by 

means of DDDU. Depress FSS(3), 

- FSS(4), FSS(5), and FSS(9), and then 
release FSS(4). 

(2) Depress PRINT switch. 

(3) When RESET light comes on, depress 

lower RELEASE switch. 

(4) Depress OPER switch. 

(5) When WL(20) comes on: release FSS(3) 

and FSS(9); depress RESET switch and 
depress FSS(13), then release FSS(13). 

(6) Change to desired standard deviations 

and change IPRINT to 3 by means of the 
the DDDU, then depress PRINT switch. 

(7) When RESET light comes on, depress 

lower RELEASE switch. 


Task 

Program initialization. 

I.C. printout of true parameters. 
Return to level (4,0). 

OPERATE to fill flight-data arrays 
with pseudo data. 

End of pseudo data fill and return to 
RESET mode. 

Adds random numbers to pseudo 
data and prints characteristics. 

Return to level (4,0). 
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Steps 


Task 


(8) Change a to offset values and change 

IPRINT to 1, then depress PRINT 
switch. 

(9) When RESET light comes on, depress 

lower RELEASE switch. 

(10) Depress FSS(10) and OPER switch. 

(11) When WL(20) comes on, release FSS(IO) 

and depress FSS(13). 

(12) Release FSS(13) at convergence. 

(13) When WL(20) comes on: depress RESET 

switch; and depress FSS(13), then 
release. 

(14) Depress PRINT switch. 

(15) When RESET light comes on, depress 

lower RELEASE switch. 


Inputs a offsets and gets I.C. 
printout. 

Return to level (4,0). 

OPERATE to obtain initial weighting. 
OPERATE automatically updating a. 

Stops run. 

Return to RESET mode. 

Routes maximum likelihood printout. 
Return to level (4,0). 


Output Listings 

Four output listings are presented to illustrate the computer printouts. (Note that 
variables in yj are referred to as algorithm variables in the output listings.) 

(1) Random number characteristics where RN and RSD denote the means and 

standard deviations, respectively, used in the program 

(2) Initial condition printout 

(a) Showing the true parameter values used to generate the pseudo test case 

(b) Showing the offset parameter values used as the starting point for the 

estimation procedure 
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(3) Maximum liklihood printout for each iteration (a total of 11) where the modified 
covariance matrix for the parameters is 



Tabulated results of the test cases are presented in reference 1. Figure 2 shows the 
CalComp plot representation of the CRT display of the converged solution and the pseudo 
data, and the control inputs. 

RANDOM NUMBER CHARACTERISTICS 


RUN* 1 NQPTS* 201 


DRM< 1 1 

RN( I ) 

DRSD< I 1 

R SO (II 

0 . 

-7.994999E-15 

5.0C0C0CE-01 

5 .O00Q00E— 01 

0 . 

3.691090E-15 

3. OOOOOCE+OO 

3.C000C0E+C0 

0 . 

-3* 866448E-1 7 

2.000000E-02 

2.COOOOOE-C2 

0 . 

• -5.1 53923E-1 6 

2.000000E-02 

2.000000E-02 
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INITIAL CONDITION PRINTOUT (TRUE PARAMETER VALUES) 
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INITIAL CONDITION PRINTOUT (OFFSET PARAMETER VALUES) 
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CONCLUDING REMARKS 


A computer program has been developed for estimating aircraft stability and con- 
trol parameters from flight test data. The maximum liklihood estimation program has 
been implemented on the Langley real-time simulation system. The control and display 
capabilities of the system allow the analyst to interact with the program. The interactive 
capability is highly desirable, as evident in the reports on the analysis of flight test data. 
Variable dimensioning allows the analyst to activate any part of the nonlinear six-degree- 
of-freedom aircraft mathematical model, select the variables in the performance index 
function, and choose which parameters are to be estimated. Although this report uses a 
particular aircraft example, it is applicable to any dynamic system fitting into the frame- 
work of the program. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., May 9, 1973. 



APPENDIX A 


PROGRAM CONTROL AND DISPLAY CAPABILITIES 

The computer program has been written in FORTRAN IV (75 000 octal locations) 
and run on the RTS system of the Control Data series 6000 digital computer complex. 

The computer program was mechanized into an iterative estimation procedure with manual 
interactive control and graphic display capabilities through the. utilization of the RTS sys- 
tem. Figure 3(a) shows a photograph of the program control station and figure 3(b) shows 
a closeup of the control panel. The components are listed below as they appear (left to 
right) in figure 3: 

Program control station: 

Graphic display unit 

Cathode ray tube (CRT) 

Interactive keyboard 
Time history recorder 
x-y plotter (not used) 

Control console 

White indicator lights (WL) 

Red indicator lights, bottom row (not used) 

Function sense switches (FSS) 

Mode control switches 

Data entry keyboard 

Digital decimal display unit (DDDU) 

Potentiometers (not used) 

Output device (typewriter) 

The CRT displays the flight test maneuver at the start of each iteration. The 
response of the equations of motion as it is computed in the digital program is plotted 
with the flight test maneuver for direct comparison. This display permits quick analysis 
of each flight test case on an iteration to iteration basis, figure 4 shows CalComp plots 
representing three CRT displays; they are part of the dynamic check. 

The analyst investigating the stability and control derivatives of the aircraft has 
direct control of the computer program through the control console. The white indicator 
lights (WL(1) - WL(39)) are used to indicate program status or diagnostics. The diagnos- 
tics are described in the LDISO array of the Display Arrays section. The function sense 
switches (FSS(l) to FSS(16)) are used to select program options (switch depressed results 
in logical true value). The options are described in the LDISI array of the Display 
Arrays section. 
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APPENDIX A - Continued 


The mode control keyboard (shown below) is used to control the running of the RTS 
computer program. 


OPER 

HOLD 

RESET 

TERM 

CHANGE 

SCAN 

RELEASE 




ERASE 


IDLE 

READ 

PRINT 

RELEASE 


Each switch (mode) is briefly described as to its use (mode nominally active when switch 
depressed): 

OPER (OPERATE) - allows normal running of parameter estimation procedure 
(integration of equations of motion and sensitivity equations) 

HOLD - holds estimation procedure at last time point (stops integration) 

RESET - initializes estimation procedure at t = 0 

TERM (TERMINATE) - terminates program at control console and transfers 
control to graphic display unit 

Activation of one mode automatically deactivates the previous one. The following modes 
are temporarily activated by the analyst during the parameter estimation study (normally 
when in RESET or HOLD modes): , 

CHANGE - changes program variable to the new value entered on the 
data entry keyboard and displayed on the DDDU 

SCAN - scans through the display addresses in conjunction with 
subroutine SCANNER 

RELEASE - releases CHANGE and SCAN modes 

ERASE - erases real-time disk file 
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APPENDIX A - Continued 


IDLE - idles the computer (no computations) 

READ - loads read overlay 

PRINT - loads print overlay 

RELEASE - releases the four preceding modes 


The data entry keyboard (shown below) is used to input new values for program 
variables. 


0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

DECIMAL 

POINT 

(-) ' 

TAB 

© 

ERASE 

CR 


The keyboard is used in conjunction with the DDDU (shown below). 


(Address field) 


Al 


a 2 



(Magnitude field) 


± 

Mi 



m 2 

m 3 

m 4 

m 5 


(Exponent field) 



The procedure for changing a floating-point number is as follows: 

(1) Enter address field - Aj A2 A3 

(2) Depress TAB 

(3) Enter magnitude field - ± . M2 M3 M 4 Mg 

(4) Depress TAB 

(5) Enter exponent field - ± Ej E2 

(6) Depress TAB 
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APPENDIX A - Concluded 


(7) Depress (CR) 

(8) Depress CHANGE 

(9) Depress RELEASE 

As the numbers and signs (plus sign assumed) are entered on the keyboard, they are dis- 
played on the DDDU. The DDDU shows the final form of the number entered by the key- 
board. Integers and logical variables are entered in a similar manner but with a different 
format. The switches ERASE and CR are used to erase the data field and character just 
entered, respectively. 

The typewriter is used to type out the new and old values of the program variables. 
The time history recorder plots the variables defined in the DAC array (Display Arrays 
section). The interactive keyboard is used to restart and exit the program from the RTS 
system. A card reader and high-speed printer are located near the program control 
station and are easily accessible to the program operator. 
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APPENDIX B 


EQUATIONS OF MOTION AND ACCELEROMETER EQUATIONS 


The mathematical model is that of a nonlinear six-degree-of-freedom rigid-body 
aircraft, in particular, a V/STOL tilt -wing aircraft. The equations of motion are 


x = F(x,3,6,V,a a ,a a ,|9,fl) 
= [F t , F 2 , . . F 8 ] t 


The state vector is 

X = [xj, X 2 , . . x 8 ] T 

= [u,v,w,p,q,r,0,(^] T 
The parameter vector is 

a = [aj, a 2 a 4 J] T 


u(0), (Cx) a _ 


,,t> 5 e,t’ 



The control deflection vector is 
® = {^a ,6 e’ 6 r] 

The equations of motion in detail are 
u = F 1 (x,a,6,V,a a ) 

T 

= -qw + rv - g sin 6 + a 1 V 2 (c xl + C x2 ) + ~ 


(Bl) 


(B2) 


(B3) 


(B4) 


(B5) 


98 



APPENDIX B - Continued 



-ru + pw + g cos 8 sin cp + ajV 2 (Cyi + C Y2 ) + 


1 - a3C Y 


0 


w =-F 3 (x,S,6,V,Q! a ) 

= -pv + qu + g cos 8 cos (p + ai v2 (c zl 



(B6) 


(B7) 


p = F 4 (x,3,6,V,0,/§) 


= a 6 F4(x,a,6,V,/3,^) + b^g^a, 6,V,/3,/3) 

= a 6 b 2 qr + I xz pq + a 4 V 2 (c a + C i2 ) + a 4 V s 2 (c z3 ) + M x 
baPq ' feV + a 4V 2 (C n i + C n2 ) - a 4 V s 2 (c n3 ) + M z 


L 
+ b 4 


q = F 5 (x;3,8,V,a a ,a a ) 

= a 7 pr + b 4 (r 2 - p 2 ) + agV 2 (c ml + C m2 ) + ^ 

r = F 6 (x,3,5,V,/5,/§) 

= biF^(x,3,6,V,/3,/3) + b 5 Fg (x,3, 


6 = F 7 (x) 

= q cos (p - r sin <p 


(B8) 


(B9) 


(BIO) 


(Bll) 
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APPENDIX B - Continued 


* - Fg(x) 

= p + (q sin 0 + r cos 0) tan 6 (B12) 

where 


a l = \ m S 

a 2 = a l(f) 

a 3 - a i(i) 

a 4 = | pSb 

a 5 = a 4 (|) 

l z 

6 ¥z • ! xz 2 

l Z - !x 
7_ I Y 

a 8 = l^ S6 

a 9 = a s(f) 

bi = IXZ 2 

Wz ' J X2 2 

b 2 = Iy - Iz 

b 3 = ! x - Iy 

. J xz 

b4= T7 

h IX 

5 ¥z ■ J xz 2 



and 


V S = Vgg + V 


C X1 


+ Cx «.( 

<*a - «a,t) 

C X2 

- r 

Cx q 2V 


C Y1 

' Hvr 

/ Cy * 

0 3 ' M +C Y 6 ( 6 r - 6 r,t) 
r 

° y 2 

= ^(° Y P P 

+ CY r r ) 

C Z1 

= ( CZ ) 

“a,t. 6 e,t 


“a - «a,t) + C Z 5e («e * 6 e,t) 

CZ2 

= C Z si 

z q 2 v 


c u = 

*fo) 

V Va.tA.t 

* c < s ( e 

" ^t) + C i 6r ( 6 r ' 6 r,t) 

C Z2 = 

■nr(V* 

%P + C ir r ) 

C ^a 

= cos i w - 

C ng sini w 
u a 

Ci3 = 

’ C ‘8 a ( a »' 

6 a,t) 

^ml 


+ Cm a a 

(“a * a a,t) + c ra^( 5 e * 6 e,t) 

Cm2 

' W ( Cm « a “ a + Cm q q ) 

C„1 = 

z (c } 

V n Va,t’ 6 r,i 


! • ^t) + C n 6 ( 6 r - 6 r,t) 
u r 

C„2 

■w( c y* 

' c n p P + c n r rj 

c* 

% 

= % a sini w + 

C n cosi w 
&a 

c n3 : 


6 a,t) 
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APPENDIX B - Continued 


The auxiliary equations are 


"'v 


• _ q sin <f) + r cos <p 
cos 6 


V = \/u2 + + w2 


or a = tan-1 2 

d. U 


• ~ w 
®a~ u 


P = sin-1 v 


P V 


17 _ 8 _ 2 2 p. 

V SS ~ D E C T 


The trim conditions are 


ALPHAT or 


a 


a,t 


tan 


-1 w (°) 
u(0) 


fit = ° 


6 a t = DABIAS 
6 e ,t = DEBIAS 
6 t = DRBIAS 


(B13) 






APPENDIX B - Continued 


The thrust and moment equations are known inputs (flight test data and constants) 
to the equations of motion — .. , „ . ' , - . - — . - » 

T x = 4 PT?e 2 °E 4 C T cos % 

T y = 0 

T Z = Cy sin i w - pr/rp D,p 4 

_ > (B14) 

M x = 2pr? E 2 D E 4 fCq> B 6B sin 

My — “4p7J E 2 rjj Cy “ P?](p Drp4 IfppCrp^ 

M z = 2prj E 2 D E 4 i Cy B 6B cos i w 


where 


c T - c T o + c Tb ab 


Crp — Crp + C r P B 

t t t t g t /3 t p 


1232 


“ “ 60 " *E 


2400 


^=- 60 "^ 


T 


The accelerometer measurements and equations were included in the parameter 
estimation algorithm to improve the extraction process. They are used with or can 
replace the linear velocities u, v, and w. The accelerometer equations were trans- 
formed to the instrument location from the center of gravity (ref. 6). 
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APPENDIX B - Concluded 


a I = 


a X,I 

a Y,I 

a Z,I 


u + qw - rv + g sin 6 


-(q 2 + r 2 )x x + (pq - r)y x + (pr + q)z x 

v + ru - pw - g cos 6 sin <p 

1 

+ g 

(pq + r)x y - (p 2 + r 2 )y y + (qr - p)z y 

w + pv - qu - g cos 6 cos (p 


(pr - q)x z + (qr + p)y z - (p 2 + q 2 )z z 


(B15) 


The accelerometer equations need only to be evaluated and not integrated. 
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APPENDIX C 


SENSITIVITY EQUATIONS AND ACCELEROMETER 
SENSITIVITY COEFFICIENTS 


The sensitivity equations for the method of quasilinearization are presented in 
detail for the equations of motion presented in appendix B. 

The sensitivity equations are 


d_ 

dt 




where 


8 _ 


_ 3 


,k=l 


9xiA a? V 3V f dx k\ 8F V da af dx k ^ 

da a Lj 9xv \9a; ; , 


J V ±F_r^k\ dF y 

\ L dxAdaJ dV L 


6 ^{daj + 8a a 3x k 


3F 9( *a| 

(9x^ 

| + ®I 

“ 

3/3 l 

('9x 2 '\ 

3 

1 + m y .sv 1 

; 8 x k'i 

9a a 9 xj ' 

V 8o V 

1 9/3 

3x 2 ' 

\ 3Q? iy 

1 9 V L 8 xJ 
k=l K 

\ 8a v 


• 3 /a _ si 

9F 9£ V 3V/ 9x kU 
8/3 8V 9x k ^aaijj 


8F 9 <*a VFJ±(*2> 


9a a 8x 3 \»ai j 9/5 9x 2 


8F 

9Q!j 




8 F 9 < *a / 9 *3\ + 3F 80 / 9 x 2 N 
8 a a 8 x 3 \ da iJ 3/3 9x 2 \ 9 Q! i/ 


8 F 

9a< 


G'(‘> ■ [«ik<‘>l 


(j,k= 1 , 2 , . . 8 ) 


The functions F 2 and F 3 do not contain a a or /3. Thus, 

8 


d f dx\ _ 


dt ySa'i 


i/ \ 9 a i/ 8a a 9x 3 3k \9o ij dp 8x 2 ^ g 2k \p a J 


' 3F + _8F_ ^a ^3 + 3F _£0_ 9f V 
,90!! + 8a a 8x 3 9a. + d p 0x 2 9aj 


(Cl) 


= G(t) 


9x 

9a 
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(i = 1, 2, . . 40) 


(C2) 



APPENDIX C - Continued 


9x _ fsu 8y _9w j^P_ _9q_ _9r_ _90_ 

9a^ ~ [aotj’ Sa^’ ao^’ da 9aq’ Sa^ da ^ J 

G(t) = [g jk (t)‘ 

F(o-i) = F^), F 2 { a i)> . . Fg^i) T 


(j,k = 1, 2, . . 8) 


(1) Sensitivity equations derived from u equation: 


'iH_y = y 

sdai ' k=l 


ax k \ 9FJ 


g lktei + ^- 


gll * aju(2C xl + C x2 ) - ajC x w 

“a 

gl2 = r + a x v(2C xl + C x2 ) 

g 13 * -q + a x w(2C xl + C x2 ) + a^ u 


gl5 = - w + a 2 VC X 


g 17 = -g cos 8 


^18 = 0 



APPENDIX C - Continued 


and 


3F 


9 «2 


I = ajV 2 


9 Fi o / \ 

— = a lV 2 ( aa - V ) 


aFj 

00!, 


= a 2 Vq 


(2) Sensitivity equations derived from v equation: 
8 

(»L\. 

dt' 


d./jv_\ _ y fok\ 8F 2 

dtlaa. J L g 2kl 9 q;. J 9 a . 

\ V k=l ' iy 1 


where 


821 “ 


-r + a 1 u^2C yl + C y2 - c y^) 

1 " a 3 C Yo 


S22 


a l v ( 2C Yl + Cy2 - C YpP) + a l VC Y^ 


1 - a 3Cy 


0 


g 23 


g 24 


p + ajW^2Cyi + Cy 2 - Cy^/3\ 


1 - a 3 c Y^ 


W + agVCy 


1 - a 3Cy 


3 


g 25 


= 0 


-u + a 3 VCy 


g 26 “ 


1 - a 3 C Y 


0 
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APPENDIX C - Continued 
-g sin 6 sin 0 

««- — i SCyJ - 


_ g cos 9 cos 0 
g 28 l _ a 3 C y ^ 


and 


9F 2 ajV"2 
ao? l9 1 " a 3 c Y 


/3 


9F 2 

aa 20 1 ' a 3 c Y^ 


9F 2 a 3 v 
8q? 21 1 ~ a 3 c Y^ 


9F2 a 3 Vp 
9«22 ~ 1 - a 3C Y ^ 


9F 2 a 3 Vr 

9a 23 1 - a 3 C Y ^ 

3F 2 ai v 8 (e r - » r , t ) 
8 «24 1 " a 3 C Y / § 


(3) Sensitivity equations derived from w equation: 



+ 


9F3 

Ba^ 


where 


(C5) 


g 31 « q + a lU (2C zl + C 22 ) - ajCg w 
g 32 = -p + aiv(2C zx + C Z2 ) 
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APPENDIX C - Continued 


g 33 * a x w(2C zl + C Z2 ) + a lCZo! a U 
g 34 = " v 


g 35 = u + a 2 VC Zq 
g 36 = 0 


g 37 = -g sin 6 cos <j> 


g 3g = -g cos 6 sin 0 


and 




da 


8 


£ ■ a l y2 ( a a - “a,t) 


9Fq 


9F 3 

da ll 


= ai v2(( 


6 e “ 6 e,t 


(4) Sensitivity equations derived from p equation: 


dt \Sa 


d / 8p ) _ 


:a 6 


k ?iV 4k+ 


4 3/3 „ V 3x k 


g 4k + 8v g 2kH8 ai ) + V 8a i + 80 9v 3 «i 


9F 4 8F 4 3(3 ®£2 


+ b, 


V /„ 8F 6 8£ V9x k \ /SFg 8F 6 8g 8F 2 

Z ^ g 6k + dp 3v e2k )\9ati) + \55f + a/ 3 8 y Sari 


U 

■I 

k=l 


g 4k -6^1 g 2k ^ 3 g 2k) (g|) + V 4 («i) ♦ b l ? 6 (“,) 


O 

'I 

k=l 




(C6) 
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APPENDIX C - Continued 

where ; 

641 ' »4 u ( 2C £X + C 12 * 2C iS S . Ci(j(3 . C; • |) 

S42 * v( 2C n ♦ Ci2 + 2C i3 S + a 4 VC i(3 

S« * v( 2C £l + C £2 + 2c £3 T - V ' fe): 

§44 = *XZ q + a 5 VC ip 
g 45 = b 2 r + *XZ P 
g 46 = b 2 q + a 5 VC *r 
g 47 = g 48 = ° 

and 

g 61 " a 4 u ( 2C nl - C n2 - 2C n3 S - - C n . ^ 

g 62 * a 4 v (^ C nl + C n2 " 2C n3 If ' ~ C n^ + ®4 vc n^ 

g 63 ~ a 4 w (^ C ni + c n2 ' 2C n3 \T " C n^ “ c n^ 

g 64 = b 3 q + a 5 VC n p 

g 65 = b 3 p ' J XZ r 


109 



APPENDIX C - Continued 


§66 = ~ha& + a 5 VC n r 

rr ft 

g 67 = g 68 " 0 

and 


sf; 

9F 6 _ 

a 4 V 2 

9 <*26 

00-34 


9F^ 

9q! 27 

9F' _ 
8a 35 

a 4 V 2 (f?-flt) 

3f; 

3Fr 


9q! 28 

D __ 

8< *36 

a 5 V^ 

9F’ 

9F’ 


8 «29 

0 ___ 
90-37 

a 5 Vp 

8F^ 

0F’ _ 

a.Vr 

da 30 

00-38 

5 

9F^ _ 
9o! 31 ~ 

8F 6 _ 
8q; 39 

a 4^ 2 ( 5 r " 5 r,t) 

da 32 

8F e ; 

30-40 

= a 4 V s 2 cos i w (s a - 

9F 4 _ 

9 «40 

9F 6 _ 
00-32 

-a 4 V s 2 sin i^ - 

9F 4 3/3 
0/3 0V 

= a 5 % 


S -Ikm 

= aeCn; 


dp 0v 

*' P 
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APPENDIX C - Continued 
(5) Sensitivity equations derived from q equation: 

v ' 

_d / V\ _ f L , \ fa k\ ffa 8F 3 \ 

dtUaJ L y 5k 9( ^ a g 3kJ^9 ai j \s ai 8d a ;. &w 9 «i/ 

' ' k — 1 

- I S 5k (||) + F 5 (a i) < C ’> 

k=l ^ x/ 

where 

g 51 ~ a 8 u ( 2C ml + C m2) " a 8 C m Q! w " a 9 C m^ “a 

a. cL ' 

g 52 = a 8 v ( 2C ml + c m2) 

g 53 * a 8 w ( 2C ml + C m 2 ) + a 8 C ™c* a u 
g 54 = a 7 r “ 
g 55 = a 9 VC mq 
g 56 = a 7 p + 2b 4 r 
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APPENDIX C - Continued 


9F 


a — = a q Vd„ 

9q! 15 y a 


9Fc 


9a 


16 


= a g Vq 


9Fr 


9a 


17 


= a 8 V 2 (6 e - 5 e 


.*) 


3F 5 8 °a 
9a a 9w 


a 9^m 



(6) Sensitivity equations derived from r equation: 


d / 9r Y b 
dtUa^y 1 


O 

I 

k=l 


g - + .£n*L V 9x k 

g 4k + a,V g 2k 


r 9F4 9F4 a/3 9F 2 ' 


9/3 9v 2k y\9qi/ \ 8Q! i 9/3 9v 9a iy 


+ b c 


I 

k=l 


9Fl 


6 W 


ax.A ,»FJ 8 p »f 2 ' 


S 6k + 9/3 9v g2k /\ 9 “i/ + \904 + 9/3 9v 9 «i 


1 ( E 6k + b l if if g 2k + b 5 f g 2b)(^) + b l n(“i) + b 5 ? K a i) 



(C8) 


where all the terms have been defined in the derivation of the sensitivity equations for 
p equation. 
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APPENDIX C - Continued 


(7) Sensitivity equations derived from 6 equation: 



where 

g 71 = g 72 = g 73 = g 74 = 0 
g 75 = cos <t> 

g 76 = " sin 0 
g 77 = 0 

g 78 = cos e 

(8) Sensitivity equations derived from $ equation: 



where 

g 81 = e 82 = %3 = 0 

hi * 1 

ggg = sin (p tan 6 
ggg = cos <p tan 9 


(C9) 


(CIO) 
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APPENDIX C - Concluded 


. _ it 

^87 cos q 

ggg = 9 tan e 


The accelerometer sensitivity coefficients were derived in terms of the sensitivity 
equations and coefficients for the equations of motion, and need only to be evaluated and 
not integrated. 

(1) Sensitivity equations derived from a^ j equation: 


aa- 


X,I _ 1 
a«i “g 


plj + q fe) + («* + z x r )(S:) + (~ 2x x q + y X P + W )(S: 


(' 2x x r + z xP - v )(£-) + g cos 8 £) + fe) + z xfe) - Txfe 


(Cll) 


(2) Sensitivity equations derived from ay j equation: 


3a 


M-i 

a^i " g 


r (j^J - p (§) + (- w + x yP - 2 yyp)fe) + ( x yP + z y r )ffir) + ( u : 2y Y r + z **)te 


+ g sin 0 sin ' g cos e cos + (4) ■ Z y(5[) + X Y 


(C12) 


(3) Sensitivity equations derived from ag j equation: 


Bar 


Z,I _ 1 


- q fe) + p fe) + ( v + x z r - 2z z p )&) + (- u + y z r - 2e z«)S) + ( x z p + y z ^)(Sr 


+ g sin 0 cos 0 g) + g cos 0 sin ^ 


(C13) 
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APPENDIX D 


VARIABLE DIMENSIONING 

The flight test runs do not necessitate the use of all the equation variables (x), 
variables in the performance index function (y), and parameters (a) for specific 
cases. These cases involve only a specific part of the program, as with an excitation of 
only the longitudinal motion of the aircraft. Variable dimensioning of the estimation 
procedure furnishes the analyst with the means of altering the program to meet the spe- 
cific needs of each flight test run; that is, the mathematical model of the aircraft dynamic 
response, the variables to be compared with flight test data, and the parameters to be 
estimated. 

Variable dimensioning of the estimation procedure is accomplished by using the 
three input arrays INTX(8), INTY(U), and INTEG, 40 (INTEG(I), 1 = 1,2,. . ., 40); the use 
of these arrays dimension x, y, and a, respectively, in the program. The elements of 
each array are entered as integers 1 or 0 to indicate whether the variables or parameters 
are active or inactive, respectively. 

The input array INTX specifies the activeness for each equation variable in the 
equations of motion and the sensitivity equations; inactive variables are treated as con- 
stants. From the input array INTX 

INTX = (1, 0, 1, . . ., 1, 0) 
k l k 2> • • ■> k IV 

” ' 

/ 

8 

IV- I INTX(K) 

K=1 

where kj, k2, . . ., kjy are element locations, INTV is generated, 

INTV = (k^, k2, . . ., kjy, 0, . . ., 0) 

which is a sequence of integers denoting the active equation variables in x. 

The input array INTY specifies the . activeness for each variable in y; inactive 
variables are ignored. From the input array INTY . 
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APPENDIX D - Continued 


INTY = (1, 0, 1, . , ., 1, 0) 

h h ’ • • -^ia 

n 

IA = ^ INTY(J) 

J=1 ' 

INTA is generated, 

INTA = (jj, j 2 , . . j IA , 0, . . o) 

which is a sequence of integers denoting the active variables in y. The integer IA1, 
where 

8 

IA1 = ^ INTY(J) 

J =i 

denotes the number of active variables in x that are active variables in y. 

The input array INTEG 4 g specifies the activeness of each parameter in the estima- 
tion procedure; inactive parameters are treated as constants. From the input array 

integ 40 

integ 40 = (0, . . 0 , 1, 0 , 1, . . 1, 0 , . . ., 0 ) 

il i 2 , . . i IP 


40 

IP = ^ INTEG(I) 

1=1 

INTP is generated, 


INTP = (ii, i 2 ,. . ., ijp, 0, . . o) 

which is a sequence of integers denoting the active parameters in a. The program is 
dimensioned for IP ^ 30. 

The resulting arrays (INTV, INTA, INTP) and numbers (IV, IA, IA1, IP) are used in 
FORTRAN DO LOOP and matrix operations. 
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APPENDIX D - Continued 
The equations of motion are reduced to the form 

*eq = Feq(x,a ( 6) = (x^, x^, . . ., 


(Dl) 


by multiplying each of the original equations of motion (eq. (la)) by their respective ele- 
ment of INTX; that is, inactive equation variables have their derivatives set to zero. 

The original sensitivity equations (eq. (6)) are reduced by using the arrays INTP 
and INTV as indices in the FORTRAN DO LOOP to the form 


d / 9 x INTV(K , ) N \ _ y 9 *WtV(K') / 9x INTV(K) \ f 8 ^INTV(K t ) 
dt \ 9a INTP(I) / 9 x INTV(K) \ 9q! INTP(I) / 9q! INTP(I) 


(K' = 1, 2, . . ., IV; I = 1, 2, . . ., IP) (D2) 

The accelerometer equations and sensitivity coefficients are handled in a similar manner 
to generate aj and ^j/^intP^)’ res P ective ly- 

The original parameter change equations (eq. (7)) reduce to 



N 

-1 

N 

A«j = 

1 R J 1 A j(‘i) 


^ A j(ti) Rj 1 Vjih) 


_i=l 


_i=l 


where 


A j(*i) 




3 o;. 

X 1 




Vjdi) - yj Oi) ■ yj(ti) 


x JVi) 

i 

—0 

*>> 


yj = 






(D3) 
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APPENDIX D - Concluded 


The vector axgq^a'jjjrppg) is searched at each t^ for the sensitivity coefficients 
associated with the active algorithm variables by using the array INTA. The selected 
sensitivity coefficients are then packed into the vector 9 ^jy/ 9a iNTP(I) * n orc * er to form 
Aj(t). If IA1 = IV, no search is made, since this condition implies that x e q = xj. The 
INTA array also selects ffj (tj) from the original vector rj’(ti). 

The covariance matrix for the parameters becomes 


r n 

2 

i=l 


Aj(ti) RJ 1 Aj(tj) 


-1 


o o a 
Pa h a h a h “il 


o o o 
“il“i 2 a i2 

o2 

i2 


p oo 

a il a ilP "ii "ijp 


p oo 

"i2"iip "i2 a ilP 


o a o p cr cr 

ai ip“il "hp “il “iip“i2 “hp i2 


(D4) 


UP 


The covariance matrix for the measurement noise becomes 


IN 

R°(N) # Estimate of Rj = j- £ %(h) 


i=l 




\\ 






°Vh Vh 


Via 

\\a 


(D5) 


J IA 


and the performance index to be minimized becomes 


j n( s J) = KH 


(D6) 
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APPENDIX E 


CALCOMP PLOT OPTION 

CalComp plots similar to the CRT displays can be obtained in overlay level (2,0). 
The main difference being that the (+) symbol is used to represent flight data points. 

Other differences are options to plot flight data only (FSS(2) true), and a variable (NPLOT) 
to alter the density of flight data points to be plotted. Examples of CalComp plots are 
shown in figures 2 and 4. 

The CalComp plot option is entered by setting IPRINT to 4 by means of the DDDU 
before exiting the CRT display loop. (Replotting of the CRT display just before entering 
the CalComp plot option causes erasure of the real-time disk file, and thereby prevents 
CalComp plots of calculated data.) After exiting the CRT loop, set FSS as shown in the 
CRT discussion of overlay level (4,0), and depress PRINT switch. When processing of the 
selected plot is completed, WL(2) will come on to signal the need for operator action. 
Additional plots can be obtained by setting the appropriate FSS before depressing FSS(14) 
and then releasing. Exiting the CalComp loop is accomplished by depressing FSS(13). 

The flow chart of the CalComp plot option follows. 
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APPENDIX E - Concluded 
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(a) Longitudinal motion. 

Figure 2.- CalComp plot of pseudo flight and calculated (converged solution) 
longitudinal motion, and control inputs. 
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Figure 3.- Operational control features. 
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(a) Longitudinal motion. 

Figure 4.- CalComp plot of pseudo and calculated 
(nonconverged solution) motion. 
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(b) Lateral motion. 
Figure 4.- Continued 
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